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Introduction 



In many fields of applied mathematics, one is led to study dissipa- 
tive evolution equations involving (i) a degenerate dissipative operator, 
and (ii) a conservative operator presenting certain symmetry proper- 
ties, such that the combination of both operators implies convergence 
to a uniquely determined equilibrium state. Typically, the dissipative 
part is not coercive, in the sense that it does not admit a spectral gap; 
instead, it may possess a huge kernel, which is not stable under the 
action of the conservative part. This situation is very similar to prob- 
lems encountered in the theory of hypoellipticity, in which the object 
of study is not convergence to equilibrium, but regularity. By analogy, 
I shall use the word hypocoercivity, suggested to me by Thierry Gal- 
lay, to describe this phenomenon. This vocable will be used somewhat 
loosely in general, and in a more precise sense when occasion arises. 

Once the existence and uniqueness of a steady state has been estab- 
lished (for instance by direct computation, or via an abstract theorem 
such as Perron- Frobenius), there are plenty of soft tools to prove con- 
vergence to this steady state. It is much more tricky and much more 
instructive to find estimates about rates of convergence, and this is the 
question which will be addressed here. 

Both hypoellipticity and hypocoercivity often occur together in 
the study of linear diffusion generators satisfying Hormander's bracket 
condition. It is for such equations that theorems of exponentially 
fast convergence to equilibrium were first established via probabilistic 
tools |36|,I47|, I40|,l41j . taking their roots in the Meyn-Tweedie theory of 
the asymptotic behavior of Markov chains. Some of these studies were 
motivated by the study of finite-dimensional approximations of ran- 
domly forced two-dimensional Navier-Stokes equations |17|, I36|, I35j : 
since then, the theory has been developed to the extent that it can 
deal with truly infinite-dimensional systems |28j . In all these works, 
exponential convergence is established, but there are no quantitative 
estimates of the rate. Moreover, these methods usually try to capture 
information about path behavior, which may be useful in a probabilistic 
perspective, but is more than what we need. 
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Analytical approaches can be expected to provide more precise re- 
sults; they have been considered in at least three (quite different, and 
complementary) settings: 

- For nonlinear equations possessing a distinguished Lyapunov 
functional (entropy, typically), robust methods, based on functional 
inequalities, time- derivative estimates and interpolation, have been de- 
veloped to establish convergence estimates in 0{t~°°), i.e. faster than 
any inverse power of time. These methods have been applied to the 
linear (!) Fokker-Planck equation |13j . the Boltzmann equation |15j . 
and some variants arising in the context of kinetic theory |Hl I21j . So 
far, they rely crucially on strong regularity a priori estimates. 

- For linear hypoelliptic equations enjoying some structural 
properties, more specific methods have been developed to prove (ide- 
ally) exponential convergence to equilibrium with explicit bounds on 
the rate. Up to now, this approach has been mainly developed by 
Herau and Nier |32j . Eckmann and Hairer |19j . Helffer and Nier |30j . 
for second-order differential operators in Hormander's form (a sum of 
squares of derivations, plus a derivation). It uses pseudo-differential 
operators, and a bit of functional calculus; it can be seen as an ex- 
tension of Kohn's celebrated method for the study of hypoellipticity 
of Hormander operators. In fact, the above-mentioned works estab- 
lish hypoellipticity at the same time as hypocoercivity, by considering 
functional spaces with polynomial weights in both Fourier space and 
physical space. After a delicate spectral analysis, they localize the spec- 
trum inside a cusp-like region of the complex plane, and then deduce 
the exponential convergence to equilibrium. Again, in some sense these 
methods capture more than needed, since they provide information on 
the whole spectrum. 

- Finally, Yan Guo recently developed a new method |26j , which he 
later pushed forward with Strain |46|, 12 7j , to get rates of convergence 
for nonlinear kinetic equations in a close-to-equilibrium regime. Al- 
though the method is linear in essence, it is based on robust functional 
inequalities such as interpolation or Poincare inequalities; so it is in 
some sense intermediate between the two previously described lines of 
research. 

The goal of this memoir is to start a systematic study of hypocoer- 
civity in its own right. The basic problem considered here consists 
in identifying general structures in which the interplay between a "con- 
servative" part and a "degenerate dissipative" part lead to convergence 
to equihbrium. 
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With respect to the above-mentioned works, the novelty of the ap- 
proach explored here resides in its abstract nature and its simplicity. 
In particular, I wish to convey the following two messages: 

1. Hypocoercivity is related to, but distinct from hypoellipticity, 
and in many situations can be established quantitatively independently 
of regularity issues, or after regularity issues have been settled. 

2. There are some general and simple techniques, based on very 
elementary but powerful algebraic tricks, by which one can often reduce 
a mysterious hypocoercive situation to a much more standard coercive 
one. 

There are three parts in this memoir: 

Part I focuses on the particular case of operators which (as in |32|, 
II 9|, I30p can be written in "Hormander form" A*A + B, where A and 
B are possibly unbounded operators on a given Hilbert space. These 
results have been applied to several models, such as the kinetic Fokker- 
Planck equation, the linearized Landau-Lifschitz-Gilbert -Maxwell model 
in micromagnetism jH], and a model problem for the stability of the 
Oseen vortices |23j . 

Part II, by far the shortest, remains at a linear level, but considers 
operators which cannot necessarily be written in the form A* A + B, 
at least for "tractable" operators A and B. In this part I shall give 
an abstract version of a powerful hypocoercivity theorem recently es- 
tablished by Mouhot and Neumann |38j . explain why we cannot be 
content with this theorem, and give some suggestions for research in 
this direction. 

In Part III I shall consider fully nonlinear equations, in a scale 
of Sobolev-type spaces, in presence of a "good" Lyapunov functional. 
In this setting I shall obtain results that can apply to a variety of 
nonlinear models, conditionally to smoothness bounds. In particular 
I shall simplify the proof of the main results in |15j . 

Though these three settings are quite different, and far from being 
unified, there is a unity in the methods that will be used: construct 
a Lyapunov functional by adding carefully chosen lower-order 
terms to the "natural" Lyapunov functional. This simple idea 
will turn out to be quite powerful. 

The method will be presented in a rather systematic and abstract 
way. There are several motivations for this choice of presentation. 
First, the methods are general enough and can be applied in various 
contexts. Also, this presentation may be pedagogically relevant, by 
emphasizing the most important features of the problem. Last but not 
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least, most of the time I really had to set the problems in abstract 
terms, to figure out a way of attacking it. 

No attempt will be made here for a qualitative study of the approach 
to equilibrium, but I believe this is a very rich topic, that should be 
addressed in detail in the future. One of the main outcomes of my 
work with Laurent Desvillettes |15j was the prediction that solutions 
of the Boltzmann equation, while approaching equilibrium, would os- 
cillate between "close to hydrodynamic" and "close to homogeneous" 
states. To some extent, this guess was in contradiction with a com- 
monly accepted idea according to which the large-time behavior should 
be dominated by the hydrodynamic regime; nevertheless these oscilla- 
tions have been spectacularly confirmed in numerical simulations by 
Francis Filbet. Further developments can be found in |22j : the re- 
sults obtained by numerical simulations are so neat that they demand 
a precise explanation. 

Research in the area of hypocoercivity is currently developing fast 
thanks to the efforts of several other researchers such as Thierry Gallay, 
Frederic Herau, Clement Mouhot, and others. I expect that further 
important results will soon be available thanks to their efforts, and 
that this memoir really can be considered as a starting point of a much 
more developed theory. 



Part I 

A'A + B 



In this part I shall study the convergence to equilibrium for degen- 
erate linear diffusion equations where the diffusion operator takes the 
abstract form A*A + B, B* = -B. 

The main abstract theorem makes crucial use of commutators, 
in the style of Hormander's hypoellipticity theorem. In its simplest 
form, it reduces the problem of convergence to equilibrium for the non- 
symmetric, non-coercive operator A* A + B, to that of the symmetric, 
possibly coercive operator y4*A+ [A, B]*[A, B]. If the latter operator is 
not coercive, then one may consider iterated commutators [[A, 5], 5], 
[[[A, B],B],B], etc. 

One of the first main results f Theorem I24j) can be informally stated 
as follows: Let A = {Ai, . . . ,Am), B* = -B, and L = A* A + B 
be linear operators on a Hilbert space Ti. Define iterated commuta- 
tors Cj and remainders Rj (1 < j < Nc) by the identities Cq = 
A, [C„ B] = Q+i + R,+i U < N,), Cn^+i = 0. // Ef=o is 
coercive, and the operators [A^Ck\, [A*,Ck\, Rk satisfy certain bounds, 
then 1 1 e"*^ 1 1 = 0{e~^^), where the "Sobolev" space Ti} is defined 
by the Hilbert norm = \\h\\'^ + ^ \\Cjh\\^ . 

The key ingredient in the proof is the construction of an auxiliary 
Hilbert norm, which is equivalent to the li} Hilbert norm, but has 
additional "mixed terms" of the form {Cjh,Cj+ih). 

Applied to the kinetic Fokker-Planck equation, these theorems will 
yield results of convergence to equilibrium that are both more general 
and more precise than previously known estimates. 

After this "abstract" framework, a "concrete" L log L framework 
will be considered, leading to results of convergence for very general 
data (say finite measures). 

My reflexion on this subject started during the preparation of my 
Cours Peccot at the College de France (Paris), in June 2003, and has 
crucially benefited from interactions with many people. The first draft 
of the proof of Theorem occurred to me while I was struggling to 
understand the results of Frederic Herau and Francis Nier |32j about 
kinetic Fokker-Planck equations. The construction of the anisotropic 
Sobolev norm was partly inspired by the reading of papers by Yan 
Guo j26j and Denis Talay |47j : although their results and techniques 
are quite different from the ones in the present paper, they were the 
first to draw my attention to the interest of introducing mixed terms 
such as Vt,/ ■ Vx/. Denis also showed me a useful trick for getting 
long-time estimates on the moments of certain hypoelliptic diffusion 
equations, which is based on the construction of an adequate quadratic 
form. 
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Apart from the above-mentioned people, I was lucky enough to have 
fruitful discussions on the subject with Bernard Helffer, Laurent Desvil- 
lettes, Luc Rey-Bellet, Jean-Pierre Eckmann, Martin Hairer, Clement 
Mouhot, Stefano OUa and Piere-Louis Lions, as well as with Christian 
Schmeiser and Denis Serre, who both suggested a relation between my 
results and Kawashima's condition in the theory of hyperbolic systems 
of conservation laws. 
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1. Notation 

1.1. Basic notation. Let 7i be a separable (real or complex) 
Hilbert space, to be thought of as L'^{fi), where /i is some equilib- 
rium measure; 7i is endowed with a norm || ■ || coming from a scalar 
(or Hermitian) product (■, ■). 

Let V be a finite-dimensional Hilbert space (say M"* or C™, depend- 
ing on whether 7i is a real or complex Hilbert space). Typically, V will 
be the space of those variables on which a certain diffusion operator 
acts. The assumption of finite dimension covers all cases that will be 
considered in applications, but it is not essential. 

Let A : — s> ® V ^ TY™ be an unbounded operator with domain 
D{A), and let B : H H he an unbounded antisymmetric operator 
with domain D{B): 

V/i, h' e D{B), {Bh, h') = -{h, Bh'). 

I shall assume that there is a dense topological vector space S m Ti 
such that S C D{A) fl D{B) and A (resp. B) continuously sends S 
into iS ® V (resp. S); this assumption is here only to guarantee that 
all the computations that will be performed (involving a finite number 
of applications of A, A* and B) are authorized. As a typical example, 
S would be the Schwartz space 5(M^) of functions f : ^ R 
whose derivatives of arbitrary order decrease at infinity faster than all 
inverse polynomials; but it might be a much larger space in case of 
need. 

For a given linear operator S, I shall denote by II^H its operator 
norm: 

\\S\\ = sup ^\ J = sup {Sh,h'). 

If there is need to emphasize that S is considered as a linear operator 
between two spaces Hi and H2, the symbol \\S\\ may be replaced by 

\\S\\ni^H2- 

The norm of an array of operators {Ai, . . . , A^) is defined as 
ll^iP' norm of a matrix-valued operator (Ajk) by a/X] ll^jfcP; 

etc. 

The identity operator X X, viewed as a linear mapping, will 
always be denoted by /, whatever its domain. Often a multiplication 
operator (mapping a function / to fm, where m is a fixed function) 
will be identified with the multiplicator m itself. 

Throughout the text, the real part will be denoted by 3?. 

1.2. Commutators. In the sequel, commutators involving A and 
B will play a crucial role. Since A takes its values in 7i®V and B is only 
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defined in 7Y, some notational convention should first be made precise, 
since [A, 5], for instance, does not a priori make sense. I shall resolve 
this issue by just tensor izing with the identity: -B] = AB — {B^I)A 
is an unbounded operator Ti — > 7i^. In a more pedestrian writing, 
[A, B] is the row of operators ([^i, B], . . . , [Am, B]). Then A"^ stands 
for the matrix of operators {AjAk)j,k, [A, [A,B]] for {[Aj, [Ak, B]])j^k, 
etc. One should be careful about matrix operations made on the com- 
ponents: For instance, [A, A*] stands for {[Aj, Al]j^k), which is an oper- 
ator n^n(^V(^V, while while [A*, A] stands for J2j[A*, Aj], which 
is an operator — > 7Y. Also note that [A, A] stands for the array 
{[Aj, Ak])j,k, and is therefore not necessary equal to 0. When there is 
a risk of confusion, I shall make the notation more explicit. 

1.3. Relative boundedness. Let S and T be two unbounded 
linear operators on a Hilbert space Ti, and let a > 0; then the operator 
S is said to be a-bounded relatively to T if D{T) C D{S), and 

WheD{S), < a||T/i||; 

or equivalently, S*S < aT*T. If S is a-bounded with respect to T for 
some a > 0, then 5* is said to be bounded relatively to T. This will be 
sometimes abbreviated into 

S 4T. 

Note that S and T need not take values in the same space. Of course, 
boundedness relative to / is just plain boundedness. 

This notion can be generalized in an obvious way into relative boud- 
edness with respect to a family of operators: An operator S is said to 
be a-bounded relatively to Ti, . . . , Tk if nD{Tj) C D{S), and 

Wh e D{S), \\Sh\\ < a{\\T^h\\ + ... + \\nh\\). 

If such an a exists, then S is said to be bounded relatively to Ti, . . . , T^, 
and this will naturally be abbreviated into 

S 4 Ti, . . . , Tfc. 

1.4. Abstract Sobolev spaces. The study of partial differential 
equations often relies on Sobolev spaces, especially in a linear context. 
If one thinks of the Hilbert space ?i as a (weighted) space, there 
is a natural abstract definition of "Sobolev norm" adapted to a given 
abstract coercive symmetric operator L = A*A: define the Ti'^-Sobolev 
norm || ■ ||-^fe by 

\\hrn.:= \\hr+j2\\^%\"- 

1=1 
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Here is a generalization: When some operators Co, . . . , Cat are given 
(playing the same role as derivation operators along orthogonal direc- 
tions in M"), one can define 
(1.1) 

N k N 

Mic^ := whf + w^MW Mn^ ■■= ii^f + E E mYhf. 

Of course, there is an associated scalar product, which will be denoted 
by or 

1.5. Calculus in M". Most of the examples discussed below take 
place in M"; then I shall use standard notation from differential cal- 
culus: V stands for the gradient operator, and V- for its adjoint in 
L^(R"'), which is the divergence operator. 

Example 1. Let x = (xi, . . . , x„) and v = {vi, . . . , w„) stand for 
two variables in M". Let A = Vi,, then is the usual Hessian operator 
with respect to the v variable, which can be identified with the matrix 
of second-order differential operators [d"^ / dvjdvk) {j,k G {1, . . . , tt,}). 
Similary, if a and b are smooth scalar functions, then [aVv, bVx] is the 
matrix of differential operators [adt,^,bdxf.]- 

The scalar product of two vectors a and b in or will be 
denoted either by (a, b) or by a • b. The norm of a vector a in R" or 
C" will be denoted simply by \a\, and the Hilbert-Schmidt norm of an 
n X n matrix M (with real or complex entries) by |M|. 

The usual Brownian process in M" will be denoted by {Bt)t>o- 
The notation H'^ will stand for the usual Sobolev space in M": 
explicitly, ||M||^fe = ^j<fc II V'''"IIl2- Sometimes I shall use subscripts 
to emphasize that the gradient is taken only with respect to certain 
variables; and sometimes I shall put a reference measure if the refer- 
ence measure is not the Lebesgue measure. For instance, ||m||^i(^) = 

2. Operators L = A*A + B 

For the moment we shall be concerned with linear operators of the 
form 

(2.1) L:=A*A + B, B* = -B, 

to be thought as the negative of the generator of a certain semigroup 
iSt)t>o of interest: St = e~^^. (Of course, up to regularity issues, any 
linear operator L with nonnegative symmetric part can be written in 
the form ()2.1|) : but this will be interesting only if A and B are "simple 
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enough".) Here below I have gathered some properties of L which can 
be expressed quite simply in terms of A and B. 

2.1. Dirichlet form and kernel of L. Introduce 

/C := KerL, 11 := orthogonal projection on /C, n"*" = 1 — 11. 

Proposition 2. With the above notation, 

(i) Wh e D{A*A) n D{B), 3? {Lh, h) = \\Ahf; 

(ii) IC = KerAnKerS. 

Proof. The proof of (i) follows at once from the identities 

{A*Ah, h) = {Ah, Ah) = \\Ahf, 3? {Bh, h) = 0. 

It is clear that Kei A fl Keri? C JC. Conversely, if h belongs to JC, 
then = ^{Lh, h) = \\Ahf, so h e Ker A, and then Bh = Lh-A*Ah = 
0. This concludes the proof of (ii). □ 

2.2. Nonexpansivity of the semigroup. Now it is assumed 
that one can define a semigroup (e~*-^)j>o, i.e. a mapping (t, h) i — > 
e~^^h, continuous as a function of both t and h, satisfying the usual 
rules e^^ = Id, e~*^*+'*''^ = e~^^e~^^ for t, s > (semigroup property), 
and 

d 



\Jh e D{L), - 



e-'^h = -Lh. 

t=o+ 

As an immediate consequence, for all h G D{A*A) fl D{B) 
1 d 



2 dt 



||e'*^/if = -^{Lh,h) = -p/if < 0. 

t=o+ 

This, together with the semigroup property, the continuity of the semi- 
group and the density of the domain, implies that the semigroup is 
nonexpansive, i.e. its operator norm at any time is bounded by 1: 

Vt>0 i|e-*^||«^^, < 1. 

2.3. Derivations in L^{fi). In most examples considered later, 
the Hilbert space Ti. takes the form L^(/ioo), for some equilibrium mea- 
sure /ioo = Poo{x) dx on W^, with density poo with respect to Lebesgue 
measure; V = M*", A = {Ai, . . . , Am), and the Aj^s and B are deriva- 
tions on M", i.e. there are vector fields aj{x) and b{x) on R" such 
that 

Ajh = ttj ■ V/i, Bh = b- Wh. 

To write things symbolically, there is an m x n matrix a = cr{x) such 
that 

A = aV. 
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Below are some useful calculation rules in that context. It will be 
assumed that everything is smooth enough: For instance p^o lies in 
C^(]R") and it is positive everywhere; and a, h axe C^. The notation a* 
will denote the transpose (adjoint) of a. 

Proposition 3. With the above notation and assumptions, 

(i) B* = -B^\/- ihp^) = 0; 

(ti) A*g = -V ■ {a*g) - {V\ogp^,a*g). 

Remark 4. As a consequence of PropositionEl^ii), the linear second- 
order operator —L = XI ~ (-^ + X] ^i^j) has the Hormander form 
(a sum of squares of derivations, plus a derivation). The form A*A + B 
is however much more convenient for the purpose of the present study 
— just as in |19j . 

Proof of Proposition El By polarization, the antisymmetry of 
B is equivalent to 

V/i e n, {Bh, h) = 0. 

But 

{Bh,h)= [ {h-Vh)hp^ = \! b-V{h^)po, 

(2.2) =-l f h'V-ihp^). 

If V ■ (bpoo) = 0, then the integral in ()2.2j) vanishes. If on the other 
hand V ■ (&Poo) is not identically zero, one can find some h such that 
this integral is nonzero. This proves statement (i). 

To prove (ii), let g : W ^ R™ and /i : M" ^ R, then {A*g,h) 
coincides with 

{g, Ah)= g- (aVh) Poo = - {cr*g) ■Vhp^ = - V ■ {(T*gp^)h 
Jr" J Jr" 

= - V-{a*g)hpoo- a*g ■ {V\ogpoo)hpoo: 
7r" jr" 

where the identity Vpoo = (V log poo)Poo was used. This proves (ii). □ 

The following proposition deals with the range of applicability for 
diffusion processes. 

Proposition 5. Let a e C2(R"; M"^*") and { g C^(R";R"), and 
let {Xt)t>o be a stochastic process solving the autonomous stochastic 
differential equation 

dXt = V2a{Xt)dBt + aXt)dt, 
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where {Bt)t>o is a standard Brownian motion in M™". Then 

(i) The law {pt)t>o of Xt satisfies the diffusion equation 

(2.3) = V-(DVp-ep), D:=a*a; 

(ii) Assume that the equation ()2.3p admits an invariant measure 
fioo{dx) = poo{x) dx (with finite or infinite mass), where poo lies in 
C^(]R") and is positive everywhere. Then the new unknown h(t,x) : = 
p(t,x)/pooix) solves the diffusion equation 

(2.4) _ = V ■ (D V/i) - - 2D V log Poo) ■ V/i, 

which is of the form dth + Lh = with L = A*A + B, B* = —B, if one 
defines 

(2.5) n:=L\fi^y, A:=aV; S := - W logpoo) ■ V. 

Proof. Claim (i) is a classical consequence of Ito's formula. To 
prove claim (ii), write 

dh 1 

137 = — V • (op^Vh + D/iVpoo - ipooh) 
at \ J 

= V ■ {DVh) + 2DVh ■ _ ^ . y/i + — [v ■ {DVp^) - V ■ {pU) 

Poo Poo '- 

As Poo is a stationary solution of ()2.3j) . the last term in square brackets 
vanishes, which leads to ()2.4j) . Define A and B by ()2.5j) . Thanks to 
Proposition El (ii), it is easy to check that 

A*Ah = -V ■ (DVh) - DV log Poo ■ V/i, 

so h indeed satisfies dth + Lh = 0. It only remains to check that 
B* = —B. By Proposition El (i), it is sufficient to check that 

V ■ (poo(e - /^V log Poo)) = -V ■ (Z^Vpoo - ^Poo) 

vanishes; but this follows again from the stationarity of poo- So the 
proof of Proposition El is complete. □ 

2.4. Example: The kinetic Fokker— Planck equation. The 

following example will serve as an important application and model. 
Consider a nice (at least C^) function : M" — *■ M, converging to 
+00 fast enough at infinity (say V{x) > i^|a;|" — C for some positive 
constants K and C). For x,v e x M", set 

g-[l/(x)+H!] 

foo{x,v) := , fj.{dxdv) = foo{x,v) dx dv, 
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where Z is chosen in such a way that /i is a probabihty measure. Define 

L := -A^ + vS/^ + v-V.^- VV{x) ■ V,„. 
The associated equation is the kinetic Fokker-Planck equation with 
confinement potential V , in the form 

(2.6) dth + v- y^h - yV{x) ■ W^h = A^h - v ■ Wyh. 

Before considering convergence to equihbrium for this model, one 
should first solve analytical issues about regularity and well-posedness. 
It is shown by Helffer and Nier |30l Section 5.2] that ()2.6p generates a 
C°° regularizing contraction semigroup in L'^ifi) as soon as V itself lies 
in C°°(M"'). To study this equation for a less regular potential V, it is 
always possible to regularize V into a smooth approximation V^, then 
perform all a priori estimates on the regularized problem, and finally 
pass to the limit as e ^ 0. The following well-posedness theorem 
justifies this procedure by forcing the convergence of the approximate 
solutions to the original solution. 

Theorem 6. Let V G C\W), iniV > -oo, and let 

E{x,v) := V{x) -\ — Poo = jji{dx dv) = p^{x,v) dv dx. 

Then, for all Hq G L'^ifJ'), equation (j2.6p admits a unique distribu- 
tional solution h = h{t, x, v) G C(M+; I?'(M^ x R^)) n L^^(^+; L'^ifi)) n 
Ll^{R+] Hl{ii)), such that /i(0, ■) = Hq. 

The proof of existence is a straightforward consequence of a stan- 
dard approximation procedure, the Helffer-Nier existence results, and 
the a priori estimate 

J h'^{t, x,v) dfx{x,v) + J J h'^{s, x,v) dn(x,v) ds 

There is more to say about the uniqueness statement, of which the 
proof is deferred to Appendix IA.20I The main subtlety lies in the 
absence of any growth condition on W; this is overcome by a local- 
ization argument inspired from |30l Proposition 5.5]. Apart from that. 
Theorem ini is just an exercise in linear partial differential equations. 
Many people (including me) would rather think of ()2.6|) in the form 



(2.7) dj + v VJ - VVix) ■ VJ = A J + V, ■ (vf), 
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in which case / at time t can be interpreted (if it is nonnegative) as 
a density of particles, or (if it is a probabihty density) as the law of a 
random variable in phase space. To switch from ()2.6p to ()2.7|) it suffices 
to set / := fooh- This however does not completely solve the problem 
because the natural assumptions for ()2.7p are much more general than 
for ()2.6|) . For instance, it is natural to assume that the initial datum 
/o for ()2.7|) is with polynomial weight; or just L^, or even a finite 
measure. Theorem [7| below yields a uniqueness result in such a setting, 
however with more stringent assumptions on the initial datum. In the 
next statement, M(]R" x R") stands for the space of finite measures 
on X M", equipped with the topology of weak convergence (against 
bounded continuous functions). 

Theorem 7. Let V e C^{W), iniV > -oo, and let E{x,v) := 

I |2 

V{x) + -i^. Then, for any fo G L'^{{1 + E)dxdv), equation ()2.7|) 
admits a unique distributional solution f = f{t,x,v) G C(R_|_;C(M^ x 
K)) n L^c{^+; L\{1 + E) dxdv)) n Ll^{^+- Hl{W^ X R-)), such that 
/(0,-) = /o. 

// moreover V^V^ is uniformly hounded, then for all finite mea- 
sure fo the equation ()2.7|) admits a unique solution f = f(t,x,v) G 
C(R+;M(R^ X R;^)). 

The proof of this theorem will be deferred to Appendix IA.20I 

3. Coercivity and hypocoercivity 
3.1. Coercivity. 

Definition 8. Let L be an unbounded operator on a Hilbert space 
Ti, with kernel /C, and let 7i be another Hilbert space continuously and 
densely embedded in K.^, endowed with a scalar product (■, and a 
Hilbertian norm || ■ ||^. The operator L is said to be A-coercive on Ti. if 

V/i G /C^ nD(L), ^{Lh,h)f^> X\\hf~, 

where 3? stands for real part. The operator L is said to be coercive on 
7i if it is A-coercive on Ti. for some A > 0. 

The most standard situation is when 7i = /C"*" ~ 7i//C. Then it is 
equivalent to say that L is coercive on /C-*- (which will be abbreviated 
into just: L is coercive), or that the symmetric part of L admits a 
spectral gap. 

Coercivity properties can classically be read at the level of the semi- 
group (assuming it is well-defined), as shown by the next statement: 
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Proposition 9. With the same notation as in Definition\^ L is 
\-coercive on H if and only if ||e~*^/io||-^ < 6~'^*||/io||7:^ for all Hq G H 
and t > 0. 

Proof. Assume by density that Hq ^Hn D{L). On one hand the 
coercivity imphes 

A \\e~'^hof^ = -2^ {Le-'^ho, e^'^ho) < -2A ||e~*^/io||', 
dt 

so by GronwalFs lemma 

II -tLr ||2 ^ p-2At||T ||2 

Conversely, if exponential decay holds, then for any ho G TinD^L), 



3? {LhQ, ho) = lim 



I/, ||2 \\f,~tLj. ||2 

I'^oll^ — ||e "•oll^ 



t^o 2t 

>.in.i.f"-""'""'°II^^A||/.||^ 

- t^o 2t " " ' 

whence the coercivity. □ 

When an operator L is in the form (j2.ip . the coercivity of L follows 
from the coercivity of A*A, at least if B has a sufficiently large kernel: 

Proposition 10. With the notation of Subsection M.IV if A*A is 
X-coercive on (Ker A)-*- and KerA C Keri?, then L is X-coercive on 

Proof. We know that /C = Ker(A) nKer(_B) = Ker(y4), so for any 
h e /C^, {Lh, h) = \\Ahf > X\\hf. □ 

Example 11. Apart from trivial examples where B = 0, one can 
consider the following operator from jl] : 

L = -{A.^-x- V^.) - (A„ - w ■ V^) + ■ - X ■ V„) 

on i:2(e-(l^l'+l^l')/2 dxdv). 

The main problem in the sequel is to study cases in which A*A 
is coercive, but L is not, and yet there is exponential convergence to 
equilibrium (i.e. to an element of /C) for the semigroup (e~*^). In view 
of Proposition ITUl this can only happen if Ker L is smaller than Ker A. 
Here is the most typical example: With the choice Ti = L^(exp(— (|f p + 
\x\'^)/2) dv dx) again, consider 

L = -(A, - V ■ V,) + {v-V,-x- V,). 

Then Ker A is made of functions which depend only on x, but KerL 
only contains constants. 
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3.2. Hypocoercivity. To fix ideas, here is a (possibly mislead- 
ing, but at least precise) definition of "hypocoercivity" in a Hilbertian 
context. 

Definition 12. Let ?i be a Hilbert space, L an unbounded opera- 
tor on Ti generating a continuous semigroup (e~*^)t>o, and Ti another 
Hilbert space, continuously and densely embedded in /C-*-, endowed with 
a Hilbertian norm || ■ ||^. The operator L is said to be A-hypocoercive 

on Ti if there exists a finite constant C such that 

(3.1) \/hoen, Vt>0 ||e-*^/io||^ < Ce-^*||/io||^. 

It is said to be hypocoercive on Ti. if it is A-hypocoercive on Ti. for some 
A > 0. 

Remark 13. With respect to the definition of coercivity in terms 
of semigroups, the only difference lies in the appearance of the con- 
stant C in the right-hand side of ()3.1|) (obviously C > 1, apart from 
trivial cases; C = 1 would mean coercivity). The difference between 
Definition |H1 and Definition El seems to be all the thinner in view of 
the following fact (pointed to me by Serre): Whenever one has a norm 
satisfying inequality (|3.ip for some constant C, it is always possible to 
find an equivalent norm (in general, not Hilbertian) for which the same 
inequality holds true with C = 1. Indeed, just choose 

N{h) := supfe^* ||e"*^/i||). 
t>o ^ ^ 

In spite of these remarks, hypocoercivity is a strictly weaker concept 
that coercivity. In particular, hypocoercivity is invariant under change 
of equivalent Hilbert norm on Ti, while coercivity is not. This has an 
important practical consequence: If one finds an equivalent norm for 
which the operator L is coercive, then it follows that it is hypocoercive. 
I shall systematically use this strategy in the sequel. 

Remark 14. It often happens that a certain space Ti is convenient 
for proving hypocoercivity, but this particular space is much smaller 
than Ti. (stated otherwise, the Hilbert norm on Ti cannot be bounded 
in terms of the Hilbert norm on 7i) : typically, Ti may be a weighted 
Sobolev space, while His a weighted space. In that situation there 
is in general no density argument which would allow one to go directly 
from hypocoercivity on Ti, to hypocoercivity on Ti. However, such 
an extension is possible if L satisfies a (hypoelliptic) regularization 
estimate of the form 

(3.2) Wto > 3C(to) < +oo; Wt > to, ||e"*^||^i^^ < C{to); 
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or, more generally, if L generates a semigroup for which there is expo- 
nential decay of singularities: 

'Vt>0, e-*^ = 5t + i?i, 
(3.3) ivto>0 3C(to)<+oo; > t^, \\St\\,^^_f^ < Cito); 



[3X>0; Vt>0, \\Rt\\K^^K^<Ce 



-At 



Such assumptions are often satisfied in realistic models. For instance, 
integral operators (generators of jump processes) usually satisfy ()3.3j) 
when the kernel is integrable (finite jump measure), and ()3.2j) when the 
kernel is not integrable. Diffusion operators of heat or Fokker-Planck 
type usually satisfy (j3.2|) . 

3.3. Commutators. If the operators A*A and B commute, then 
so do their exponentials, and e~*^ = e~*^*'^e~*'^. Then, since B is 
antisymmetric, e~*^ is norm-preserving, and it is equivalent to study 
the convergence for e~*^ or for e~^^*^. On the other hand, if these 
operators do not commute, one can hope for interesting phenomena. 

Proposition 15. With the notation of Sub section ] l.li in particular 
L = A*A + B, define recursively the iterated commutators 

Co '■= A, Ck := [Ck-i,B], 

and then K' := (IkyoKeT Ck- Then /C C /C', and JC' is invariant for 

Proof. Assume that /C C KerCo fl . . . n KerCj. Then, for all 

Helens, 

Ck+ih = CkBh - BCkh = CkBh = ~CkA*Ah = 0. 

Thus K, C KerCfc+i. By induction, /C is included in the intersection /C' 
of all KeiCj. 

Next, if /i e /C' n 5, then Lh = Bh, so CkLh = CkBh = Ck+ih + 
BCkh = 0; since k is arbitrary, in fact Lh G /C', so L leaves /C' invariant, 
and therefore so does e~*^. □ 



In most cases of interest, not only does /C' coincide with /C, but in 
addition /C' can be constructed as the intersection of just finitely many 
kernels of iterated commutators. Thanks to the trivial identity 

k / k 

f|KerQ = Ker K^qC, 

j=o \i=o 
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the condition that K,' is the intersection of finitely many iterated com- 
mutators may be reformulated as 

(3.4) There exists N^e^ such that Ker C*C^ = Ker L. 

Example 16. For the kinetic Fokker-Planck operator ()2.()|1 . Nc = I 
will do. 

If the goal is to derive estimates on the rate of convergence, it is 
natural to reinforce the above condition into a more quantitative one: 

(3.5) ^^C^Cfc is coercive on /C"*". 

fc=0 

Condition ()3.5|) is more or less an analogue of Hormander's "rank r" 
bracket condition (as explained later, r = 2Nc + 1 is the natural con- 
vention) , but in the context of convergence to equilibrium and spectral 
gap, rather than regularization and elliptic estimates. There is how- 
ever an important difference: Here we are taking brackets always with 
B, while in Hormander's condition, brackets of the form, say, [Ai, Aj] 
would be allowed. This modification is intentional: in all the cases of 
interest known to me, there is no need to consider such brackets for 
hypocoercivity problems. A basic example which will be discussed in 
Appendix IA.19I is the following: The differential operator 

L ■= -{x^dy*dy + (9^*9^), 

although not elliptic, is coercive (not just hypocoercive) in L^(7)/]R, 
where 7 is the gaussian measure on M^. For this operator, brackets of 
the form [dxixdy] play a crucial role in the regularity study, but they 
are not needed to establish lower bounds on the spectral gap. 

Remark 17. It was pointed out to me by Serre that, when Ti is 
finite-dimensional, condition (j3.5|) is equivalent to the statement that 
Ker A does not contain any nontrivial suhspace invariant by B. In the 
study of convergence to equilibrium for hyperbolic systems of conser- 
vation laws, this condition is known as Kawashima's nondegener- 
acy condition j34|, I29|, I45j . It is not so surprising to note that the 
very same condition appears in Hormander's seminal 1967 paper on 
hypoellipticity p. 148] as a necessary and sufficient condition for 
a diffusion equation to be hypoelliptic, when the constant matrices A 
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and B respectively stand for the diffusion matrix and the hnear drift 
function.^ 

Taking iterated commutators may rapidly lead to cumbersome ex- 
pressions, because of "lower-order terms". In the present context, this 
might be more annoying than in a regularity context, and so it will be 
convenient to allow for perturbations in the definition of Ck, say 

[Ck, B] = Ck+i + Rk+i, 

where -R^+i is a "remainder term", chosen according to the context, 
that is controlled by Co, . . . , C^. An easy and sometimes useful gener- 
alization is to set 

[Ck, B] = Zk+iCk+i + -Rfc+i, 

where the Z^'s are auxiliary operators, typically multipliers, satisfying 
certain identities. 

Once the family (Co, . . . , CatJ is secured, one can introduce the 
corresponding abstract Sobolev IH} norm as in (jl.lj) . This norm will 
be used on 7i, or (more often) on /C-*-. On the latter space we may also 
consider "homogeneous Sobolev norms" such as 

(3.6) \\h\\%:=Y.\\CM?- 

j=o 

Note that, with the above assumptions, the orthogonal space to 
the kernel JC in Ti.^ does not depend on whether we consider the scalar 
product of H or that of H^. (See the proof of Theorem |211 below.) So 
a natural choice for Ti. will be 

Ti. = Ti.^ /}C, 

which is /C"*" equipped with the Ti.^ norm. 

4. Basic theorem 

In this section linear operators satisfying a "rank-3" condition (c = 
2 in ()3.5p ) are considered. Although this is a rather simple situation, 
it is already of interest, and its understanding will be the key to more 
complicated extensions; so I shall spend some time on this case. Here 

^At first sight, it seems that both problems arc completely difFcrent: 
Kawashima's condition is applied to systems of unknowns, while Hormander's ex- 
ample deals with scalar equations. The analogy becomes less surprising when one 
notices that for such a diffusion equation the fundamental solution, viewed as a 
function of time, takes its values in the finite-dimensional space of Gaussian distri- 
butions, so that the equation really defines a system. 
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it will be assumed for pedagogical reasons that the operators A and C 
commute; this assumption will be relaxed in the next section. 

Theorem 18. With the notation of Subsection \l.l] consider a lin- 
ear operator L = A*A + B (B antisymmetric) , and define C := [A, B]. 
Assume the existence of constants a, (3 such that 

(i) A and A* commute with C; A commutes with A (i.e. each Ai 
commutes with each Aj); 

(a) [A, A*] is a-bounded relatively to I and A; 

(Hi) [B,C] is j3-bounded relatively to A, A? , C and AC; 

Then there is a scalar product ((■, ■)) on IH} /K,, which defines a norm 
equivalent to the IH} norm, such that 

(4.1) V/i G 1-1}/1C, {{h, Lh)) > K{\\Ahf + \\Chf) 

for some constant K > 0, only depending on a and (3. 
If, in addition, 

A*A + C*C is K-coercive 

for some k > 0, then there is a constant A > 0, only depending on a,P 
and K, such that 

\/hen}/lC, {{h,Lh))> X{{h,h)). 

In particular, L is hypocoercive in Ti^/lC: 

\\e~'^\\nWK^n^/Jc<ce~^' (c < +oo), 

where both A and c can be estimated explicitly in terms of upper bounds 
on a and (3, and a lower bound on n. 

Before stating the proof of Theorem ^1 I shall provide some re- 
marks and further explanations. 

Remark 19. Up to changing a and /9, it is equivalent to impose 
(ii) and (iii) above or to impose the seemingly more general conditions: 
(ii') [A, y4*] is a-bounded relatively to /, A and A*, 
(iii') [B, C] is /^-bounded relatively to A, A^, A*A, C and AC- 

Indeed, 

{A*h, A*h) = {AA*h, h) = {A*Ah, h) + J2 ([^»' A*]^' 

i 

SO 

\\A*h\\^ < WAhW^ + \\[A, A*]h\\ \\h\\. 
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Then assumption (ii') imphes 

||A*/if < ||A/if + a(||/if + \\Ah\\ \\h\\ + \\A*h\\ \\h\\^ 

< \\Ahf + a(^\\hf + \\Ah\\ + ip*/if + ^||/if, 

and then it follows that A* is bounded relatively to / and A, so that 
(ii) holds true. This also implies that A*A is bounded relatively to 
and A, so (iii') implies (iii). 

Remark 20. Assumption (ii) in Theorem IT^ can be relaxed into 
(ii") [A, A*] A is relatively bounded with respect to A and A"^; 
where by convention 

\\[A,A*]Ahf = J2\\j2i^^^^*M^^f- 

i 3 

Remark 21. Here is a crude heuristic rule explaining a bit the 
assumptions (i) to (iii) above. As is classical in Hormander's theory, 
define the weights w{0) of the operators involved, by 

w{A) = w{A*) = l, w{B)=2, w{[0i,02]) =w{0i)+w{02). 

Then rules (i) to (iii) guarantee that certain key commutators can be 
estimated in terms of operators whose order is strictly less: for instance, 
the weight of [B, C] is 2 + 3 = 5, and assumption (iii) states that it 
should be controlled by some operators, for which the maximal weight 
is 4. (This rule does not however explain why I is allowed in the 
right-hand side of (ii), but not in (iii); so it might be better to think 
in terms of Assumption (ii") from Remark 1201 rather than in terms of 
Assumption (ii).) 

Remark 22. In particular cases of interest, it may be a good idea 
to rewrite the proof of Theorem^l taking into account specific features 
of the problem considered, so as to obtain better constants A and C. 

4.1. Heuristics and strategy. The proof of Theorem ITHl is quite 
elementary; in some sense, the most sophisticated analytical tool on 
which it rests is the Cauchy-Schwarz inequality. The argument con- 
sists in devising an appropriate Hilbertian norm on Ti^/JC, which will 
be equivalent to the usual norm, but will turn L into a coercive oper- 
ator. One can see an analogy with a classical, elementary proof of a 
standard theorem in linear algebra [3 pp. 147-148]: If the real parts of 
the eigenvalues of a matrix M are all positive, then e~*^ (expo- 
nentially fast) as t — * oo. 
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Define 

(4.2) ([h,h)) = \\hf + a\\Ahf + 2b^{Ah,Ch) +c\\Ch\\\ 

where the positive constants a, b, c will be chosen later on, in such a 
way that 1 3> a 3> 6 ^ c. (The constant c here is not the same as the 
one in the conclusion of Theorem 1181) 

By polarization, this formula defines a bilinear symmetric form on 
Ti^. By using Young's inequality, in the form 



2b{Ah,Ch) < 2b\\Ah\\ \\Ch\\ < b^^WAhf + b^^\\Ch\\ 

one sees that the scalar products ((■,■)) ^^'^ define equivalent 

norms as soon as 6 < ^/ac, and more precisely 

(4.3) 

min(l, a, c) I 1 = I < {{h, h)) < max(l, a, c) I 1 H — = ) 

\ V ac/ \ vac/ 

In particular, the scalar products ((■, ■)) and (■, define equivalent 
norms. 

In spite of their equivalence, the scalar products ((■, ■)) and (■, 
are quite different: it is possible to arrange that L is coercive with 
respect to the former, although it is not with respect to the latter. 
Heuristically, one may say that the "pure" terms H^/ip and 

||C/i|p will mainly feel the influence of the symmetric part in L, but 
that the "mixed" term {Ah, Ch) will mainly feel the influence of the 
antisymmetric part in L. The following simple calculations should help 
understanding this. Whenever Q is a linear operator commuting with 
A (be it /, y4 or C in this example), 

^ WQe-'^'^hf = -2\\QAh\\\ 

but on the other hand 
d 



<^Ae~'^h,Ce-'^h) = -{ABh,Ch) - {Ah,CBh). 

t=o 



dt 

Pretend that B and C commute, and this can be rewritten 

-{ABh,Ch) - {Ah,BCh) = -{ABh,Ch) - {B*Ah,Ch) 

= -{ABh,Ch) + {BAh,Ch) 

= -{[A,B]h,Ch) = -\\Ch\\\ 

where the antisymmetry of B has been used to go from the first to the 
second line. This will yield the dissipation in the C direction, which 
the symmetric part of A was unable to provide! 
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4.2. Proof of Theorem llHL Introduce the norm ()4.2j) . By Propo- 
sition ^[ any h E K, = KerL satisfies Ah = 0, Ch = 0, in which case 
((/i, h')) = {h, h')y^i = {h, h'). In particular, the orthogonal space /C"*" is 
the same for these three scalar products. So it makes sense to choose 

il = v}/ic. 

Let us compute 
2 at 

if we can bound below this time-derivative by a constant multiple of 
((e~*^/i, e~^^lifj , then the conclusion of Theorem^Jwill follow by Gron- 
wall's lemma. By semigroup property, it is sufficient to consider t = 0, 
so the problem is to bound below ^{{h,Lh)) by a multiple of ([h,h)). 
Obviously, 

(4.4) 3? {{h, Lh)) = 3? {h, Lh) + a (I) + 6 (II) + c (III), 
where 

(I) := 3ft {Ah, ALh), (II) := 3ft {ALh, Ch) + 3ft {Ah, CLh), 

(III) ■=^{Ch,CLh). 

ByPropositionEfi), 3ft(/i,L/i) = \\Ahf. For each of the terms (I), (II), (III), 
the contributions of A*A and B will be estimated separately, and the 
resulting expressions will be denoted (I)^, (I)^, (11)^, (H)^, etc. For 
consistency with the sequel, I shall introduce the notation 

(4.5) R2:=[C,B]. 

Moreover, to alleviate notation, I shall temporarily assume that ?i is a 
real Hilbert space; otherwise, just put real parts everywhere. 
First of all, 

(I)^ = {Ah, ABh) = {Ah, BAh) + {Ah, [A, B]h) 

= 0+{Ah,Ch) > -\\Ah\\\\Ch\\, 

where the antisymmetry of B was used. Then, 

(I)^ = {Ah,AA*Ah) = {A^h,A^h) + {Ah, [A,A*]Ah), 

to be understood as 

{AjAih, AiAjh) + {Aih, [Ai, A*j\Ajh). 
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This can be rewritten 
J2 WAAjhW^ + {[Aj, Ai]h, AiAjh) + {A,h, [Ai, A*]Ajh) 

= \\A^hf + {[A, A]h, A^h) + {Ah, [A, A*]Ah). 

In the present case it is assumed that [A, A] = 0, so the second term 
vanishes. Then from the Cauchy-Schwarz inequahty we have 

(I)^ > p^/if - \\Ah\\ \\[A,A*]Ah\\. 

Next, 

(II) ^ = {ABh,Ch) + {Ah,CBh) 

= {ABh, Ch) + {Ah, BCh) + {Ah, [C, B]h) 
= {ABh, Ch) - {BAh, Ch) + {Ah, R2h) 
= {[A,B]h,Ch) + {Ah,R2h) 
>\\Chr-\\Ah\\ WR^hW; 

(11)^ = {Ah, CA*Ah) + {AA*Ah, Ch) 

= {Ah, A*CAh) + {A*A^h, Ch) + {[A, A*]Ah, Ch) 
= {A^h, CAh) + {A^h, ACh) + {[A, Ch) 
= 2{A^h, CAh) + {Ch, [A, A*]Ah) 
> -2\\A'^h\\\\CAh\\ - \\Ch\\\\[A,A*]Ah\\. 

(Here the commutation of C with both A and A* was used.) 
Finally, 

(III) ^ = {Ch,CBh) = {Ch,BCh) + {Ch, [C,B]h) 

= + {Ch,R2h) 
> -\\Ch\\ \\R2h\\; 



(III)^ = {Ch,CA*Ah) = {Ch,A*CAh) = {ACh,CAh) = \\CAh 
(here again the commutation of C with A and A* was used). 
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On the whole, 
(4.6) 

^{{h,Lh)) > \\Ahf 

+ a(^\\A^hf - \\Ah\\ \\[A,A*]Ah\\ - \\Ah\\ \\Ch\\ 

+ b[\\Ch\\^-\\Ah\\ \\R2h\\-2\\A^h\\\\CAh\\ - \\Ch\\\\[AA*]Ah\\ 

+ c(^\\CAhf- ||C/i||||i?2/^ 
The assumptions of Theorem imply 

\\[AA*]y\\<a{\\y\\ + \\Ay\\), 

\\R2h\\ < f3{\\Ah\\ + \\A^h\\ + \\Ch\\ + \\CAh\\). 

Plugging this into (|4.6|) . follows an estimate which can be conveniently 

recast as 

^{{h,Lh))>{X, mX)K4, 
where X is a vector in and m is a 4 x 4 matrix, say upper- diagonal: 

X := (^\\Ah\\,\\A^h\\,\\Ch\\,\\CAh\\^, 

-{aa + bp) -(aa + bP) -{a + ba + b(3 + c(3) -b(3 

a -iba + cp) -2b 

b-c(3 -cP 

c 



m := 



[f^ij\l<i,j<i 



If the symmetric part of m is definite positive, this will imply in- 
equality ()4.H) . Then the rest of Theorem follows easily, since the 
K-coercivity of A*A + C*C implies 

\\Ahr+\\chr>^-{\\Ahr+\\chr)+^ 



min(l, k) ,, , 
> h 



|2 



2 

So it all boils down now to choosing the parameters a, b and c in such 
a way that the symmetric part of m is positive definite, and for this it 
is sufficient to ensure that 

Vz, rriii > 0; 

In the sequel, the statement "the symmetric part of mi is greater than 
the symmetric part of 1712" will be abbreviated into just "mi is greater 
than 7712" ■ 
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(4.7) l>a>h>2c. 
Then m can be bounded below by 
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1 - 2Ma -2Ma 
a 



If now it is further assumed that 

1 



(4.8) 



a < 



4M' 



-4Ma -Mb 

-2Mb -2Mb 

b-Mc -Mc 
c 



c < 



2M' 



then the latter matrix can in turn be bounded below by 



1/2 


-2Ma 


-AMa 


-Mb 





a 


-2Mb 


-2Mb 








b/2 


-Mc 











c 


fhij 


< \/fhi 


imjj/2 < 


{fhii + 



By imposing 
(4.9) 
it will follow 

^ mijXiXj > ^ fhiiX'^ ^''^^ rhuXf. 

ij i i 

(The 3 in 3/4 is because each diagonal term should participate in the 
control of three off-diagonal terms.) To ensure ()4.9j) . it suffices that 

2Ma < 4Ma < a/ — , Mb<J-, 2Mb < \ 

- V 8' - V 16' -1/8' - V 8 ' 



Mc < 




2Mb < 

All these conditions, including ()4.8p . are fulfilled if 

1 a2 52 I 



T' ^ - 256 M2' 



Lemma IA.16I in Appendix IA.22I shows that it is always possible to 
choose a, b, c in such a way. This concludes the proof of Theorem | 
□ 



34 



Remark 23. There are other possible ways to conduct these cal- 
culations. In an early version of this work, the last term (IH)^ was 
rewritten in three different forms to create helpful terms in ||/lC*/i|p 
and ||/l*C/i|p, at the cost of requiring additional assumptions on [C, C*]. 

5. Generalization 

Now I shall present a variant of Theorem which covers more 
general situations. 

Theorem 24. Let H he a Hilbert space, let A : H ^ H"- and 
B : Ti. ^ H. be unbounded operators, B* = —B, let L := A*A + B and 
K, := KerL. Assume the existence of &'N and (possibly unbounded) 
operators Co, Ci, . . . , Cnc+i, Ri, ■ ■ ■ , Rnc+i o.^'^d Zi, . . . , ^at^+i such that 

Co = A, [C„B] = Z^+,C,+i+R,+i (0<j<iVc), CAr,+i = 0, 

and, for all k ^ {0, . . . , N^}, 

(i) [A,Ck\ is bounded relatively to {Cj}o<j<fe and {CjA}o<j<fc-i; 

(ii) [CkiA*] is bounded relatively to I and {Cj}Q<j<k', 

(Hi) Rk is bounded relatively to {Cj}o<j<k-i and {CjA}o<j<fc_i. 
(iv) There are positive constants Xj, Aj such that Xjl < Zj < Ajl. 

Then there is a scalar product ((■, ■)) on IH} , which defines a norm 
equivalent to the H} norm, 



HI 



fe=0 

such that 

(5.1) When^/}C, ^{{h,Lh))>Kj2\\Cjhf 

j=0 

for some constant K > 0, only depending on the bounds appearing 
implicitly in assumptions (i)-(iii). 

If, in addition, 

CjCj is K-coercive 

for some k > 0, then there is a constant A > 0, only depending on K 
and K, such that 

\/h G n^/lC, 3? ((/i, Lh)) > A {{h, h)). 
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In particular, L is hypocoercive in Ti^/IC: There are constants C > 
and A > 0, explicitly computable in terms of the bounds appearing 
implicitly in assumptions (i)-(iii), and k, such that 

This result generalizes Theorem in several respects: successive 
commutators are allowed, remainders Rj+i and multiplicators ^j+i are 
allowed in the identity defining Cj+i in terms of Cj, and the operators 
Co, . . . ,Cnc are not assumed to commute. 

Remark 25. The same rule as in Remark 1211 applies to Assump- 
tions (i)-(iii). 

Remark 26. Theorem lA. 121 in Appendix IA.21I will show that the 
same structure assumptions (i)-(iv) imply an immediate regularization 
effect Ti H}. This effect extends the range of application of the 
method, allowing data which do not necessarily lie in IH} but only in 

n. 

Proof of Theorem [211 The proof is an amplification of the proof 
of Theorem ^1 Let 

(5.2) ((/i, h)) := ||/if + 5^(afc||Cfc/if + 23?6fc(Cfc/i, Ck+ih) ^ 

k=0 

where {ak}o<k<Nc+i and {bk}o<k<N^ are families of positive coefficients, 

satisfying 

(5.3) 

{do < 6, bk < S Ofc, Ofc+i < S bk, 
al<5bk-ibj, (l</e<A^e), bl < 5 aua^+i (0 < /e < iVj. 

The small number 5 > will be chosen later on, and the existence of 
the coefficients Ofc, bk is guaranteed by Lemma [A . 1 61 again. 

Since Cn^+i = 0, the last term in ()5.2|) . with coefficient ^at^, does 
not play any role. For k < Nc — 1, the inequality b^ < ^^0^0^+1 implies 

{Ckh,a^,h)\ < ^-Y\\Ckhr+^-^\n^^hr. 

Hence, for 6 small enough, 

N. 



k=0 

So the norm defined by ()5.2|) is indeed equivalent to the norm. 
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The next observation is that the space /C is the same, whether the 
orthogonahty is defined with respect to the scalar product in 7i, the 
one in Ti.^ or the one defined by ()5.2j) . To show this it is sufficient to 
prove 

(5.4) heKeiL^ VA; e {0,...,iV, + l}, Ckh = 0. 

This will be achieved by finite induction on k. Let h G Ker L; by 
Lemma El Ah = 0, Bh = 0; so ()5.4|) is true for k = 0. Assume 
now that Cjh = for j < /c; it is obvious that also CjAh = for 
j < k; then our assumption on -Rfc+i implies that Rk+ih = 0. So 
Ck+ih = CkBh — BCkh — Rk+ih = 0. This concludes the proof of ()5.4|) . 

To prove (j5.H) it is obviously sufficient to establish 

(5.5) 5J iih, Lh)) > -WAhf + ^ -jWCk^ihr). 

k=0 

As in the proof of Theorem ^1 one can compute, with obvious 
notation, 

(5.6) 5J {{h, Lh)) = WAhf + $^{a4(I)^ + (I)^] + 264(11)^ + (11)^]}. 

k=0 

To alleviate the notation, assume for a moment that we are working 
in a real Hilbert space, so there is no need to take real parts (otherwise, 
just put real parts everywhere). Explicit computations yield, for k < 

(I)^ = {Ckh, CkBh) = {Ckh, [Ck,B]h) + {Ckh, BCkh) 
= {Ckh, [Ck,B]h)+0 
= {Ckh, Zk+iCk+ih) + {Ckh, Rk+ih) + 
> — Afc+i||Cfc/i|| ||Cfc+i/i|| — ||Cfe/;,|| 



(I)^ = {Ckh,CkA*Ah) = {Ckh,AVkAh) + {Ckh, [Ck,A*]Ah) 

= {ACkh, CkAh) + {Ckh, [Ck, A*]Ah) 
= {CkAh, CkAh) + {[A, Ck]h, CkAh) + {Ckh, [Ck, A*]Ah) 
> WCkAhf - \\CkAh\\ \\[A,Ck]h\\ - \\Ckh\\ \\[Ck,A*]Ah\\ 
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and, ioT k<N^-l (when k = N^, (11)^ = 0), 



37 



(11)^ = {CkBh,Ck+,h) + {Ckh,Ck+iBh) 

= {CkBh,Ck+ih) + {Ckh,BCk+ih) + {Ckh, [Ck+i,B]h) 
= {CkBh,Ck+ih) - {Ckh,BCk+ih) + {Ckh, [Ck+i,B]h) 
= {[Ck,B]h,Ck+ih) + {Ckh, [Ck+i,B]h) 

= {CkBh,Ck+ih) — {BCkh,Ck+ih) + {Ckh, Zk+2Ck+2h) + {Ckh, Rk+2h) 

= {[Ck,B]h,Ck+ih) + {Ckh,Zk+2Ck+2h) + {Ckh,Rk+2h) 

= {Zk+iCk+ih, Ck+ih) + {Rk+ih,Ck+ih) + {Ckh, Zk+2Ck+2h) + {Ckh,Rk+2h) 

> Afc+i llCfc+l/lll^ — ||C/fc+i/;,|| ||i?/c+i/l|| — Ai;+2||C'fc^|| ||C'a;+2^|| — ||C'/fc^|| ||-Rfc+2^ 



(11)^ = {Ckh,Ck+iA*Ah) + {CkA*Ah,Ck+ih) 

= {Ckh, [Ck+i,A*]Ah) + {Ckh,A*Ck+iAh) + {A*CkAh,Ck+ih) 

+ {[Ck,A*]Ah,Ck+ih) 
= {Ckh, [Ck+i,A*]Ah) + {ACkh,Ck+iAh) + {CkAh, ACk+ih) 

+ {[Ck,A*]Ah,Ck+ih) 
= {Ckh, [Ck+i,A*]Ah) + {CkAh,Ck+iAh) + {[A,Ck\h,Ck+iAh) 

+ {CkAh, Ck+iAh) + {CkAh, [A, Ck+i]h) + {[Ck, A*]Ah, Ck+ih) 
> -\\Ckh\\ \\[Ck+i,A*]Ah\\ - \\CkAh\\ \\Ck+iAh\\ - \\Ck+iAh\\ \\[A,Ck]h\\ 

- \\CkAh\\ \\Ck+iAh\\ - \\CkAh\\ \\[A,Ck+i]h\\ - \\Ck+ih\\ \\[Ck,A*]Ah 



The next step is to use the quantities a^- HCjA/ip and 6j||Cj+i/i||^ 
to control all the remaining terms. For this I shall apply Young's 
inequality, in the form XY < eX^ + {^C{e) = 6-^4). In the 

computations below, the dependence of C on the constants Aj will not 
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be recalled. 
(5.7) 

au[{l)\ + (I)^,] + hu[{ll)\ + (11)^] > auWCkAhW" + hu\k+i\\Ck+ih\\'' 



5.8) 




Ok-l 


5.9) 


- ehu-i\\Ckh\\'' 


C \\Rk+ih\\ 

Ok-l 


5.10) 


-eakWCkAhW 


-Cak\\[A,Ck]hf 


5.11) 




-C^\\[C,,A*]Ahr 

Ok-l 


5.12) 




-Cb,\\R,+,hf 


5.1oj 


-ehk-i\\Cuhf 


Ok-l 


5.14) 


- ehk-i\\Ckh\\'' 


-C-^WRk^.hf 

Ok-l 


5.15) 


-ehu-i\\Cuhf 


-C-^\\[Ck+,,A*]Ahf 

Ok-l 


5.16) 


-eakWCkAhf 


- C^\\Ck+iAh\\' 


5.17) 


- eau+A\Cu+iAh\\'' - C-^\\[A,Cu]hf 


5.18) 


-eak\\CkAh\\^ 


- C^WCk+iAhf 
ak 


5.19) 


-eakWCkAhf 


-c^\\[A,Ck+i]hr 

ak 


5.20) 


-£6fc||C,+i/if 


-C^\\[C,,A*]Ahf, 



ak 



with the understanding that ()5.12|) to ()5.20p in the above are not 
present when k = Nc- The problem is to show that each of the terms 
appearing in hnes ()5.8|) to ()5.20p can be bounded below by 

(5.21) -s[\\Ahr + J2MC,Ahr + &,||C,+i/.f )), 

j 

as soon as 6 is small enough. This is true, by construction, of all the 
terms appearing on the left in these lines; so let us see how to control 
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all the terms on the right. In the sequel, the notation u v means 
u < rjv, where 77 > becomes arbitrarily small as 5 0. 

- For lines (|5.8|) . (|5.13p . (|5.16|) and (|5.18p . it is sufficient to impose 
^2 52 p 



The first and the third of these inequalities are true by construction; 
as for the second one, it follows from 

bl -C ttkttk+i < \/bk-ibk \/bkbk+i =^ bl < bk-ibk+i- 



- For lines dlH^ and (jHUlD, we know that \\Rk+ih\\ is con- 

trolled by a combination of . . . , ||Cfc/i||, 

||CfcA/i||; hence it is sufficient to bound the coefficients appearing in 
front of Rk+ih (resp. Rk+2h) by a small multiple of ak (resp. afc+i). So 
these terms are fine as soon as 

«fc 7 bl 

- — < afc, hk < afc, - — -C Cfc+i. 

The second of these inequalities is true by construction, while the ffist 
and third one follow from 

al < bk-ibk < bk-ittk, bl < bk-ibk+i < 6^-10^+1. 



- For lines (jOTl|l . (jOTjl and we know that || [A, Ck]h\\ is con- 

trolled by ||Ci/i||, ||Cfc/i||. By a reasoning 

similar to the one above, it is sufficient to ensure 

52 ^2 
ak < — ^ < — < 

The ffist and third of these inequalities are true by construction, while 
the second one follows from 

bl < afcOfc+i -C fefc-iOfc+i. 



- For lines (PTTTll . (UnHll and (PT^ . we know that || [C^, is 
controlled by \\y\\, \\Ay\\, \\Ciy\\, \\Cky\\, so || [Cfc, is con- 

trolled by ||Ci/l/i||, . . . , By a reasoning similar 

to the one above, it is sufficient to ensure 

bl bl 

Ofc-l Ofc-l ^k 
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The first and third of these inequalities are true by construction, while 
the second one follows from 



Putting all together, for all rj there is a 5 such that each of the 
"error terms" which appeared in the estimates above can be bounded 
below by ()5.21|) . Then, if Ne stands for the number of error terms, 

'^fc II ^ ^ 7 ii2 , '^k 



5J {{h, Lh)) > (1 - N,e) \\Ah\\' + J^(^^\\C,Ahf + -j\\C,+ih\\ 



k=0 

which implies fl5.5|) for 6 small enough. The proof of ()5.1|) is now 
complete, and the end of Theorem 1241 follows easily, as in the proof of 
Theorem UHl □ 



I shall conclude this section with a simple generalization of Theo- 



rem 



Theorem 27. With the same notation as in Theorem \24\ define 
Cj+i/2 = CjA. Then the conclusion of Theorem \24\ still holds true if 



Assumptions (i) to (Hi) are relaxed as follows: There exists a constant 
M such that 

(V) \\[A,Ck]h\\ < M^||C«/if \\Cf3hf-\ where the sum is over 
all couples of indices (a, (3) such that 6a + {1— 6)[3 < A; + 1/2, (a, /?) 7^ 
(A; + 1/2, A; + 1/2); 

(11') \\[Ck,A*]Ah\\ < M^||C„/i|n|C^/if-^, where the sum is 
over all couples of indices such that 9a + {1 — 9)(3 < k + 1, 

{a,p) ^ {k + l,k + l); 

(Hi') \\Rkh\\ < ||Cq,/i||^ IIC^/ill"'^"^, where the sum is over all 

couples of indices {a,P) such that 9a + {1 — 9)[3 < k, (a,/?) 7^ {k,k). 

In statements (i') to (iii')? may vary from one couple {a, (3) to the 
other. The conditions on admissible couples can be understood 

more easily if one remembers Remark |^ then the weight w{Ca) is 
2a + 1. If one formally attributes to [y4,Cfc], [^*,C'a;]^ and Rk the 
weights 2/c+2, 2A; + 3 and 2A; + 1, and decides that the weight of a formal 
product \\Oxh\\^ llOa/iH^"^ is 9w{0i) + (1 - 9)w{02), then Conditions 
(i') to (iii') mean that each of the operators [A, Cfc], [y4*,Cfc]y4 and Rk 
can be bounded in terms of lower weights. For instance, an estimate 
hke 

||i?5^|| < M^\\CM \\CM\ 

is admissible, since [w{C'i) + w{C^)]/2 = (7 + ll)/2 < 11 = w{R^). 
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Proof of Theorem EH The strategy is the same as in Theo- 
rem |22 but now one should use Young's inequahty in the form 

aV"^ <ea + Ch, 

and note that, if {um) is given by Lemma fA.16| then 

777 -L ri 

<^-, (m,n)^(£,£) 



m ^n- 



Then all the estimates entering the proof of Theorem 1241 can be adapted 
without difficulty. □ 

6. Hypocoercivity in entropic sense 

In this section I shall consider the problem of convergence to equi- 
librium for solutions of diffusion equations in an LlogL setting. This 
represents a significant extension of the results already discussed, be- 
cause in many cases of interest, after a finite time the solution auto- 
matically belongs to L log L{fi), where /i is the stationary solution, but 
not to L'^ifi).^ 

For that purpose, I shall use the same information-theoretical func- 
tional as in the theory of logarithmic Sobolev inequalities: First, the 
Kullback information (or Boltzmann H functional, or Shannon in- 
formation), 



H^i{v) = J hloghdfi, u = hfi; 
and secondly, the Fisher information 



2 



Recall that a probability measure fi on satisfies a logarithmic Sobolev 
inequality if there is a constant A > such that 

for all probability measures u on (with the convention that H^{v) = 
Iu{f^) = +00 if u is not absolutely continuous with respect to fi). 

The main difference with the classical theory is that I shall distort 
the Fisher information by using a suitable field of quadratic forms; 
that is, replace / \ Vh\'^/hdfi by / {S'Vh,'Vh)/hdfi, where x —>■ S{x) is 

^Of course this docs not contradict the fact that it will be locally C°°. The 
most basic illustration is the case of the linear Fokker-Planck equation dth ~ 
Ayh — V ■ Vi,h: The theory of hypercontractivity tells us that the semigroup at 
time t is regularizing from LP to only after a time \og{{q — — l))/2, which 
is finite only if p > 1 . 
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a function valued in the space of quadratic forms, in such a way that 
S{x) > kIn for some k > 0, independently of x. It turns out that the 
same algebraic tricks which worked in a Hilbertian context will also 
work here, at the price of more stringent assumptions on the vector 
fields: The proofs will be based on some slightly miraculous-looking 
computations, which may be an indication that there is more structure 
to understand. 

Here below is the main result of this section. Note carefully that 
this is not expressed in terms of linear operators in abstract Hilbert 
spaces, but in terms of derivation operators on M^. (The theorem 
might possibly be generalized by replacing by a smooth manifold.) 
So as not to be bothered with regularity issues, I shall assume here that 
the reference density is rapidly decaying and that all coefficients are C°° 
and have at most polynomial growth; but of course these assumptions 
can be relaxed. I shall also assume that the solution is smooth if the 
initial datum is smooth. 

Theorem 28. Let E e C^(]R^), such that is rapidly decreas- 
ing, and fi{dX) = e~^^^^ dX is a probability measure on M^. Let 
(y4j)i<j<m and B be first-order derivation operators with smooth co- 
efficients. Denote by A* and B* their respective adjoints in L'^{fi), 
and assume that B* = —B. Denote by A the collection (Ai, . . . , Am), 
viewed as an unbounded operators whose range is made of functions 
valued in M™. Define 

m 

L = A* A + B = Y^ A*Aj + B, 
i=i 

and assume that e~^^ defines a well-behaved semigroup on a suitable 
space of positive functions (for instance, e~^^h and log(e~*^/i) are C°° 
and all their derivatives grow at most polynomially if h is itself with 
all derivatives bounded, and h is bounded below by a positive constant). 
Next assume the existence of N^. > 1, derivation operators Cq, . . . , Ctv^+i, 
. . . , Rjy^+i; and vector-valued functions Zi, . . . , Z^^+i (all of them 
with coefficients, growing at most polynomially, as their partial 
derivatives) such that 

Co = A, [Cj,B] = Z,+,Cj+i+Rj+i (0<j<iVc), Cjv.+i=0, 
and 

(i) [A,Ck] is pointwise bounded relatively to A; 
(a) [Ck,A*] is pointwise bounded relatively to /, {Cj}o<j<ki 
(Hi) Rk is pointwise bounded with respect to {Cj}o<j<k-i; 
(iv) there are positive constants Aj,Aj such that Xj < Zj < Aj; 
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(v) [A,Ck]* is pointwise bounded relatively to I, A. 

Then there is a function x S{x), valued in the space of quadratic 
forms on M^, uniformly bounded, such that if one defines 

S{h) := Jh\oghdfi + 1 i^^X^l^ d;., 

then one has the estimate 

—S(e < -a- -da, 

dt ^ ' - h ^ 

for some positive constant a > 0, which is explicitly computable in 
terms of the bounds appearing implicitly in conditions (i)-(v). 

If furthermore 

(a) there is a positive constant A such that '^kC^Ck > XIn, point- 
wise on M^; 

(b) fi satisfies a logarithmic Sobolev inequality with constant K; 

then S{x) is uniformly positive definite, and there is a constant n > 
such that 

dt 

In particular, 

/^((e-*^/i/i) = 0{e-^'), H,{{e-'''hfi) = 0{e-^'), 

and all the constants in this estimate can be estimated explicitly in 
terms of the bounds appearing implicitly in conditions (i)-(v), and the 
constants A, K . 

Remark 29. The matrices S{x) will be constructed from the vector 
fields entering the equation, by linear combinations with constant coef- 
ficients. I expect that for more degenerate situations it will be useful 
to use varying coefficients. 

Remark 30. A major difference between the assumptions of The- 
orem 121] and the assumptions of Theorem |2H1 is that the latter impose 
pointwise bounds on M^, in the following sense. First, A is an m-tuple 
of derivation operators (A.t)i<i<m; each of which can be identified with 
a vector field cxj, in such a way that Aif = ai ■ V/; so a = (c3"j)i<j<m 
can be seen as a map valued in (m x A^) matrices. Then each commu- 
tator Ck is also an m-tuple of derivation operators (C'fc,j)i<j<m, so that 
Cfej has been obtained from the commutation of Ck-ij with B. Then 
[Ai, Ckj] is represented by a vector field ^ij,k', and Assumption (i) says 
that \^i,j,k{x)\ is bounded, for all x, by c||(t(x)||, where c is a constant. 
The other pointwise conditions are to be interpreted similarly. Let us 
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consider for instance Assumption (ii). Since Ai is a derivation, the 
adjoint of takes the form —Ai + a^, where is a function; so the 
adjoint of A is of the form {gi, ...,gm) ^ - I] ^ifl'i + Yl (^i9i- Then 
the commutator of C with A* is the same as the commutator of C 
with A, up to an array of operators that are the multiphcation by the 
functions Cjai. So the second part of Assumption (ii) really says that 
the functions Cja^ are all bounded. Finally, note that in Assumption 
(a), each Ck is an m-tuple of derivations, so it can be identified to a 
function valued in m x matrices; and ^ C^Cfc to a function valued 
in N X N matrices, which should be uniformly positive definite. 

Remark 31. Theorem lA . 1 51 below will show that the Assumptions 
of Theorem |2H1 entail an immediate "entropic" regularization effect: 
If the initial datum is only assumed to have finite entropy, then the 
functional S becomes immediately finite. This allows to extend the 
range of application of the method to data which are only assumed to 
have finite entropy. In the case of the Fokker-Planck equation I shall 
show how to relax even this assumption of finite entropy and treat 
initial data which are only assumed to have finite moments of large 
enough order. 

The key to the proof of Theorem |2H1 is the following lemma, which 
says that the computations arising in the time-differentiation of the 
functional S are quite the same as the computations arising in Theo- 
rem 1^ provided that A and Ck commute. 

In the next statement, I shall use the notation 



for the time- derivative of the functional along the semigroup gener- 
ated by the linear operator —S. More explicitly. 



Moreover, when no measure is indicated this means that the Lebesgue 
measure should be used. 



Lemma 32. Let ^{dX) = e~^WrfX, A = {Ai,...,Am), B and 
L = A* A + B be as in Theorem\M Let C = (C^, . . . , C„) and C = 
{C[, . . . , C^) be m-tuples of derivation operators on (all of them 
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with smooth coefficients whose derivatives grow at most polynomially). 
Then, with the notation f = he~^ , u = log h, one has 

(6.1) ) [ h loghdfi = 0; 



dt 



B 



(6.2) -(^^^^^lh\oghdfi = l\^dfi= I f\Au\' 

where by convention l^up = 



(6.3) 



d\ f {Ch,C'h) ^ f {Ch, [C',B]h) ^ f {[C,B]h, C'h) ^ 

dt)j -^'^=j — h — ^^+y — h — 



(6.4) = J f{Cu, [C\B]u) + J f{[C,B]u, C'u), 

where by convention {[C, B]u, C'u) = '^j{[Cj,B]u){CjU); 



+ f{[C,A*]Au, C'u) + J f{Cu, [C',A*]Au) 
+ (^J f{CAu, [AC]u) + J f{[A,C]u, C'Au) 

+ j fQA,C,C'{u), 

where by convention {Cu, [C',A*]Au) = Y.ijiCj'^)i{Cj, A*]Aiu) , etc. 
and 

(6.6) Qa,c,c'{u) := [A, CY{Au ® C'u) + [A, C']*{Au ® Cu) 

:= ^[A„C,]*(^*wC» + [A,,C'^*iA,uC,u). 

Remark 33. If = [Aj,Cj] = for all then obviously 

Qa,c,c' vanishes identically. The same conclusion holds true if A = C = 
C", even if [AjjAj] is not necessarily 0. Indeed, [Aj, Ai\*{AjU Aiu) = 
— [Ai, Aj]*{AiU Aju), so ^^j[Ai, Aj]*{AiU Aju) = by the symmetry 
i ^ j. I don't know whether there are simple general conditions for the 
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vanishing of Qa,c,c'i that would encompass both [A, C] = [A,C] = 
and A = C = C a.s particular cases. 

Remark 34. One of the conclusions of this lemma is that the 
derivatives of the quantities / h log hdfi and / {Ch, C'h)/hdfi can be 
computed just as the derivatives of the quantities / h"^ dfi and / {Ch, Ch) dfi, 
if one replaces in the final result the measure fi hj f = h ^, and in 
the integrand the function h by its logarithm, as long as the quantity 
Qa,c,c' vanishes. In the special case when C = C = A, this principle 
is well-known in the theory of logarithmic Sobolev inequalities, where 
it is stated in terms of Bakry and Emery's calculus" . As in the the- 
ory of r2 calculus, Ricci curvature should play a crucial role here, since 
it is related to the commutator [A, A*]. (In a context of Riemannian 
geometry, this is what Bochner's formula is about.) 

Proof of Lemma E21 The proofs of ()6.H) and ()6.2|) are easy and 
well-known, however I shall recall them for completeness. The proof 
of (j6.3|) will not cause any difficulty. But the proof of (j6.5|) will be 
surprisingly complicated and indirect, which might be a indication that 
a more appropriate formalism is still to be found. 

As before, I shall assume that the function h is very smooth, and 
that all the integration by parts or other manipulations needed in the 
proof are well justified. To alleviate notation, I shall abbreviate e~*^h 
into just h, with the understanding that the time dependence is im- 
plicit. Also / = he~^ and u = log h = log f + E will depend implicitly 
on the time t. Recall that the Lebesgue measure is used if no integra- 
tion measure is specified. 

First, by the chain-rule, 

— ^ J h\oghd^ = J (log h + 1) (Bh) dfi = J B{hlogh)dfi 

{B*l){h\ogh) d/i, 

and this quantity vanishes since B* = —B is a derivation. This 
proves ()6.H) . 
Next, 

jh\oghdfi = j {log h+l){A*Ah)dfi= j{A{logh+l), Ah)dfi 

— — , Ah) da. 
h ' 



This proves ()6.2I 
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To prove (jfi.Hj) . it suffices to remark that (a) the integrand can be 
written as a quadratic expression of \fh: Indeed, by chain rule, 

I (S^KCh) ^ 4 1 ^^^^ 

and that (b) the evolution equation for ^JTi along B is the same as for h: 
Indeed, dth+Bh = implies dtVh+By/h = 0. So to compute the time- 
derivative in ()6.3|1 the problem reduces to a quadratic computation: 

J {CBVh, C'Vh) dfi = j {BCVh, C'Vh) dfi+ j ([C, B]Vh, C'Vh) dfi. 

Then the first term in the right-hand side vanishes since B is antisym- 
metric. Formula (|6.3p follows upon use of the chain-rule again. 

Now it only remains to establish (j6.5|) . Before starting the compu- 
tations, let us recast the equation dth + A*Ah = in terms of / = he~^ . 
It follows by Proposition with poo = and -8 = 0, that 

dJ = V- {DiVf + / VE)) = V ■ P/Vm), 

with the diffusion matrix D = A* A, or more rigorously a*a, where cr 
is such that Ah = (t(V/i). In particular, if v and w are two smooth 
functions, then 

(6.7) {DWv, Ww) = {Av, Aw). 

Another relation will be useful later: by explicit computation, if (yf is a 
vector-valued smooth function, then 

A*g = V-{a*g)-{aVE, g); 

it follows that, for any real-valued smooth function u, 

(6.8) V ■ (DVu) - {DVE, Vn) = A*Au. 
Next, by chain-rule. 



h 

So the left-hand side of ()6.5p is equal to 
(6.9) - / V ■ (/ DVu) {Cu, C'u) - j f (^ CV-{DfVu) ^ 

C'V ■ {DfVu 



f{Cu, 



f 
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The three terms appearing in the right-hand side of ()f).9j) will be 
considered separately. First, by integration by parts and (j6.7|) . 

(6.10) - y" V ■ ifDVu) {Cu, C'ti) = j f (pVu, V{Cu, C'u) 

(6.11) = / jIau, A{Cu, C'u) 



For the second term in ()6.9j) . we use the identity 
V-iDfVu) = fViDVu)+{DVu, V/) = fViDVu)+f{DVu, Vlog/). 
So 



j f (^C{DVu, Vlog f)^N, C'u 



= -J f{cV ■ (DVu), C'u) - J f (c{DVu, Vu)^N, C 

+ j f(c{DVu, VE)^N,C'u 

By combining the first and third integrals in the expression above, then 
using ()6.8|1 and ()6.7j) again, we find that 

= - j f(^C{W-{DVu)-{DVE, Vu)^n), C'u)^^- J f (c{DVu, Vu)^n,C'u 

= J f {CA*Au, C'u) - J f {C\Au\^, C'u) 
(6.12) 

/ {[C, A*]Au, C'u) + [ f {A*CAu, C'u) -2 [ f {{CAu) ■ (Au), C'u). 



(In the last term, the dot is just here to indicate the evaluation of the 
matrix CAu on the vector Au. Also {A*CAu, C'u) should be under- 
stood as Y^iji^iCjAiU, C'ju) .) 
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Now the second integral in ()fi.l2j) needs some rewriting. By using 
the chain rule as before, and the definition of the adjoint, 

(6.13) 

/ {AVAu, C'u) = J {AVAu, hC log h) dfx 



{AVAu, C'h) 
{CAu, AC'h) dfi 

(6.14) = J {CAu,[A,C']h)dii+ J {CAu, C'Ah)dfi. 
The first term in ()6.14|1 can be rewritten as 

(6.15) [ {CAu, [A,C']ti)hdfi= [ f{CAu, [A,C']u). 



As for the second term in ()6.14|) . since C is a derivation, it can be 
recast as 



{CAu, C'ihAu)) dfi = J {CAu, hC'Au) dfi + J {CAu, (C'h) ® Au) dfi 
(6.16) = [ f {CAu, CAu) + [ f {{CAu) ■ (Au), C'u). 



Note that there is a partial simplification with the last term of ()6.12j) 
(only partial since the coefficients are not the same). 

Of course, the expressions which we obtained for the second term 
in (|6.9p also hold for the third term, up to the exchange of C and C 
After gathering all these results, we obtain 

(6.17) 

- J f{Cu, C'u) =Jf {Au, A{Cu, C'u)) 

(6.19) f {CAu, [A,C']u) + j f {[A,C]u, CAu) 

(6.20) +2 j f{CAu, CAu) 

(6.21) -iff {{CAu) ■ Au, C'u) + / / {{CAu) ■ Cu, Au) 
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The terms appearing in ()fi.l8j) , (jfi.l9j) and ()fi.2()|l coincide with some 
of the ones which appear in (j6.5|) , so it only remains to check that the 
ones in ()6.17p and ()6.2H) add up to ()6.6p . By using the identity 

A{Cu, C'u) = {ACu) ■ {C'u) + (ACu) ■ (Cu), 

we see that the sum of ()(i.l7j) and ()(i.2H) can be recast as 



/ (^Au, {{AC - CA)u, C'u)^ + j f (^Au, {{AC - C'A)u, Cu)^ 

fl^Au, {[A,C]u, C'u)) + j f(^Au, {[A,C']u, Cu)); 
or, more exphcitly: 

(6.22) J2 [ fiA,u)i[A,,C,]u)iC;u) + J2 [ / 



It remains to check that ()6.22|1 can be transformed into ()6.6|1 . Consider 
for instance the first term in ()fi.22j) . for some index Since [Ai, Cj] 

is a derivation, 



f{A,u){[A,C,]u)iC;u) = I {A,u)i[A,C,]h){C;u)d^^ 

h[Ai,Cj]*{AiuCjU) dfi 



f[A,C,r{A,uqu). 

This concludes the proof of Lemma D 

Proof of Theorem EHl Here I shall use the same conventions 
as in the proof of Lemma 1221 The functional S will be searched for in 
the form 

S{h) = j fu + Y,(^ak j f\Cku\^ + 2bk J f{CkU, Ck+iu)^ . 

In other words, the quadratic form S in the statement will be looked 
for in the form 

where Ck is identified with a function valued in m x matrices. 
If the inequalities ()5.3|) are enforced, then for 6 small enough 

{S{x)^,0m>Kj2\Ck{x)em; 

then 5* will be a nonnegative symmetric matrix. 
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Next, we consider the evolution of 8 along the semigroup. As re- 
called in Lemma IS21 

d f f\Ah\''^ 
— - / niogndjj, = + / — ; — djjL. 
at . I h 



Next, 

d f \Ckh 
"di 



j^^d^, = 2{{V)\ + {l)%)^ 



d f {Ckh, Ck+ih) k , fJJ^k 

-d/i = (II)^+ (11)^, 



. dt J h 

where the subscript A indicates the contribution of the A* A operator, 
and the subscript B indicates the contribution of the B operator. The 
goal is to show that these terms can be handled in exactly the same 
way as in Theorem |23] Everything can be controlled in terms of the 
quantities 

(6.23) [ f\C,Au\' and / f\Cku\' = [ rf/i. 



h 

(These integrals play the role that the quantities ||Cfc74/i|p and ||C/t/i|p 



were playing in the proof of Theorem | 

The terms (I)^ and (11)^ are most easily dealt with. By Lemma [ 
we just have to reproduce the result of the computations in the proof 
of Theorem |2S and divide the integrand by h. So in place of 

{Ckh, Zk+iCk+ih) dfi + j {Ckh, Rk+ih) dfi. 



we have 



/T\k f {Ckh, Zk+iCk+ih) f {Ckh, Rk+ih) 
Wb= / r dfi+ / d/i. 



h 

Then we proceed just as in the proof of Theorem |211 By Cauchy- 
Schwarz inequahty (applied here for vector- valued functions). 



— r — djj \ / d/i. 



Then \Rk+ih\ can be bounded pointwise in terms of \Ah\, . . . , \Ckh\, so 
/ \Rk+ih\'^ /hd^i can be controlled in terms of / \Cjh\^/hdfi for j < k. 
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The treatment of (11)^ is similar: 



{Zk+iCk+ih,Ck+ih) f {Rk+ih, Ck+ih) 

— h — ^^+y — h — 

J h . h 



h \ h \ h 



Then once again, one can control the functions \Rk+2h\ by \Cjh\ for 
i <k + l. 

Now consider the terms coming from the action of A* A. Let us first 
pretend that the extra terms Qa,c,c' in ()6.5|) do not exist. Then by 
Lemma |221 again, 



(I)^ = J f\CkAu\' + J f{[Ck,A*]Au,Cku) + j f{CkAu, [A,Ck]u). 
By Cauchy-Schwarz inequality (for vector- valued functions). 



(!);;> j f\CkAu\'-Jj f\[Ck,A*]Au\^J I f\Cku\'^ 



f\CkAu\\l I f\[A,Ck]u\'. 



Then Assumption (iii) implies 



(6.24) J I f\[Ck,A*]Au\' <cU I f\A^u\^ + Jj f\ACM' 



Finally, 

(11)^ = 2 j f{CkAu,Ck+iAu)+ j f{[Ck,A*]Au, Ck+iu)+ j f{CkU, [Ck+uA*]Au) 
+ [ f{CkAu,[A,Ck+i]u)+ [ f{[A,Ck]u,Ck+iAu). 
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and this can be bounded below by a negative multiple of 



flCkAuW / f\Ck+iAu 




f\CkAu\^^ J f\[A,C,+,]u\'-^ J \[A,Ck]u\'^ J 

then one can apply ()6.24|) (as it is, and also with k replaced by A; + 1) 
to control the various terms above. 

All in all, everything can be bounded in terms of the integrals ap- 
pearing in (j6.23|) . and the computations are exactly the same as in the 
proof of Theorem |211 then the same bounds as in Theorem |21] will 
work, provided that the coefficients and bk are well chosen. The 
result is 

(6.25) ^£(,,,<_A'/<£(£»ZM£)>,^M. 



dt J h{x) 

Now let us see what happens if Assumptions (a) and (b) are en- 
forced. By assumption (a), we have ^ |Cfc(x),^p > A|^p, where A > 0; 
so there exists k > such that 

\Vh\ 



£ih) > J fu + K J f\Vu\^ = J h\oghdi2 + K J 



h 



|2 



Thus S will dominate both the Kullback information i7^(/i/i), and the 
Fisher information I^{hfi). 

Then, since S is uniformly positive definite, 

/ > A / ^ * - 

As a consequence, by Assumption (b), 

J h{x) ^^'"^^) ^ ^H^ihfi) = K j fu 

for some k > 0. So the right-hand side of (j6.25p controls also j fu, 
and in fact there is a positive constant k such that 

j^S{h) < -KS{h). 

Then we can apply Gronwall's inequality to conclude the proof of The- 
orem! 



It remains to take into account the additional terms generated by 
Qa,c,c' in (ESI)- More precisely, in (I)^ we should consider / / QA,Ck,Cki^); 
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and in (11)^ we should handle / f QA,Ck,Ck+i{'^)- So the problem is to 
bound also these expressions in terms of the quantities ()f).23j) . 
We start with the additional term in (I)^, that is, 

(6.26) j fQ^^c„c,{u)= j f[A,CknAu®Cku). 

By assumption [A, Ck]* is controlled by / and A, so there is a constant 

c such that 

(6.27) 

[A,Ck]*{Au^Cku) < c (^J f \A{Au Cku)\ + j f\Au®CkU 

Next, by the rules of derivation of products, 

A{AuCku) = {A^u){Cku) + {Au){ACku) 

= {A^u){Cku) + {Au){CkAu) + {Au){\A,Ck]u). 

Here as in the sequel, I have omitted indices for simplicity; the above 
equation should be understood as Ai{AiuCkju) = {AiAiu){Ck,ju) + 
{A£u){Ck,jAiu) + {Aeu){[Ai,Ckj]u). Since by assumption [A, C^] is 
controlled by {Cj}o<j<k, there exists some constant c such that the 
following pointwise bounds holds: 



\A{Au0Cku)\ < c \A^u\ \Cku\ + \Au\ \CkAu\ + + ^ \Au\ \Cj 

V 0<j<k 

Plugging this in ()6.27|) and then in ()6.26|) . then using the Cauchy- 
Schwarz inequality, we end up with 



u\ 



fQA,Ck,Ck{u 



<c{Jf \A^u\\ / f\Cku\^ + \ f \Au\\ / f\CkAu\^ 



All these terms appear in (I)^ with a coefficient. So they can be con- 
trolled in terms of ()6.2H|1 . with the right coefficients, as in the proof of 
Theorem 1241 if 

ttk -C max{^/ao bk^i, y/a^, 1, max a/^T, \/&fe-i). 

These conditions are enforced by the construction of the coefficients 
and {hj). 
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Now we proceed similarly for the additional terms in (11)^. By 
repeating the same calculations as above, we find 



fQA,C„C,+,{u 



<c{Jf \Au\\ / / \Ck+M' + \ f l^M' 



1 1 f\Au\^J I f\c,^,Au\^+ J2 Jj /i^^Py / f\CM'+ 



j<k+i 



All these terms come with a coefficient bk, and they are properly con- 
trolled by (jFT^ if 

bk < max(v^, V o-o max ^/bj^i, y/cn). 

j<k+l 

Again, these estimates are enforced by construction. This concludes 
the proof of Theorem |2H1 D 

7. Application: the kinetic Fokker— Planck equation 

In this section I shall apply the preceding results to the kinetic 
linear Fokker-Planck equation, which motivated and inspired the proof 
of Theorem CHI as well as previous works |13|. I32[ I30j . 

The equation to be studied is (|2.6p . which I recast here: 

dh 

(7.1) + ^ ■ '^^^ ~ Vy(x) ■ Vyh = A^h - V ■ Vyh] 

and the equilibrium measure takes the form 

fi{dxdv) = 7(t;)e"^('') dvdx, -f{v) = J e"^ = 1. 

Let n := A := V^, B := v ■ V ^ - \/V{x) ■ V^,. Then 

takes the form dh/dt + Lh = 0, with L = A*A + B, B* = -B. The 
kernel /C of L is made of constant functions, and the space Ti^ = H^{n) 
is the usual L^-Sobolev space of order 1, with derivatives in both x and 
V variables, and reference weight /i: 

ll^llwi = / (\Vvh{x,v)\'^ + \Vxh{x,v)\'^j fi{dxdv). 

By direct computation, 

[A, A*] = 1, C:= [A, B] = V,, [A, C] = [A\ C] = 0, 
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[B, C] = V V(x) ■ V„. 

7.1. Convergence to equilibrium in H^. In the present case, 
assumptions (i)-(iii) of Theorem ITHl are satisfied if 

(7.2) VV is relatively bounded by {/, V^.} in ^^(e"^). 

By Lemma IA.18I in Appendix IA.221 this is true as soon as there exists 
a constant c > such that 

(7.3) |W| < c(l + |Vy|). 

The other thing that we should check is the coercivity of A*A+C*C, 
which amounts to the validity of a Poincare inequality of the form 



(7.4) 



J (iv„/ip+iVx./inf//i 



> K 



h djj 



1 2 



Since /i is the tensor product of a Gaussian distribution in (for 
which the Poincare inequality holds true with constant 1) and of the 
distribution in M^, the validity of ()7.4|) is equivalent to the validity 
of a Poincare inequality (in M") 

1 2 



(7.5) 



|V,/i(x)|2e-^(") dx>\ 



he 



This functional inequality has been studied by many many authors, 
and it is natural to take it as an assumption in itself. Roughly speak- 
ing, inequality ()7.5j) needs V to grow "at least linearly" at infinity. In 
Theorem lA.ll in Appendix lA. 191 1 recall a rather general sufficient con- 
dition for ()7.5j) to be satisfied; it holds true for instance if ()7.3|1 is true 
and |VV^(x)| — ^ oo at infinity. Then Theorem^! leads to the following 
statement: 

Theorem 35. Let V he a C"^ potential in , satisfying condi- 
tions ()7.3|1 and ()7.5|) . Then, with the above notation, there are con- 
stants C > and A > 0, explicitly computable, such that for all 
ho e H\fi), 



-tL 



hn 



hQ dji 



< Ce 



h 



0||//i(/^)- 



Remark 36. Conditions ()7.3|) and (j7.5p morally mean that the po- 
tential V should grow at least linearly, and at most exponentially fast 
at infinity. These conditions are more general than those imposed by 
Helffer and Nier |30j ^ in that no regularity at order higher than 2 is 

^This comparison should not hide the fact that the estimates by HelfFer and 
Nier were aheady remarkably general, and constituted a motivation for the genesis 
of this paper. 
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needed, and there is no restriction of polynomial growth on V. Here is 
a more precise comparison: Helffer and Nier prove exponential conver- 
gence under two sets of assumptions: on one hand, |30| Assumption 
5.6]; on the other hand, |301 Assumption 5.7] plus a spectral gap con- 
dition which is equivalent to ()7.5|) . Both these assumptions 5.6 and 5.7 
contain |30| eq.(5.17)], which is stronger than ()7.3|) . Finally, the spec- 
tral gap condition is not made explicitly in |30| Assumption 5.6], but 
is actually a consequence of that assumption, since it implies flA.19.1|) . 

Proof of Theorem We already checked all the assumptions 
of Theorem except for the existence of a convenient dense subspace 
S. If V is C°°, it is possible to choose the space of all C°° functions on 
M" X M" whose derivatives of all orders vanish at infinity faster than any 
inverse power of (1 -|- |Vy|)(l + \v\). (Note that the operators appear- 
ing in the theorem preserve this space because |V^V^| is bounded by a 
multiple of 1 -|- |Vy|.) Then there only remains the problem of approx- 
imating by a C°° potential, without damaging Condition (|7.3p . This 
can be done by a standard convolution argument: let Vs := V * ps, 
where Pe{x) = e~"'p{x/e), and p is C°°, supported in the unit ball, 
nonnegative and of unit integral. Then, for all e > 0, 

|W,(x)| < sup \V^V{y)\<C sup {l + \VV{y)\). 

\x-y\<£ \x-y\<e 

But (Q implies that log(l + IVV^P) is L-Lipschitz (L = 2C), so 

' y^-'^ i + \vv{yw-' ■ 

In particular, by (|7.3|) . |V^V(y)| can be controlled in terms of | 
for y close to x. It follows 

\x-y\<e^ \VV{x) - VV{y)\ < C(l + \VV{x)\)ee^'. 

As a consequence, 

\VV,{x) - VV{x)\ < C{1 + \VV{x)\)ee^'. 

From this it is easy to deduce that 1 -|- [^^(x)! > (1 — C'e){l + 
\'W{x)\), for some explicit constant C. Then, satisfies the same 
condition ()7.3|) as V, up to replacing the constant C by some constant 
C{e) which converges to C as 5 ^ 0. 
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All in all, after replacing by V^, we can apply the first part of 
Theorem UHl to get 

(7.6) 

< {{ho, ho)), 

where he{t) = e~^^''{hQ — J Hq), and i^' is a constant independent of e. 

By the uniqueness theorem of Appendix IA.201 converges to 

h{t) = e~^^h, in distributional sense as e — 0. Also, J e~^=7 — > 
J ho e~^7. Since the left-hand side is a convex functional of h and 
converges locally uniformly to V , inequality (|7.6|) passes to the limit 
as e — i> 0. The Poincare inequality for and the definition of the 
auxiliary scalar product guarantee the existence oi K' > Q such that 

{{h{t), h{t))) + K' f {{h{s), h{s))) ds < {{ho, h)). 
Jo 

The exponential convergence of (ih{t), h{t))) to follows, and the the- 
orem is proved. □ 

7.2. Explicit estimates. As a crude test of the effectiveness of 
the method, one can repeat the proof of Theorem^] on the particular 
example of the Fokker-Planck equation, taking advantage of the extra 
structure to get more precise results. Using [A, A*] = I, one obtains 

(7.7) {{h, Lh)) > \\Ahf + a\\A^hf + b\\Chf + cWCAhf - {E), 

{E) := a(p/if \\Ch\\)+b{\\Ah\\ \\R2h\\+2\\A'^h\\ \\CAh\\ + \\Ah\\ \\Ch\\) 

+ c\\Ch\\ \\R2h\\. 

Moreover, i?2^ = ~[B,C] = V'^V ■ A; to simplify computations even 
more, assume that |VV| < M (in Hilbert-Schmidt norm, pointwise on 
M"). Then 

{E) < a{\\Ah\\^ + \\Ah\\ \\Ch\\) + b(^M\\Ahf + 2\\A^h\\ \\CAh\\ + \\Ah\\ \\Ch\\^ 

+ cM\\Ah\\ \\Ch\\ 
= {a + bM)\\Ahf + {a + b + cM)\\Ah\\ \\Ch\\ + 2b\\A^h\\ \\CAh\\ 

<{a + bM + l/^)\\Ahf + {a + b + cM)^\\Chf + a\\A^hf + -\\CAhf. 
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Since < c, the last two terms above are bounded by the terms in 
WA'^hf and \\CAhf in ((TTD; so 

{{h, Lh)) > [1 - (a + hM + 1/4)] \\Ahf +[h-{a + h + cMf] ||C/if 

On the other hand, taking into account the spectral gap assumption 
on A*A + C*C, 

ih,h)) < {2a + K~^)\\Ah\\'^ + {2c + K-^)\\Ch\\'^. 

So the proof yields a convergence to equilibrium in like 0(e~'^*), 
where 

l-{a + hM + \) h-{a + h + cMf \ 
2a + ' 2c + y ' 

and the supremum is taken over all triples (a, 6, c) with < ac. 

In the particular (quadratic) case where V'^V is the identity, one 
has M = 1, K = 1; then the choice a = b = c = 0.05 yields A = 0.025, 
which is off the true (computable) rate of convergence to equilibrium 
A = 1/2 (see 01 p. 238-239]) by a factor 20. Thus, even if the method 
is not extremely sharp, it does yield quite decent estimates."^ Note that 
the coefficients a, b, c chosen in the end do not satisfy c <^ 6 ^ a! 

7.3. Convergence in L^. Theorem ESI is stated for initial 
data. However, it can be combined with an independent regularity 
study: Under condition (|7.3|) . one can show that solutions of (|7. 1|) 
satisfy the estimate 

(7.8) 0<t< 1^ \\fit,-)\\HH,) < ^11/(0, OIIl^M- 

A proof is provided in Appendix IA.21I Combined with Theorem ESI 
this estimate trivially leads to the following statement: 

Theorem 37. Let V be a C"^ potential in M", satisfying condi- 
tions ()7.3|) and (j7.5|) . Then, with the above notation, there are con- 
stants C > and A > 0, explicitly computable, such that for all 

ho e L^ifi), 

< Ce"^*||/io||L2(^). 



A := sup min 

(a,b,c) 



t > 1 



-t.L 



hn 



ho djjL 



a comparison, the bounds by Herau and Nier |32l formula (4)] yield a lower 
bound on A which is around 10"''. 
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7.4. Convergence for probability densities. Write e~^(x)7(f) = 
^-E{x,v)^ and set / = e~^h, then the Fokker-Planck equation ()7.1|) be- 
comes the kinetic equation for the density of particles: 



The previous results show that there is exponential convergence to 
equilibrium as soon as 



As an integrability estimate, this assumption is not very natural for 
a probability density; as a decay estimate at infinity, it is extremely 
strong. The goal now is to establish convergence to equilibrium under 
much less stringent assumptions on the initial data, maybe at the price 
of stronger assumptions on the potential V. 

An obvious approach to this problem consists in using stronger 
hypoelliptic regularization theorems. For instance, it was shown by 
Herau and Nier |32j that if the initial datum in ()7.1|) takes the form is 
only assumed to be a tempered distribution, then the solution at later 
times lies in and in fact takes the form V e~^g, where g is C°° 

with rapid decay. Similar results can also be shown by variants of the 
method exposed in Appendix IA.21I for instance one may show that 
if the initial datum belongs to a negative L^-Sobolev space of order k 
then for positive times the solution belongs to a positive L^-Sobolev 
space of order k', whatever k and k'. In particular, this approach works 
fine if the initial datum for ()7.9|1 is a probability measure /o satisfying 



However this still does not tell anything if we assume only polynomial 
moment bounds on /q. 

In the next result (apparently the first of its kind), this problem will 
be solved with the help of Theorem |2H1 that is, by using an entropy 
approach. 

Theorem 38. Assume that V is C~ with | V^V(x)| < Cj for allj > 
2; and that the reference measure satisfies a logarithmic Sobolev 
inequality. Let /o be a probability measure with polynomial moments of 
all orders: 



(7.9) 



dt 



+ V ■ VJ - WV{x) ■ = V, ■ (V,,/ + fv). 
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Then the solution to ()7.9j) is C°° in x and v for all positive times, and 
converges to exponentially fast as t ^ oo, in the sense that 

J fit,x,v)\og dxdv = 0(e-"*) (t > 1), 

with explicit estimates. 

Remark 39. A well-known sufficient condition for to satisfy 
a logarithmic Sobolev inequality is V = W + w, where V^iy > 
K > 0, and w is bounded (this is the so-called "uniformly convex + 
bounded" setting). 

Proof of Theorem EHl Since by assumption the function V\^(x) 
is Lipschitz by assumption, it can be shown by standard techniques that 
the Fokker-Planck equation admits a unique measure-valued solution. 
So it is sufficient to establish the convergence for very smooth initial 
data, with rates that do not depend on the smoothness of the initial 
datum, and then use a density argument. 

Since V'^V is bounded, the transport coefficients appearing in (|7.9|) 
are Lipschitz (uniformly for (x, v) G M" x M"), and it is easy to show by 
classical estimates that all moments increase at most linearly in time: 

It is shown in App endix lA.21l that /(t, ■) also belongs to all Sobolev 
spaces (in x and v) for t > 0; in fact, estimates of the form 

ll/(t,-)llH£HKMsxMs) = 0(r^('='^)) o<t<i 

will be established there. Then by elementary interpolation, /(t, ■) lies 
in all weighted Sobolev spaces for all t G (0, 1): That is, 

||/(t,-)|k. := ||/(t,x,t;)(l + \v\' + \x\Y'\\h^ < +oo. 

It is shown in |49| Lemma 1] that /(/) < C||/||j|/fc for k and s large 
enough (depending on n), where /(/) stands for the Fisher information, 
J /|V(log/)p//. So f(t, ■) has a finite Fisher information (in both x 
and v variables) for all t > 0. Since also /(t, ■) has all its moments 
bounded and IV-E] = 0(1 + |x| + \v\), we have in fact 

J /|v(iog/ + E)|' = o(r^) o<t< 1 

for some 7 > 0, where the time variable is omitted in the left-hand 
side. So from time t = to > on, the solution / has a finite relative 
Fisher information with reference measure n{dx dv) = e~'^*^^'^^ dxdv. 
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Then we can apply Theorem 1281 with A = Vu, B = v - Vx — W{x) ■ 
V,, Ci = [A, B] = V,, R, =0,C2 = 0, i?2 = V^V{x) ■ V,, = I. 
Assumptions (i), (ii), (iii) and (v) in Theorem EH are automatically 
satisfied, and Assumption (iv) is also satisfied since V'^V is bounded. 
(This is the place where the boundedness of the Hessian of V is crucially 
used.) Since the reference measure fi is the product of e~^^^^ dx (which 
satisfies a logarithmic Sobolev inequality by assumption) with 7(f) dv 
(which also satisfies a logarithmic Sobolev inequality), fi itself satisfies 
a logarithmic Sobolev inequality. So Theorem |2H1 yields the estimate 

I /(log/ + E) = 0(e-^(*-*°)) 

for t > to- In words: The relative entropy of the solution with respect 
to the equilibrium measure converges to exponentially fast as t — cxd. 
This concludes the proof of Theorem EHl D 

8. The method of multipliers 

A crucial ingredient in the treatment of the Fokker-Planck equa- 
tion was the use of the mixed second derivative CAh = V^^xh to 
control the error term [B,C]h = (V'^V) ■ V^h. There is an alterna- 
tive strategy, which does not need to use CA: It consists in modify- 
ing the quadratic form (j4.2j) thanks to well-chosen auxiliary operators, 
typically multipliers. In the case of the Fokker-Planck equation, this 
method leads to less general results; it is however of independent inter- 
est, and can certainly be applied to many equations. In this section I 
shall present a variant of Theorem ITHl allowing for multipliers, and test 
its applicability to the Fokker-Planck equation. Some extensions are 
feasible, but I shall not consider them. 

Let again A and B be as in Subsection 11.11 and C = [A, B], R2 = 
[C,B]; assume that [A,C] = for simplicity. Let M, TV be two self- 
adjoint, invertible nonnegative operators such that 

[B, M] = 0, [B, N] = 0, [M, N] = 

(these conditions can be somewhat relaxed by imposing only an ade- 
quate control on the commutators, but this leads to cumbersome cal- 
culations). Instead of (j4.2p . consider the quadratic form 

(8.1) ([h, h)) = ||/if + a\\MAh\\^ + 2b{MAh, NCh) + c||iVC/i|p. 

By straightforward variants of the calculations performed in the proof 
of Theorem ^[ one obtains 

iKLh)) = WAhf + aWMA^hf + h\\^/MNChf + cWNCAhW - {E), 
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where 

-(E) := a(^{MAh, MCh)+{MAh, [M, A*]A'^h)+{[A, M]Ah, MA^h) 

+ {MAh,M[A,A*]Ah)^ 

+b(^{MAh, N[B, C]h)+2{MA^h, NCAh) + {[A, M]Ah, NCAh) + {MAh, [N, A*]CAh) 
+ {[A, N]Ch, MA^h) + {NCh, [M, A*]A^h) + {NCh, M[A, A*]Ah)^ 

+c(^{NCh, NR^h) + {[A, N]Ch, NCAh) + {NCh, [N, A*]CAh)'^ . 

It would be a mistake to use Cauchy-Schwarz inequality right now. 
Instead, one should first "re-distribute" the multipliers M and on the 
two factors in the scalar products above. For instance, {MAh, MCh) is 
first rewritten {Ah, M^Ch) since it should be controlled by (inter alia) 
not \\MAhf. To obtain the correct weights, one is sometimes 
led to introduce the inverses M-^ and A^~^ In the end, 

(E) < a(^\\Ah\\ \\M^Ch\\ + \\Ah\\ \\M[M, A*]A^h\\ + \\[A, M]Ah\\ \\MA^h\\ 

+ \\Ah\\ \\M'[A,A*]Ah\\^ 

+b(^\\Ah\\ \\MN[B,C]h\\+2\\MA^h\\ \\NCAh\\ + \\[A, M]Ah\\ \\NCAh\\ 
+ \\^h\\ \\M[N, A*]CAh\\ + \\(VMN)Ch\\ \\(^N/M)[M, A*]A'^h\\ + \\ [A, N]Ch\\ \\MA'^h\\ 

+ \\y/MNCh\\ \\y/MN[A,A*]Ah\\^ 

+ c(^+\\(VMN)Ch\\ \\(N'^/yM^/^)R2h\\ 

+ \\[A,N]Ch\\ \\NCAh\\ + \\(^^MN)Ch\\ \\(^/N/M)[N, A*]CAh\\y 

Of course, in the above y^N/M stands for N-^^^M'-^^'^ , etc. 

Repeating the scheme of the proof of Theorem it is easy to 
see that (E) can be controlled in a satisfactory way as soon as, say 
(conditions are listed in order of appearance and the notation of Sub- 
section [Ol is used), 

^/MN, [M, A*] =<; /, [A, M] ^ /, M^[A, A*] =5; /, 

MN[B, C] ^ A, [A, M] ^ /, M[N, A*] 4 N, [A, N] ^ VMN, 
(^N/M)[M, A*] ^ M, y/MN[A, A*] 4 I, (N^'^ / M'^'^)[B , C] ^ A, 
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[A, N] ^ VMN, ( ^/N/M) [iV, A*] ^ N. 

For homogeneity reasons it is natural to assume = M^. Then 
the above conditions are satisfied if 

(8.2) M^[A, A*] 4 /, M^[B, C] ^ A, 

(8.3) 

[M, A] ^ /, [M, A*] ^ /, [M^, A] ^ M^, [M^, A*] =^ M^. 

If these conditions are satisfied, then one can repeat the scheme 
of the proof of Theorem with an important difference: instead of 
+ , it is only \\Ahf + \\{^/MN)Chf = WAhW^ + \\M'^Chf 

which is controlled in the end. This leads to the following theorem. 

Theorem 40. With the notation of Subsection \l.A assume that 

[C,A] = 0, [C,A*]=0, 

and that there exists an invertible nonnegative self-adjoint bounded op- 
erator M onTi, commuting withB, such that conditions ()8.2|) and ()8.3p 
are fulfilled. Define 

(8.4) {{h, h)) = \\hf + a\\MAhf - 2b{M'^Ah, Ch) + cUM^C/if . 

Then, there exists K > only depending on the bounds appearing 
implicitly in ()8.2|) and ()8.3|) . such that 

^ {{h, Lh)) > K{\\Ahf + \\M''Ch\\^). 

If in addition 

(8.5) A*A + C*M^C admits a spectral gap k > 0, 

then L is hypocoercive on H^/JC: there exists constants C > and 
A > 0, explicitly computable, such that 

||e~*'^||-Hi//c^7^i/'c — Ce~^\ 

As usual, it might be better in practice to guess the right multipliers 
and re-do the proof, than to apply Theorem 1^ directly. It is also clear 
that many generalizations can be obtained by combining the method of 
multipliers with the methods used in Theorem |211 Rather than going 
into such developments, I shall just show how to apply Theorem EUl on 
the Fokker-Planck equation with a potential V G C^(M"'). In that case, 
[A*, A] = I and [B, C] = (VV)A. When VV is bounded and A*A + 
C*C is coercive, there is no need to introduce an auxiliary operator M: 
the choice M = J is sufficient to provide exponential convergence to 
equilibrium. But a multiplier might be useful when is unbounded. 
Assume, to fix ideas, that V behaves at infinity like 0(|xp"''") for some 
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a > 0, and |V^l^| like 0(|x|"); then it is natural to use an operator 
M which behaves polynomially, in such a way as to compensate the 
divergence of V . In the rest of the section, I shall use this strategy 
to recover the exponential convergence for the kinetic Fokker-Planck 
equation under assumptions ()8.6|) and ()8.7p below. 

Let M be the operator of multiplication by m(x, f), where 

m(x, v) :- 



Vo + V{x) + ^ j 



and Vo is a constant, large enough that Vq + V is bounded below by 1. 
Since Bm = and B is a. derivation, it is true that B commutes with 
M. Assume that Vq + is bounded below by a multiple of 1 + Ixl""*"^; 
then 

4 1 K 
m < 7i — < 



for some constant > 0, and then m^(V^y) is bounded, so that 
M^[B,C] is relatively bounded by A. Finally, condition ()8.3|) reduces 
to 

|Vt,m| < Km, 

which is easy to check. To summarize, conditions (|8.2j) and (|8.3p are 
fulfilled as soon as there exist constants C > and K > Q such that 
(8.6) 

VxgM", ^(x) > i^lxp+^-C, |VV(x)| < (a> 

To recover exponential convergence under these assumptions, it re- 
mains to check the spectral gap assumption ()8.5p ! This will be achieved 
under the following assumption: there exists a potential W , and con- 
stants C > 0, -fT > such that 

(8.7)VxgM", \V{x)-W{x)\<C, V^W{x)> K{l + \x\)-''. 
From (|8.6|) there exists K > {] such that 

K \ 1 

=: ^{y)'^{x). 



v>^v. > $(v:^v..). 

Now I claim that V*\E'Vx (where is a shorthand for the multiplication 
by ^) is coercive in L'^{e~^ dx), or in other words that there exists 
K > such that 
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or equivalently, that there is a constant C such that for all f & L'^{e 

J[f{x) - f{y)]^e-^'-^^e-^'-y^ dxdy<C J ^(a;)|V./(x)| V^(") dx. 
Indeed, with C standing for various positive constants, one can write 

[fix) - fivY^^e-^^'^^e-^^yUxdy <C J [f{x) - /(y)]^ g-^^'^e-^^^) dxdy 

<C j {(y^W{x))-'VJ{x),VJ{x))e-'^^^Ux 



< 



C y"(l + |a;|)"|V./(x)pe-^(^) dx, 



where the passage from the first to the second line is justified by the 
Brascamp-Lieb inequality j7l Theorem 4.1]. 

Now it is possible to conclude: the operator A*A = V* V^, is coercive 
on L'^i'y), 7 standing for the Gaussian distribution in the v variable, 
and the operator V*\E'Vx is coercive on L^(e~^). Theorem IA.2I in 
Appendix lA. 191 shows that y4*y4 + $V*\E'Vx is coercive on where 
fi is the equilibrium distribution for the Fokker-Planck equation. By 
monotonicity, A*A + C*M^C also admits a spectral gap; this was the 
last ingredient needed for Theorem EU] to apply. 

9. Further applications and open problems 

A very nice application of Theoreml24lwas recently done by Capella, 
Loeschcke and Wachsmuth on the so-called Landau-Lifschitz-Gilbert- 
Maxwell model arising in micromagnetism. Under certain simplifying 
assumptions, the linearized version of this model can be written 

' dtm = J{h — m) 
(9.1) <dth = -V AV Ah- J{h-m) 



y.h = -V-m, 

where m : — *• stands for the (perturbation of the) magnetization, 
and /i : — >• for the (perturbation of the) magnetic field; moreover, 
J is the usual symplectic operator J[xi,X2,X3] = [—X2,Xi]. Obviously, 
the system ()9.H) is dissipative but strongly degenerate, since the dissi- 
pation term — VAVAh only acts on h, and not even on all components 
of h. This case turns out to be particularly degenerate since one needs 
three commutators to apply Theorem |^ For further details I refer the 
reader to the preprint by Capella, Loeschcke and Wachsmuth jU]. 
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Still, many issues remain open in relation to the hypocoercivity of 
operators of the form A* A + B. I shall describe four of these open 
problems below. 

9.1. Convergence in entropy sense for rapidly increasing 
potentials. In the present paper I have derived some first results of 
exponential convergence to equilibrium for the kinetic Fokker-Planck 
equation based on an entropy method (Theorem I38|). While these re- 
sults seem to be the first of their kind, they suffer from the restriction 
of boundedness imposed on the Hessian of the potential. It is not clear 
how to relax this assumption in order to treat, say, potentials that 
behave at infinity like a power of that is higher than 2. A first pos- 
sibility would be to try to adapt the method of multipliers, but then 
we run into two difficulties: (a) Entropic variants of the Brascamp-Lieb 
inequality do not seem to be true in general, and are known only under 
certain particular restrictions on the reference measure (see the discus- 
sion by Bobkov and Ledoux jHl Proposition 3.4]); (b) It is not clear 
that there is an entropic analogue of Theorem IA.3I Both problems (a) 
and (b) have their own interest. 

Another option would be to try to relax the local conditions (i)- 
(iv) into global (integrated) boundedness conditions, so as to have an 
analogue of Lemma IA.18I where the reference measure would be the 
solution / of the Fokker-Planck equation. This is conceivable only if 
/ satisfies some good a priori estimates for positive times. 



9.2. Application to oscillator chains. One of the motivations 
for the present study was the hope to revisit the works by Eckmann, 
Hairer, Rey-Bellet and others on hypoelliptic equations for oscillator 
chains, modelling heat diffusion [201 HHl HB UHl HH- So far I have 
obtained only very partial success in that direction. If we try to apply 
Theorem to the model, as it is described e.g. in the last section 
of |19j . we find that the assumptions of Theorem |21] apply as soon as 

(a) the "pinning potential" Vi and the "interaction potential" V2 
have bounded Hessians; 

(b) the Hessian of the interaction potential is bounded below by a 
positive constant; 

(c) the second derivatives of the logarithm of the stationary density 
are bounded; 

(d) the stationary measure satisfies a Poincare inequahty. 
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Let us discuss these assumptions. Assumption (a) is a bit too re- 
strictive, since it excludes for instance the quartic double-well poten- 
tials which are classically used in that field; but it would still be admis- 
sible for a start; and hopefully this restriction can be relaxed later by a 
clever use of the method of multipliers. By the way, it is interesting to 
note that such assumptions are not covered by the results in |19j which 
need a superquadratic growth at infinity. Next, Assumption (b) is not 
so surprising since (as far as I know) it has been imposed by all authors 
who worked previously on the subject.^ But it is a completely open 
problem to derive sufficient conditions for Assumptions (c) and (d), 
except in the simple case where the two temperatures of the model are 
equal. This example illustrates an important remark: The range of ap- 
plication of Theorem 121] (and other theorems of the same kind) will be 
considerably augmented when one has qualitative theorems about the 
stationary measure for nonsymmetric diffusion processes. For instance, 

- When does the stationary measure satisfy a Poincare inequality? 

- Can one derive bounds about the Hessian of the logarithm of its 
density? 

The first question was addressed recently in papers by Rockner 
and Wang (see for instance |43j ) in the context of elliptic equations, 
and it looks like a challenging open problem to extend their results to 
hypoelliptic equations. The second question seems to be completely 
open; of course it has its intrinsic interest, since very little has been 
known so far about the stationary measures constructed e.g. in |20j . 



9.3. The linearized compressible Navier— Stokes system. An 

extremely interesting instance of hypocoercive linear system is pro- 
vided by the linearized compressible Navier-Stokes equations for per- 
fect gases. In this example, the noncommutativity does not arise be- 
cause of derivation along noncommuting vector fields, but because of 
the noncommutativity of the space where the unknown takes its values. 

Obtained by linearizing the nonlinear system of Section around 
the equilibrium state (1, 0, 1), the linearized compressible Navier-Stokes 



^More generally, as pointed out to me by Hairer, all existing results seem to 
require that the interaction potential does dominate the pinning potential. 
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system reads as follows: 

(dtp + V ■u = 0; 



(9.2) 



< 9tM + V(p + ^) =/iAM + /i(l-|) V(V 



u 



N ^ 



U 



Here is the dimension, (p, u, 9) are fluctuations of the density, ve- 
locity and temperature respectively, p > is the viscosity of the 
fluid and k > the heat conductivity. So it is natural to define 
7^ = L2(fi;M X X 



where Vt C M is the position domain, and 



the target space 



X 



(p, u,d) 



is equipped with the Euclidean norm 

2 o , ,o 



P 



\u\ 



9' 



which is (up to a factor —1/2) the quadratic approximation of the usual 
entropy of compressible fluids. 

Let h = (p, u, 6); it turns out that (j9.2p can be written in the form 
dth + Lh = 0, where L = A* A + B, B* = -B, and A, B are quite 
simple: 



(9.3) 



Ah = (o, ^/2jI{Vu}, Ve 
Bh= (v -u, Vp + Ve, jfV -u 



Here I have used the notation 



{Vn}, 



1 /dui 



2\dxj 



+ 



9uj 
dxi 



V ■ u 

N 



for the traceless symmetrized (infinitesimal) strain tensor of the fluid. 

The system fl9.2|) is degenerate in two ways. First, the diffusion on 
the velocity variable u does not control all directions: In general it is 
false that / |{Vm}P controls the whole of / |VmP (see the discussion 
in |15j for instance: one needs at least an additional control on the 
divergence). Secondly, there is no diffusion on the density variable p. 
This suggests to consider commutators between A : h ^ {0,0, V6') and 
B. After some computations one gets (in slightly sketchy notation) 



IA,B]^C\ + Ri; 



0,0, VV-M 



i?i = -(0,V'^,0); 



[C, B] = C2+R2; C2h = -(0, 0, VAp); R^h = -(0, -V'V-n, VA^). 



A^ 
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So the commutator Ci controls the variations of the divergence of u, 
while the iterated commutator C2 controls the variations of the den- 
sity p. However, if we try to apply Theorem |21] in this situation, we 
immediately run into problems to control the remainder R2, and need 
to modify the strategy. This problem is tricky enough to deserve a 
separate treatment, so I shall not consider it in this memoir. 

9.4. A model problem arising in the study of Oseen vor- 
tices. All the material in this subsection was taught to me by Gallay. 
Oseen vortices are certain self-similar solutions to the two-dimensional 
incompressible Navier-Stokes equation, in vorticity formulation \24\ 
12 5j . The linear stability analysis of these vortices reduces to the spec- 
tral analysis of the operator S + aB in L^(M^), where 



(9.4) 



buj = —Auj H uj . 

16 2' 



Buj = BS[G] ■ Vuj + 2BS[G^/^lo] ■ VG^/^; 
here BSfcu] is the velocity field reconstructed from the vorticity u: 

271 J^2 \x-y\ 



BSM(x) = - / \^ '> uj{y)dy, 



and is obtained from v by rotation of angle 7r/2; moreover G is a 
Gaussian distribution: G{x) = e~'^' ''^/(47r); and a is a real parameter. 

The spectral study of S* + aB turns out to be quite tricky. In the 
hope of getting a better understanding, one can decompose uj in Fourier 
series: uj = J^nez^ni^)^''"'^ ^ where (r, 9) are standard polar coordinates 
in M^. For each n, the operators S and B can be restricted to the vector 
space generated by e*""^, and can be seen as just operators on a function 
uj{r): 

{S^uj){r) = -dluj - + ^) - (1 - 7^ 
{BnUj) {r) =in {ipuj - gVtn) ; 

here g{r) = e~''^''^/47r, ip{r) = (1 — e~'"^/^)/27rr^, and fln{i^) solves the 
differential equation 

2 

Ti T 

-irn')' + —n = -uj. 

r 2 

The regime |a| — 00 is of physical interest and has already been the 
object of numerical investigations by physicists. There are two families 
of eigenvalues which are imposed by symmetry reasons; but apart from 
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that, it seems that all eigenvalues converge to infinity as |a;| ^ oo, and 
for some of them the precise asymptotic rate of divergence 0(|a|^/^) has 
been established by numerical evidence. If that is correct, this means 
that the "perturbation" of the symmetric part S by the antisymmetric, 
lower-order operator aB is strong enough to send most eigenvalues to 
infinity as |a| — > oo. Obviously, this is again a manifestation of a 
hypocoercive phenomenon. 

To better understand this stability issue, Gallagher and Gallay sug- 
gested the following 

Model Problem 41. Identify sufficient conditions on / : M ^ M, 
so that the real parts of the eigenvalues of 

La : UJ I — > {—dluj + x^uj — u;) + iafu 

in L^(R) go to infinity as \a\ oo, and estimate this rate. 

Here is how Gallagher and Gallay partially solved this problem. Set 
H = L^(M;C), A = d^uJ + XLJ, Boj = {iaf)u. Then Coj = iafu, so the 
operator A* A + C*C is of Schrodinger type: 

{A* A + C*C) u = {-dluj + x^uo -uo) + a^f^uj, 

and the spectrum of A* A + C*C can be studied via standard semi- 
classical techniques. For instance, if f'{xY = +x'^)^, /c G N, 
then the real part of the spectrum of A* A + C*C is bounded below 
like Odap''), with u = min(l, 2/k). Then a careful examination of the 
proof of Theorem ^1 yields a lower bound like Odal*^) on the real part 
of the spectrum of A* A + B. 

This estimate is optimal for A* A + C*C, but it is in general not 
so for A* A + B. For instance, if f{x) = 1/(1 + x^), then u = 1/4, 
but numerical simulations suggest that the growth is like This 
might indicate a fundamental limitation of the techniques developed in 
this part, and motivate the development of a refined analysis. 



Part II 

The auxiliary operator method 



In this part I shall present an abstract hypocoercivity theorem ap- 
plying to a linear operator L whose symmetric part is nonnegative, but 
which does not necessarily take the form A* A + B. Still it will be useful 
to decompose L into its symmetric part S and its antisymmetric part 
B. Of course, we could always define A to be the square root of S, but 
this might be an extremely complicated operator, and the assumptions 
of the A* A + B Theorems might in practice be impossible to check. 
Important applications arise when the operator S contains an integral 
part, as in the linearized Boltzmann equation. 

A classical general trick in spectral theory, when one studies the 
properties of a given linear operator L, consists in introducing an aux- 
iliary operator which has good commutation properties with L. Here 
the idea will be similar, with just an important twist: We shall look 
for an auxiliary operator A which "almost commutes" with 5* and 
"does not at all" commute with S, in the sense that the effect of the 
commutator [A, B] will be strong enough to enforce the coercivity of 
S+[A,B]*[A,B]. 

With this idea in mind, I had been looking for a hypocoercivity 
theorem generalizing, say. Theorem ^[ but stumbled on the problem 
of practical verification of my assumptions. In the meantime, Clement 
Mouhot and Lukas Neumann found a theorem which, while in the 
same spirit of Theorem has some important structural differences. 
The Mouhot-Neumann theorem is quite simple and turns out to be 
applicable to many important investigated in |38j : so in the 

sequel I shall only present their approach, with just slight variations and 
a more abstract treatment. Then I shall discuss the weak points of this 
method, and explain why another theory still needs to be developed, 
probably with slightly more sophisticated tools. At the time of writing, 
Frederic Herau has made partial progress in this direction, but still did 
not manage to get things to work properly. 
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10. Assumptions 



In the sequel, 7i is a separable Hilbert space on M or C, 5* is a 
nonnegative symmetric, possibly unbounded operator Ti ^ Ti and B 
is an antisymmetric, possibly unbounded operator Ti ^ Ti. Then 
A = {Ai, . . . , Afn) is an array of unbounded operators Ti ^ Ti. All of 
these operators are defined on a common dense domain. I shall actually 
ignore all regularity issues and be content with formal calculations, to 
be considered as a priori estimates. 

The same conventions as in Section Q will apply. Some of the as- 
sumptions below will involve \fS U for various operators U ; of course, 
this is not rigorous since is in general not invertible. To make sense 
of these assumptions, one can either consider them as a priori estimates 
for a regularized problem in which S is replaced by an invertible ap- 
proximation (something like S + e/, and one tries to get estimates 
which are independent of e); or supply them with the condition that 
\/S is invertible on the range of U (a trivial case of application is when 



The object of interest is the semigroup generated by the operator 



The next hypocoercivity theorem for L will make crucial use of the 
commutator of A and B. I shall write 



where Z is bounded from above and below, and R is some "remainder" . 

Now come a bunch of commutator conditions which will be used in 
Section ^2 Later in Section I shall make some simplifying assump- 
tions which will drastically reduce the number of these conditions; but 
for the moment I shall keep the discussion at a general level. 



U = 0). 



L = S + B. 



[A,B] = ZC + R, 



(Al) 



either [C, S] ^ ^/S 

or Vs'\c,S]4y/S 




either {A ^ VSA, C, y/SC, v^) and ^ v^, ^C) 

or v^"^[C,L] =^ 75, ^^SC 



(A3) ^[A\C] ^ v^, v^C, C, ^/SA 



(A4) (v^A* ^ v^A, C, v^) and (v^C* ^ v^C, v^) 



r either {A* 4 VSA, C, ^/SC, v^) and {[C*,S] 4 v^, v^C) 
(A5) I or (v^A* v^A, C, ^/SC, S) and {^^'\c\ S] 4 v^, v^C) 
[or [C*,S] = 

(A6) R^y/S, VSC. 

(A7) There exist constants k,c > sucli tliat for all h ^Ti, 
{Ah,ASh) > K{SAh,Ah) -c({Sh,h) + {SCh,Ch) + ||C/if). 



Here is a simple, but sometimes too restrictive, sufficient condition 
for (A7) to hold (the proof is left to the reader): 

J either {A 4 VSA, C, VSC, VS) and {[S, A] 4 C,VSC,VS) 
[or ^^~\S,A] 4 v^, v^C, C 

Remark 42. Some of the assumptions (A1)-(A7) can be replaced 
by other assumptions involving the commutator [A, S]. I did not men- 
tion these alternative assumptions since they are in general more tricky 
to check that the ones which I chose. In case of need, the reader can 
easily find them by adapting the proof of the main theorem below. 



11. Main theorem 

Theorem 43 (hypocoercivity for L = S* + 5). With the same no- 
tation as in Section{T^ assume that (A1)-(A7) are satisfied. Further 
assume that 
(11.1) 

\i) 3/€,c>0; V/iG(KerL)^, 

{ASh,Ah) > K\\Ahf - c{{Sh,h) + {SCh,Ch) + \\Chf); 

^ (ii) S + A*SA + C*SC + A* A + C*C is coercive on (Ker L)^. 

Then there are constants c, A > 0, depending only on the constants 
appearing implicitly in (A1)-(A7) and (jll.ip . such that 

where Ti C (KerL)-*- is defined by the Hilbert norm 
WhWl = \\hf + \\Ahf + \\Chf. 
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Remark 44. Although Condition ()ll.l|l (i) formally resembles As- 
sumption (A7), I have preferred to state it together with (jll.l|) fii) 
because its practical verification often depends on a control of by 
{Sh, h) + \\Chf. 

Proof of Theorem HSl The proof is quite similar in spirit to 
the proof of Theorem so I shall be sketchy and only point out the 
main steps in the calculations. 

First note that (KerL)-*- is stable under the evolution by e^*'^. In- 
deed, ifke KeiL, then {d/dt){e-'^h, k) = {Le~'^h,k) = {e-'^h,L*k), 
so it is sufficient to show that L*k = 0. But Lk = implies {Sk, k) = 
{Lk, k) = 0, so k E Ker S* (here the nonnegativity of S is essential), so 
keKeiB also, and L*k = {S - B)k = 0. 

Next let 

J^[h) = \\hf + a\\Chf + 2b^{Ch,Ah) + c\\Ah\\^, 

where K stands for real part, and a, 6, c will be chosen later in such 
a way that l^a^6^c>0, Vb, b ^ ^/ac. In particular, 
jF(/i) will be bounded from above and below by constant multiples of 

~; so to prove the theorem it is sufficient to establish the estimate 
{—d/dt)T{e~^^h) > const. J^(e~*^/z). Without loss of generality, we can 
do it for t = only. In the sequel, I shall also pretend that 7i is a real 
Hilbert space, so I shall not write real parts. 

By direct computation, 

(11.2) -A ^(e-*^/,) ={Sh,h) 

+ a{CSh,Ch) + a{CBh,Ch) 
+ b{CLh, Ah) + b{Ch, ALh) 
+ c{ASh, Ah) + c{ABh, Ah). 

Now we shall estimate ()11.2p line after line. 

(1) The first line of pi.2|) is kept unchanged. 

(2) The second line of ()11.2|) is rewritten as follows: 

(11.3) a{CSh, Ch) = a{SCh, Ch) + a([C, S]h, Ch). 

Then the second term in the right-hand side of (jll.3p is estimated from 

below, either by -a\\[C,S]h\\ \\Ch\\, or by -a\\^fS~^[C,S]h\\ \\\fSCh\\- 
By Assumption (Al), these expressions can in turn be estimated from 
below by a constant multiple of 



'.(\\VSh\\ \\Ch\\ + \\Vsh\\ \\Vsch\\^. 
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(Here I used the identity {Su,u) = W^/SuW^.) 

(3) The treatment of the third hne of pi.2|) is crucial; this is where 
the added coercivity from the commutator [A, B] will show up. To 
handle the first term in this line, we write 

{CLh, Ah) = {LCh, Ah) + {[C, L]h, Ah) 

= {SCh, Ah) + {BCh, Ah) + {[C, L]h, Ah) 

= {SCh, Ah) - {Ch, BAh) + ([C, L\h, Ah). 

When we add this to the second term of the third line, {Ch,ALh) = 
{Ch, ABh) + {Ch, ASh), we obtain 

{Ch, {AB - BA)h) + {SCh, Ah) + {[C, L]h, Ah) + {Ch, ASh) 

= {Ch, {ZC + R)h) + {SCh, Ah) + {[C, L]h, Ah) + {Ch, ASh) 

> K\\Chf + {Ch, Rh) + {SCh, Ah) + {[C, L]h, Ah) + {Ch, ASh). 

So there are four "error" terms to estimate from below: 

(11.4) {Ch,Rh), {SCh, Ah), {[C,L]h,Ah), {Ch,ASh). 

- To estimate the first term in ()11.4j) . just write 

{Ch,Rh) > -\\Ch\\ \\Rh\\ 

and apply Assumption ( A6) ; it follows that there is a lower bound by 
a constant multiple of 

-b\\Ch\\ {\\\^h\\ + \\VsCh\\). 

- To estimate the second term in ()11.4|) . use the Cauchy-Schwarz 
inequality: 

{SCh, Ah) > - \\^/SCh\\ \\^/SAh\\. 

- To estimate the third term in (jll.4p . write either 

{[C,L]h,Ah)>-\\[C,L]h\\ \\Ah\\ 

or 

{[C,L]h,Ah) > - \\^/s'\c,L]h\\ \\y/SAh\\ 
and apply Assumption (A2). It results a lower bound by a constant 
multiple of 

- b{\\VSAh\\ + \\Ch\\ + \\^/SCh\\ + v^/i||) (||v^/i|| + \\VsCh\\) 

-b {\\^/Sh\\ + \\^/SCh\\)\\^/SAh\\). 
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- The fourth term in ()11.4j) is a bit more tricky: 

{Ch.ASh) = {A*Ch,Sh) 

= {[A*,C]h,Sh) + {CA*h,Sh) 
= {[A*,C]h,Sh) + {A*h,C*Sh) 
= {[A*, C]h, Sh) + {A*h, SC*h) + {A*h, [C\ S]h). 
This gives rise to three more terms to estimate: 

(11.5) {[A*,C]h,Sh), {A*h,SC*h), {A*h,[C*,S]h). 

- To handle the first term in ()11.5|1 . write 

{[A*,C]h,Sh) = {^/S[A*,C]h,^/Sh) > -\\y/S[A*,C]h\\ \\sqrtSh\\; 

then apply Assumption (A3) to bound ||v^[^*, C'J/iH. The result is a 
lower bound by a constant multiple of 

-6||v^/i||(||v^/i|| + \\^/SCh\\ + \\Ch\\ + \\v^Ah\\). 

- To bound the second term in ()11.5p . write 

{A*h,SC*h) = {y/SA*h,\^C*h) > -\\\^A*h\\ \\y/SC*h\\; 

then apply Assumption (A4) to bound these two norms separately. 
The result is a lower bound by a constant multiple of 

-b{\\^/SAh\\ + \\Ch\\ + \\\^Ch\\ + \\y/Sh\\) {\\\^Ch\\ + \\\^h\\). 

- To bound the last term in pi.5|l . one possibility is to write 

{A*h, [C*,S]h) > -\\A*h\\ \\[C* , S]h\\; 
another possibility is 

{A*h, [C*, S]h) = {^A*h, ^~\c\ S]h) > -||v^.4*;t|| \\^/S~\c*, S]h\\. 

Then one can apply Assumption (A5) to control these terms. In the 
end, this gives a lower bound by a constant multiple of 

-b{\\^^Ah\\ + \\Ch\\ + ||v^C/i|| + ||v^/i||)(||v^/i|| + \\^Ch\\). 

(4) Finally, the fourth line of pi.2p is handled as follows: 

(11.6) {Ah, ASh) + {Ah, ABh) = a{Ah, ASh) + p{Ah, ASh) 

+ {Ah, BAh) + {Ah, [A, B]h), 

where > and a + (3 = 1. The first term a{Ah, ASh) is esti- 
mated by means of Assumption (A7); the second term [3{Ah, ASh) by 
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means of Assumption (jll.ip fi): altogether, these first two terms can be 
bounded below by a constant multiple of 

c{\\y/SAh\\ + p/if ) - c{\\\^h\\ + \\y/SCh\\ + \\Ch\\^). 

Then the third term {Ah, BAh) in ()11.6|) vanishes; and the last term 
{Ah, [A, B]h) is bounded below by -\\Ah\\ \\ZCh\\ - \\Ah\\ \\Rh\\, which 
in view of Assumption (A6) can be bounded below by a constant 
multiple of 

-c\\Ah\\ \\Ch\\ - c\\Ah\\ {\\VSh\\ + \\^Ch\\). 



Gathering up all these lower bounds, we see that 

d 



2dt 



T > const. {X, mX), 



where 



X= (\\VSh\\, \\VSCh\\, \\Ch\\, \\VSAh\\, \\Ah\\^, 

m is the 5 X 5 matrix 

l-Mh- Mc -Ma - Mb -Ma - Mb -Mb -Mc 

a-Mb-Mc -Mb -Mb -Mc 

m= b-Mc -Mc 

c 

c 



and M is a large number depending on the bounds appearing in the 
assumptions of the theorem. 

Then by reasoning as in Section 0] and using Lemma (A. 161 we can 
find coefficients a,b,c > and a constant n > such that 
(11.7) 

--^ J^{e-'^h) > K({Sh,h) + {SCh,Ch) + \\Chf+{SAh,Ah)+\\Ahf 
2 at V 

By Assumption pi.l|) fii). this implies the existence of k,',k," > 
such that 



d 



2dt 



t=o 



J^{e-'^h) > K'(||/if + {Sh,h) + {SCh,Ch) + \\Chf + {SAh,Ah) + \\Ahf^ 



> K"j^{h). 

This concludes the proof. 



□ 
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12. Simplified theorem and applications 

In this section I shall consider a simplified version of Theorem |331 

Corollary 45. Let A = [Ai, . . . , Am), B, S be linear operators on 
a Hubert space Ti, and let C = [A,B]. Assume that 

A* = -A, B* = -B, C* = -C, S* = S> 0; 

[C,A]=0, [C,B] = 0, [C,S] = 0. 

Further assume that there exists k, c > such that for all h G (Ker A fl 
KeiB)^, 

(12.1) {Ah, ASh) > K{{SAh, Ah) + \\Ahf) 

- c{{Sh,h) + {SCh,Ch) + \\Chf); 

and that 

(12.2) S + C*C is coercive on (Ker An KerB)^. 
Then there exists A > such that 

where Ti C (KerL)-*- is defined by the Hilbert norm 
\\hf~ = \\hf + \\Ahf + \\Chf. 

Proof of Corollary The assumptions of the theorem triv- 
ially imply assumptions (A1)-(A6) from Sectional Assumption 112. ll 
is equivalent to the conjunction of (A7) and (jll.ip (i). Finally, (jl2.2j) 
is obviously stronger than (jll.ljl fii). □ 

Now let us make the link with the Mouhot-Neumann hypocoer- 
civity theorem |38l Theorem 1.1]. Although the set of assumptions in 
that reference is not exactly the same as in the current section, we shall 
see that under a small additional hypothesis, the assumptions in |38j 
imply the present ones. 

In HHI, the Hilbert space H is ^^(T^xM;;^), and A = V„, B = vV^, 
C = Vx', and the operator S only acts on the velocity variable v, so we 
have indeed A* = —A, B* = —B, C* = —C, and C commutes with A, 
B and S. The kernel of L is similar to the kernel of S (up to identifying 
V h{v) with (x, u) — > h{v)), and contains constant functions. Since 
C*C = —Ax has a spectral gap. Condition (jl2.2p is equivalent to the 
fact that S has a spectral gap in L^(R"), which is Assumption H.3 
in j38j . So it only remains to check (|12.1|) . which will be true as soon 
as 

(12.3) (VAV,5/i) >/s:((^VAV,/i) + ||V,/if) - c\\hf. 
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It is assumed in |38j that S, viewed as an operator on L^(M"), 
can be decomposed into the difference of two self-adjoint operators: 
S = A — K, where A is positive definite and 

(12.4) {Vyh, V^A/i) > K{Wyh, AWyh) - c||/if ; 

(12.5) V(5>0, 3c((5)>0; (V„/i, V„K/i) < V^,/if + c(5)||/if . 

Let us further assume that K is compact relatively to A, in the sense 
that 

Ve > 0, 3c{e) > 0; K <eA + c{e)I, 
or equivalently (since A = S + K) 

(12.6) V£>0, 3c(e)>0; K<eS + c{e)I. 

By using p2.4p . (jl2.5p and ()12.6|) . and denoting by c and n various 
positive constants, one easily obtains 

{Vvh,VvSh) = {\/vh,Vr,Ah) - {Vyh,V.uKh) 

> K{Vvh,AVvh) - c||/i|p - {V^h,V^Kh) 

> K{{V^h,AV^h) + ||V„/if ) - c||/if - {K/2)\\Vyhf - c\\hf 

> K{{\/vh,A\/-ah) + ||V^/if ) - c\\hf 

> K{{V.,h, SV„h) + II V^/if ) - c(||/if + (V^/i, KV^h)) 

> K{{V^h, SV^h) + II V^/if ) - c\\hf. 

This estabhshes fjl2.3p . 

Assumption ()12.6p is not made in |38j . but it is satisfied in all the 
examples discussed therein: linear relaxation, semi-classical relaxation, 
linear Fokker-Planck equation, Boltzmann and Landau equations for 
hard potentials. So all these examples can be treated by means of 
Theorem US I refer to [38] for more explanations and results about 
all these models. Mouhot and Neumann also use these hypocoercivity 
results to construct smooth solutions for the corresponding nonlinear 
models close to equilibrium, thereby simplifying parts of the theory 
developed by Guo, see e.g. |26j . 

13. Discussion and open problems 

Although it already applies to a number of interesting models. 
Theorem HHl suffers from several shortcomings. Consider for instance 
the case when S* is a bounded operator (as in, say, the linearized 
Boltzmann equation for Maxwellian cross-section), and there is a force 
term —W{x) ■ V^, in the left-hand side of the equation. Then the 
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higher derivative term in [C, L] is — V^y(x) ■ V„, which certainly can- 
not be bounded in terms of S and C; so Assumption (A2) does not 
hold. It is likely that Theorem 1321 rarely applies in practice when 
[C,B] = [[A,BIB]^Q. 

Other problems are due to Assumption (A4). This assumption will 
not hold for, say, C = Vx in a bounded domain fl C M"; indeed, in 
a slightly informal writing, C* = —C + a ■ dS, where a is the outer 
unit normal vector on dfl and dS is the surface measure on dfl. So the 
computation used in the proof of Theorem 1321 does not seem to give 
any result in such a situation.^ 

A last indication that Theorem 1431 is not fully satisfactory is that it 
does not seem to contain Theorem II 81 as a particular case, although we 
would like to have a unified treatment of the general case L = S + B 
and the particular case L = A* A + B. In fact, as the reader may 
have noticed, the choices of coefficients in the auxiliary functionals 
appearing respectively in the proof of Theorem ^1 and in the proof of 
Theorem 1321 go in the opposite way!! Indeed, in the first case it was 
\hf + a\Ahf + 1h{Ah, Ch) + c\\Chf with a > > c, while in the 
second case it was \\hf + a\\Chf + 2h{Ch, Ah) + c\\Ah\\^. 

Some playing around with the functionals suggests that these prob- 
lems can be solved only if the auxiliary operator A is comparable" to 
say in terms of order of differential operators. So if S is bounded, 
then also A should be bounded. This suggests to modify the Mouhot- 
Neumann strategy in the case when S is bounded, by choosing, instead 
of A = V^, something like A = {I - + v ■ V^y^^^V^. (I wrote 

— V ■ Vy rather than A„, because in many cases known to me, the 
natural reference measure is the Gaussian measure in M".) Then com- 
putations involve nonlocal operators and become more intricate. I shall 
leave the problem open for future research. 



^By the way, at present there seems to be no really satisfactory treatment of 
bounded domains in linearized kinetic theory, apart of course from the case of a 
periodic box. 
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Fully nonlinear equations 



In this part I shall consider possibly nonlinear equations, and I shall 
not depend on "exact" commutator identities. To get significant results 
under such weak structure assumptions, I shall assume that I deal with 
solutions that are very smooth, uniformly in time. Moreover, I shall 
only prove results of convergence like 0{t~°°), that is, faster than any 
inverse power of t. 

As in Remark the assumption of uniform smoothness can be 
relaxed as long as one has good estimates of exponential decay of sin- 
gularities, together with a stability result (solutions depart from each 
other no faster than exponentially fast). However, I shall not address 
this issue here. 

At the level of generality considered here, the rate 0{t~°°) can- 
not be so much improved, since some cases are included for which 
exponential convergence simply does not hold, even for the linearized 
equation. In many situations one can still hope for rates of convergence 
like 0(e~'^*^), as in the close-to-equilibrium theory of the Boltzmann 
equation with soft potentials |27j . If a linearized study suggests con- 
vergence like 0{e~^^) or 0(e~''**^) for a particular nonlinear model, then 
one can try to obtain this rate of convergence by putting together the 
present nonlinear analysis (which applies far from equilibrium) with a 
linearization procedure (close to equilibrium) and a subsequent linear 
study. 

This part is strongly infiuenced by my collaborations with Laurent 
Desvillettes on the convergence to equilibrium for the linear Fokker- 
Planck equation |13j and the nonlinear Boltzmann equation |15j . The 
method introduced in these papers was based on the study of second- 
order time differentiation of certain functionals; since then it has been 
successfully applied to other models jHl Our scheme of proof had 
several advantages: It was very general, physically meaningful, and 
gave us the intuition for the strong time-oscillations between hydrody- 
namic and homogeneous behavior, that were later observed numerically 
with a high accuracy |22j . On the other hand, our method had two 
major drawbacks: First, the heavy amount of calculations entailed by 
the second-order differentiations (especially in the presence of several 
conservation laws); and secondly, the particularly tricky analysis of the 
resulting coupled systems of second-order differential inequalities. 

The approach will I shall adopt in the sequel remedies these draw- 
backs: First, it only uses first-order differentiation; secondly, it confines 
many heavy computations into a black box that can be used blindly. 
The price to pay will be the loss of intuition in the proof. 
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The main result is a rather abstract theorem stated in Section 
and proven in Sectional Then I shall show how to use this abstract re- 
sult on various examples: the compressible Navier-Stokes system (Sec- 
tion the Vlasov-Fokker-Planck equation with smooth and small 
coupling (Section [T7|): and the Boltzmann equation fSectionlTHj). 

In the case of the Vlasov-Fokker-Planck equation to be considered, 
the coupling is simple enough that all the smoothness bounds appearing 
in the assumptions of the main theorem can be proven in terms of just 
assumptions on the initial data. In the other cases, the results will be 
conditional (depend on the validity of uniform regularity estimates). 

The hard core of the proof of the main result was conceived dur- 
ing the conference "Advances in Mathematical Physics" in the honor 
of Carlo Cercignani (Montecatini, September 2004). It is a pleasure 
to thank the organizers of that meeting (Luigi Galgani, Maria Lampis, 
Rossanna Marra, Giuseppe Toscani) for helping to create a fruitful and 
pleasant atmosphere of work. The main results were first announced 
two weeks later, in an incomplete and preliminary form, at the Confer- 
ence "Mathematical Aspects of Fluid and Plasma Dynamics" (Kyoto, 
September 2004), beautifully organized by Kazuo Aoki. During the 
Summer of 2006, for the purpose of various lectures in Porto Ercole, 
Trieste and Xining, I rewrote and generalized the main theorem, and 
added new applications. Additional thanks are due to Kazuo for an 
important remark about the treatment of the Boltzmann equation with 
Maxwellian diffusive boundary condition. 



88 



14. Main abstract theorem 

The assumptions in this section are expressed in a rather abstract 
formahsm. "Concrete" examples wiU be provided later in Sections ITBl 
tod 

14.1. Assumptions and main result. The theorem below in- 
volves five kinds of objects: 

- a family of normed spaces (X^, || ■ ||s)s>o; the index s can be thought 
of as a way to quantify the regularity (smoothness, decay, etc.); 

- two "differential" operators B and C, such that B is "conservative" 
and C is "dissipative" ; 

- a "very smooth" solution t ^ f{t) of the equation 

dtf + Bf = Cf, 

with values in a subset X of the intersection of all the spaces X^; 

- a Lyapunov functional S, which is dissipated by the equation 
above, and admits a unique absolute minimizer /oo; 

- a finite sequence of "nested nonlinear projections" (Jij)i<j<j] one 
can think that Uj is the projection onto the space of minimizers of S 
under J — j constraints, and in particular Uj is the map which takes 
everybody to /oo- 

The goal is to prove the convergence of f(t) to the stationary state 
/oo, and to get estimates on the rate of convergence. 

I shall make several assumptions about these various objects. Even 
though these assumptions may look a bit lengthy and complicated, 
I tend to believe that they are satisfied in many natural cases. The 
following notation will be used: 

- If A is an operator, then the image of a function f hj A will be 
denoted either by A{f) or simply by Af. 

- The Frechet derivative of A, evaluated at a function /, will be 
denoted by A'{f) or A'j-; so A'j ■ g stands for the Frechet derivative of 
A evaluated at / and applied to the "tangent vector" g. 

- The notation ||A'(/)||x^y stands for the norm of the linear oper- 
ator A'{f) : X ^ Y, i.e. the smallest constant C such that ■ 
glW < C\\g\\x for all geX. 

- Similarly, the second (functional) derivative of A, evaluated at a 
function /, will be denoted by A"{f) or A'^; so A'jr ■ {g, h) stands for the 
Hessian of A evaluated at / and applied to the two "tangent vectors" 
g and h. The notation ||A"(/)||x^y stands for the smallest constant 
C such that p'(/) ■ {g, /i)||y < C||^||x||/i||x for all g,heX. 
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Assumption 1 (scale of functional spaces). (X*, || ■ ||s)s>o is a non- 
increasing family of Banach spaces such that 

(i) is Hilbert; its norm \\ ■ \\q will be denoted by just \\ ■ \\; 
(a) The injection X"^ C X"^ is continuous for s' > s; that is, there 
exists C = C{s,s') such that 

(14.1) ll/lls<C||/||.,. 



(Hi) The family {X^)s>o is an interpolation family: For any Sq, Si > 
and 9 e [0, 1] there is a constant C = C{sq, si, 9) such that 
(14.2) 

s={i-9)s, + 9s,^ v/GX^onx^S 11/11. <c 11/11^;^ ll/llf,; 

One may think of s as an index quantifying the regularity of /, say 
the number of derivatives which are bounded in a certain norm. In the 
sequel, I shall sometimes refer informally to s as an index standing for 
a number of derivatives, even if it is not necessarily so in general. 

Assumption 2 (workspaces). X and Y are two sets such that X c 
Y C r\s>oX^ ; moreover, Y is convex and bounded in all spaces X^ . 

Assumption 3 (solution). / e C(M+;X") n Ci((0, +oo); X^) for 
all s; moreover f{t) G X for all t. (In particular f is bounded in all 
spaces X'^ .) 

In the sequel, the notation /o will be a shorthand for /(O). 
Assumption 4 (equation). / solves the equation 

(14.3) % + ^^ = 

where 

(i) B, C are well-defined on Y and valued in a bounded subset of 
X^ for all s; 

(a) For any s there is s' large enough such that B' is bounded X'^ 
X", uniformly on Y ; 

(Hi) C is Lipschitz X^ — > X", uniformly on Y , for s large enough. 

In short, B and C satisfy a "Lipschitz condition with possible loss of 
derivatives" . If X'^ is a Sobolev space of order s on a bounded domain, 
then any reasonable differential operator of finite order, with smooth 
coefficients, will satisfy these assumptions. 

Assumption 5 (stationary state), foo is an element of X, satisfy- 
ing Bfoo = C/oo = 0. 
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Assumption 6 (projections). (n.j)i<j<j are nonlinear operators 
defined on Y , with U j{Y) = {/oo}- fHj sends everybody to the station- 
ary state.) Moreover, for all j G {1, . . . , J}, 

(i) Iij{X) C CoHj = 0; 

(ii) Ilj/oo /oo) 

(Hi) For any s there is s' large enough such that (Jlj)' and ijlj)" 
are hounded X'^' X'^ , uniformly on Y . 

The last of these assumptions morally says that Uj is with pos- 
sible loss of derivatives. 

Assumption 7 (Lyapunov functional). S : Y ^ R is on Y 
viewed as a subset of X^ for s large enough. For all f one has £{f ) > 
Eiliif) > £{foo), and more precisely 

(i) For any e G (0, 1) there is > such that for all f eY , 

(14.4) £U) - £{Jlif) > Ke 11/ - njf 

(ii) For any e G (0, 1) there are K^, Cs > such that for all f eY , 



(14.5) \\U^f - /oof +^ < £{UJ) - £{f^) < \\Yi,f - /o^f-. 



Note that Hi/ and / are bounded uniformly, so these bounds be- 
come more and more stringent when e decreases. 

Assumption 8 (Key hypocoercivity assumptions). 

(i) C alone is dissipative, strictly out of the range of Hi: For any 
e > there is a constant K,. > such that for all f E X, 

(14.6) -£\f) ■ (Cf) > [£{f) - £{UJ)] 

(ii) C — B is dissipative just as well: For any s > there is > 
such that for all / G X, 



(14.7) V{f) := -£'{f) ■ {Cf - Bf) > [£{f) - £{UJ)] 



l+e 



(Hi) For any k < J and for any e > there is a constant > 
such that for all f E X , 



(14.8) V{f) + (Id - n,);,^.^ ■ {BU.f) > WiUk- Uk+i)f 



1 2+e 



Remark 46 (Simplified assumptions). In many cases of applica- 
tion, B is conservative, in the very weak sense that £'{f) ■ (Bf) = 0; 
then Assumption IHJii) trivially follows from Assumption EJi). Also 
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most of the time, Assumption IHl^iii) will be replaced by the stronger 
property 

(14.9) (Id - n,)^^^ ■ {BUjf) ' > - n,+o/ 



I 2+e 



In the sequel, I shall however discuss an important case where none of 
these simplifications holds true (Boltzmann equation with Maxwellian 
diffuse boundary condition). 

Remark 47 (Practical verification of the key conditions) . Often the 



Ilj's are nested projectors, in the sense that nj+illj 
becomes 



(Id - Hfc)' ■ (fi^7) >KJ|(id -nfc+i)^?l 



n^+i. Then flTOjl 



9 e Uk{X). 



{dt)t=o{9t - ^kgt] 



\gt - UkOtl 



t=o 



So the recipe is as follows: (a) Take g G Ilk{X), let it evolve according 
to dtg + Bg = 0; (b) compute dtijlkg) at t = 0; (c) check that \\Bg + 
dt(Jlkg)\\ controls \\g — Ilk+ig\\^~^'^ for any e > 0. 

Remark 48 (Connection with earlier works). To make the connec- 
tion with the method used in |15j . note that ii g = Hkg at t = 0, then, 
since || ■ || is Hilbertian, 

- f_ 
~ dt^ 

In view of this remark. Assumption 114.81 can be understood as a very 
abstract reformulation of the property of "instability of hydrodynamic 
description" introduced in jl5j . 

Now comes the main nonlinear result in this memoir: 

Theorem 49. Let Assu'mptions[I\to{^be satisfied. Then, for any 
f3 > there is a constant Cp, only depending on the constants appearing 
in these assumptions, on [3 and on an upper hound on S{fo) — £{foo), 
such that 

Vt>0, S{f{t))-S{U<Cpt~^. 
As a consequence, for all s > 0, 

||/(t)-/oo||. = 0(t-~). 

14.2. Method of proof. To estimate the speed of approach to 
equilibrium, the first natural thing to do is to consider the rate of decay 
of the Lyapunov functional £. From the assumptions of Theorem! 
if e > is small enough then 

(14.10) ±[£[f) - £{f^)] = -V{f) < -K,[£{f) - £{Tl,f)\'^\ 
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(I have omitted the exphcit dependence of font.) But the differential 
inequahty (jl4.1(J|) cannot in general be closed, since S{f) — S(Jlif) 
might be much smaller than S{f) — £{foo)- It may even be the case 
that / = Hi/, yet / 7^ /oo (the dissipation vanishes). So this strategy 
seems to be doomed. 

In |13|, I15j we solved this difficulty by coupling the differential 
inequality ()14.10|) with some second-order differential inequalities in- 
volving other functionals. Here on the contrary, I shall modify the 
functional £ by adding some "lower-order" terms. So the proofs in the 
present paper are based on the following auxiliary functional: 

(14.11) L(/)= [^(/)-£(/oo)]+^a,((Id-n,)/,(Id-n,)^-(i?/)), 

i=i 

where (■, ■) denotes the scalar product in and aj > (1 < j < J— 1) 
are carefully chosen small numbers, depending on smoothness bounds 
on /, and also on upper and lower bounds on £{f) — £{f^). 

The coefficients aj will be chosen in such a way that L(/) is always 
comparable to £{f) — £{foo)', still the time-derivatives of these two 
quantities will be very different, and it will be possible to close the 
differential inequalities defined in terms of L. 

When the value of £{f) — £{foo) has substantially decreased, then 
the expression of L should be re-evaluated (the coefficients aj should be 
updated), so L in itself does not really define a Lyapunov functional. 
But it will act just the same: On any time-interval where £{f) —£{foo) 
is controlled from above and below, one can choose the coefficients aj 
in such a way that {d/dt)L{f) < —KL{fY~^^, for any fixed 5 > 0. 
This will be sufficient to control the rate of decay of L to 0, and as a 
consequence the rate of decay of £ to its minimum value. 

Complete proofs will be given in the next section. It is clear that 
they enjoy some flexibility and can be slightly modified or adapted in 
case of need. 

15. Proof of the Main Theorem 

Theorem 1^ will be obtained as a consequence of the following more 
precise result: 

Theorem 50. Let AssumptionsU\ to El be satisfied, and let E > 
be such that 

(15.1) ^<£{f)-£{foo)<E. 
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Let further 

(15.2) i(/) = [£:(/)-£:(/oo)]+5^a,((id-n,)/,(id-n,)'^-(5/)), 

i=i 

where (aj)i<j<j-i positive numbers; let ao = 1. Then, 

(i) For any e G (0, 1), there is a constant K > 0, depending only 
on e and on the constants appearing in AssumptionsU\to{^ (but not on 
E) such that if aj < KE^ for all j , then 

v/ex, |<£(/)<^. 

(a) There are absolute constants eQ,k > 0, and there are constants 
K,K' > 0, depending only on e, on an upper bound on S{f) — £{foo) 
and on the constants appearing in Assumptions Ql to such that, if 
Q < e < Sq and 

Oj+i < ttj- < Ka)t\ E^\ 

aj 

for all j E {0, . . . , J — 2}, then 

V/ e X, L'if) ■ (Cf -Bf)<- aj^, K' E'+\ 

Remark 51. Lemma lA.171 in Appendix IA.22I shows that Condi- 
tions (i) and (ii) can be fulfilled with aj_i > KiE^^, where £ only 
depends on J and k. 

Remark 52. In concrete situations, the explicit form of L might 
be extremely complicated. In the case of the Boltzmann equation, to 
be considered later on, the formula for L requires eight lines of display. 

Before starting the proof of Theorem 1501 let me make some remarks 
to facilitate its reading. First of all, when uniform bounds in the 
spaces are taken for granted, a bound from above by, say, ||/ — /oo||s is 
better if the exponent a is higher^ this is somewhat contrary to what 
one is used to when working on smoothness a priori estimates. 

In all the sequel the exponents s, s' and the constants C, C", K, 
K', etc. may change from one formula to the other. These quantities 
can all be computed in terms of an upper bound on S{fo) — S{foo), 
the exponents and constants appearing in Assumptions Q to El (and 
for given e, they only involve a finite number of these constants and 
exponents). As a general rule, the symbols C, C, etc. will stand for 
constants which should be taken large enough, while the symbols K, 
K', etc. will stand for positive constants which should be taken small 
enough. 
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Finally, I shall frequently use the following fact: // < Cs' for 
all s' > 0, then for any s and any 6 there exists a constant C , only 
depending on Cg' for some s' large enough, such that 

(15.3) 11^711. <C||^7||^r^ 

To see this, it suffices to use fll4.2|) with sq = 0, si = s/6, 9 = S. In 
other words, it is always possible to replace the norm in some X'^ by 
the norm in any other X^, up to a arbitrarily small deterioration of the 
exponents. 

Proof of Theorem EOl To prove (i), it is sufficient to show that 
there exists C such that 

l-£ 



(15.4) ((Id - n,)/, (Id - n,)'^ ■ (Bf)) < c [£U) - £Uoo)] 

for all / G X. Indeed, it will follow from (jl5.ip that 

(Id - n,)/, (Id - n,)} . {Bf))\ < ^ [SU) - ^(/oo)]; 

then if aj < KE^ , the definition of C (formula ()15.2|) ) will imply 

(1 - TJKC) [£{f) - £{f^)] < L(/) < (1 + rjKC) [£{f) - £{f^)]. 

Then the conclusion will be obtained by choosing, say, K = 1/ (2^~^^ JC) . 
(Here C is the same constant as in p5.4j) .) 

To prove ()15.4j) . I shall first apply the Cauchy-Schwarz inequality, 
and bound separately ||(Id — nj)/|| and ||(Id — 11^)^ ■ {Bf)\\. 

Bound on ||(Id -Ilj)f\\: 

By Assumption Efii), / - lijf = (/ - /oo) - (n^/ - lijfoo), so 

||/-n,/|| < ||/-/^|| + ||n,-/-n,/oo||. 

By Assumption ini^iii) and the convexity of Y , Ilj is Lipschitz X^ —>■ X^ 
for some s large enough; so 

||n,/-n,/oo|| <c||/-/ooL. 

Both / and /oo belong to Y, so by Assumption El they are bounded in 
X^ for all s', and we can apply the interpolation inequality ()15.3|) : 

||/-/oo||.<C7||/-/^||i-i. 

Then by Assumption [Z(i)-(ii), 

||/-/ooir^^<c[£:(/)-£:(/oo)]^-t. 

All in all, 

(15.5) ||/-n,/||<c[£(/)-^(/oo)]^-i 
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Bound on ||(Id - n^)'^ ■ (5/) || ; 

By Assumption infiii), there are constants C and s such that 

||(id-n,)'r W)ll<c||i?/|U 

By Assumption m^i), Bf is bounded in all spaces X^' , so by interpola- 
tion, 

\\Bf\U<C\\Bfr-i. 

It follows from Assumption Efii) and the convexity of Y that B is 
Lipschitz —>■ X^ on Y; in view of Assumption EJ (-B/oo = 0), this 
leads to 




The end of the estimate is just as before: 
All in all, 

\\{ld-U,)y{Bf)\\<C[S{f)-S{U]'^-'^. 
This combined with (jl5.5|) establishes (jl5.4j) . 

Now we turn to the proof of (ii), which is considerably more tricky. 

Let 

(15.6) V{f):=-C{f)-{Cf-Bf). 
The argument will be divided in three steps. 

Step 1: The estimates in this step are mainly based on regularity 
assumptions. 

By direct computation, 

j-i 

-v{f) = -v{f) + J2 - n,)'^ ■ {Cf - Bf), (Id - n,)'^ • (Bf))) 

i=i 
j-i 

+ ^a,((id - n,)/, (Id - n,)'; ■ (c/ - 5/, 5/)) 
j-i 

+ 5^a,((id - n,)/, (Id - n,)'^ . (5} ■ {Cf - Bf))) 
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J-1 



-D(/)-^a,||(Id-n,)^-(S/)|P 
j-i 

+ 5^a,((id - n,.)'^ ■ (c/), (Id - n,)'^ ■ (Bf) 

J-i 

+ ^ a,((id - n,)/, (Id - n,)'; ■ (c/ - 5/, s/) 



j-i 

+ J]a,((id - n,)/, (Id - n^o; • (5; • {Cf - Bf)) 

Then by Cauchy-Schwarz inequality, 
(15.7) 

-V{f) < -V{f) - 5^a,||(Id - n,); • {Bf)\\' 



j-i 



+ 5^ a, II (Id - n,)'^ ■ (c/) II II (Id - n,)'^ ■ (5/) II 
+ J] a, II (Id - n,)/|| II (Id - n,)'; ■ (Cf - Bf, Bf)\\ 

J-1 

+ 5^ a, II (Id - n,)/|| II (Id - n^o; ■ {B'f . {Cf - Bf)) I 



By applying the inequality ab < (a^ + &^)/2, with a = \\ (Id — 11^)^ ■ 
(5/) II and b = ||(Id — 11^)^ ■ (C/)||, we see that the second and third 
terms in the right-hand side of (jl5.7|) can be bounded by 



j-i , J-i 



(15.8) -l^a,.||(Id -n,)'^-(5/)||%i5^a,||(Id -n,)'^-(C/)| 



Then we apply the Hilbertian inequality 

_,i„f<_M£ + ll,_,l 
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with a = (Id -n^YfiBf) and b = (Id -Ujy^^fiBUjf), to bound (fT^ 

by 

j=i j=i 

1 -^"^ 

+ 2E«^-|Kid-n.)/-(c/)|r. 

It follows, after plugging these bounds back in ()15.7j) . that 
(15.9) 

^ j-i j-i 

< _ - 5^ a, II (Id - n,);,^,^ ■ (i?n,/)||' + J2a, (/?)„ 

j=i j=i 

where 
(15.10) 

(Rh ■■= I II (M - n,)ii^.; ■ (i?n,/) - (Id - n,)'^ ■ {Bf)f 

+ ^||(id-n,)'r(c/)ir 

+ ||(id -n,)/||(||(id -n,)';-(c/-5/,5/)|| 

+ ||(id-n,)}.(i?}.(c/-i?/))||). 

Now I shall estimate the various terms in (|15.1(jp one after the other. 

First line of (I15.1()|l : 

First, 
(15.11) 

||(id-n,)[,^^-(i?n,/)-(id-n,)}-(i?/)|| < ||(id-n,)i,^^.(i?/-i?n,/)|| 

+ ||[(n,)n,/-(n,)^]- 

By Assumption ini^iii), 

II (Id - U^Yn^j . {Bf - BU,f)\\ < C\\Bf - BH.fl. 

(Here I use the fact that Ilj{X) C Y.) Also, from Assumptions El and |3 
(and again Ilj{X) C Y), Bf and BUjf are bounded in all spaces X'^' , 
so by interpolation 

\\Bf - BU.fl < C\\Bf - BU^f\\'~t 
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As a consequence of Assumption IHii) and the convexity of Y, B is 
Lipschitz continuous —>■ on Y, so 

115/ - i?n,/f -f <c\\f- n,/||r^ <c\\f- n,/f-, 

where the last inequahty is obtained again from interpolation. This 
provides a bound for the first term on the right-hand side of ()15.1H1 

Next, as a consequence of Assumption IHl^iii), (11^)' is Lipschitz con- 
tinuous on Y, in the sense that for all f,g&Y, 

||[(n,)/-(n,);]-/^||<c||/-^?|M|/.||.. 

Combining this with Assumption EJ we find 

||[(n,)'n,/ - (n,)/] ■ iBf)\\ < c\\u,f - f\u \\Bf\u 

<C'\\f-U,fl<C"\\f-U,f\\'-^. 

This takes care of the second term on the right-hand side of ()15.1HI . 
The conclusion is that the first line of ()15.10p is bounded by 0{\\f — 
n,/||i--), for any £G (0,1). 

Second line of fllS.lOj) .- 

First, by Assumption El^iii), 

||(id-n,y;-(c/)||<c||c/||,. 

By Assumption m^i), Cf is bounded in all spaces X^' , so by interpola- 
tion: 

\\cf\u<c\\cf\r'^. 

By Assumption infi), CUjf = 0; and by Assumption |3fiii), C is Lipschitz 
X' X^ on Y] so 

lic/f-f = lie/ - cn,/||i-i < c||/ - n,/||^^ < c"||/ - n./f-. 

The conclusion is that the second line of p5.10|) can be bounded just 
as the first line, by 0(||/ - YijfW^-'), for any e e (0, 1). 

Third and fourth lines of ()15.10|) .■ 

By Assumption infiii), 

||(id-n,)';.(c/-i?/,i?/)|| < c'||c/-5/|U|s/||. < ^'(iis/ii^+lic/ 

The second term ||C/||^ can be bounded by 0(||/ — flfc/jp"^^), as we 
already saw; by taking k = J we get a bound like 0(||/ — /oojP"^)- 
As for the first term ||i?/||^, we saw before that it is also bounded 
like 0(11/ — /oo||^~^)- In the sequel I shall only keep the worse bound 

0(11/ -/ooll^"^). 
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Next, by Assumption IHfiii) again, 

||(Id - n,)'^ ■ ■ {Cf - Bf))\\ < C\\B'f ■ {Cf - Bf)\U. 

From Assumption El^ii), 

\\B'f ■ {Cf - Bf)\U < C\\Cf - BfW,, < C{\\Cf\U^ + 115/11,,), 

and as before this can be controlled by 0{\\f — /oo||^~^)- 

The conclusion is that the third and fourth lines of p5.10|) can be 
bounded by C\\f - U,f\\ \\f - for any e e (0, 1). 

Gathering all these estimates and replacing e by e/2, we deduce 
that the expression in p5.1U|) can be bounded as follows: 

(R), < c (11/ - n,/f + 11/ - n,/|| 11/ - /oof -f ). 

As we already saw before, 

||/-n,/ir-f <c||/-/oo|^-^ 

so actually 

(15.12) {R)^<c\\f-U^f\\\\f-U\'-^. 
The temporary conclusion is that 

(15.13) -V{f) < -V{f) - -$^a, II (Id -Yi.y^^f ■ (i?n,/)|f 



4 



+ (Esii/-n,/ii)ii/-/ooir- 



Step 2: This step uses Assumption |H1 crucially. 
By triangle inequality, 

||/-n,/|| = ||no/-n,/|| < Iin./-n,+i/||; 

0<fc<j-l 

so 

a,||/-n,/||< ^^■( E linfc/-n,+,/|f 

l<j<J-l i<j<J-i 0<fc<j-l 

= E ( E %)linfc/-n,+i/|| 

0<fc<J-2 fc+l<i<J-l 
0<fe<J-2 
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where the last inequahty follows from the fact that the sequence {cLj)o<j<j-i 
is nonincreasing. 

Renaming /c as j , plugging this inequality back in p5.13|) , we arrive 

at 

(15.14) 

< _ 1 ^ a,||(id -n,)i,^,^-(5n,/) 

0<i<J-2 

+ c «.+i||n,-/-n,+i/||||/-/oo||i- 

0<j<J-2 

Since a/j < 1, we can write 



(15.15) 



> 



2 2(J-1) 



1 



k=l 
J-1 



On the other hand, k > j =^ aj > ak, so 



(15.16) 



j-i j-i 



5^a,||(id -n,)^^,^-(5n,/) 

i=i 



J- 



TEE«i(id-n.)n,/-(i^n,/) 



j=i k=i 
j-i j-i 



i=i fc=i 



From (Uni, (fTK7[H|l and (imTH) . 



-W) < 



— ^ (^(/) + Ell (Id - n,)^,, ■ (5n,/) 

^ ^ fc=l 7 = 1 



J-2 



c5^a,+i||n,/-n,+i/||||/-/ooir- 
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At this point we can apply Assumption IHfii)-(iii) and Assumption [3 

and we get constants C such that 



:i5.i7) -v{f) < - K([S{f) - ^(ni/)]i+i + 11/ - njii 

j-i 

-KY,(^kmkf-iik+if\\ 



\2+e 

"fc ||J-J-fcJ — J-J-fc+lJ I 

k=l 
J-2 

+ c^a,+ii|n,/-n,wil ll/-/ooir-^ 

i=o 

By applying Young's inequality, in the form 

aXY^-' <b ' ' ' ^ ^ 



2+^ \b^j m) ' 

with a = aj +i X = WHjf - Uj+if\\, F = ||/ - /^||, 6 = Kaj in the 
last line of (jl5.17p . and get 

(15.18) -v{f) < - K[S{f) - ^(ni/)]i+f - K\\f - ni/|p+^ 

J-1 



'^1 + 1 \ ,, „ „ ,, (l-e)(2+e) 



, a, 

0<j<J-2 \ 

Since a^+i < 1, we can bound trivially a^^f by Oj+i- Moreover, for 
£ < £o small enough, we have {1 — e){2 + e) / {1 + e) >2 — ie and, so 

Taking into account once again the fact that aj > aj_i for all j, (jl5.18j) 
implies, with the convention Ho/ = f, Qq = 1, 

(15.19) -P(/) < - K[Sif) - £:(ni/)]^+f 

\0<J<J-1 / 

+ c sup f^y"||/-/oo|^-^^ 

Next, 

ll/-/ooir+^ = l|no/-n,/|p+^<c ll^,/-^,+,/f+^ 

o<i<J-i 
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SO from (jl5.19j) we deduce 

(15.20) - v{f) < -K[S{f) - £:(ni/)]i+f -Kaj^^Wf - 



+ C sup f^y"||/-/oof-^^ 



Step 3: Now a few complications will arise because we only have a 
control from below of S{f) — £{foo) in terms of ||/ — /oo||; so the fact 
that £{f) — S{foo) is of order E does not imply any lower bound on 
11/ /oo||, and then ||/ — /oo||^~^^ might be much, much higher than 
11/ — /oolP"*""^- To solve this difficulty, a little additional detour will be 
useful. 

From Assumption ini^ii)-(iii) and interpolation, 

(15.21) \\uj - = iinj - njoof <c\\f- /ooir+^; 

on the other hand, 

(15.22) 11/ - u\'-'^ < c{\\f - ujf-'^ + iinj - /oof-'O- 

By using 1)15.211) and p5.22|) in p5.19p . and replacing the exponent 
1 + e/4 by the worse exponent 1 + 2e (which is allowed since S{f) — 
^(Hi/) < S{f) — £{foo) is uniformly bounded), we obtain 

(15.23) - VU) < -K[£{f) - ^(ni/)]i+2e _ Kaj.^n.f - /^of 

+ C5T^ (11/ - njip-^^ + pi/ - /oof -'^), 

where 

:= max — — . 

o<i<J-i aj 

Then from Assumption (I7j)(i)-(ii) (both the upper and the lower 
bounds are used in (ii)), 

(15.24) 

- VU) < -K[£if) - ^(ni/)]i+2e _ Kaj^^[£iU,f) - £iU]'+'^ 

+ C6Th [g(f) - £(ni/)]i-3e + C6Th [g(uj) - £iU]'-''- 

Let us distinguish two cases: 
First case: £{UJ) - £{f^) < £{f) - £{Tl^f). 
Then 

£{f)-£{u<mf)-mif)i 
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and in particular 

(15.25) S{f)-S{UJ)>^. 

In that case we throw away the second negative term in ()15.24j) . and 
bound the last term by the but-to-last one: 

-V{f) < -K[£{f) - £{Ii^f)Y+'^ + C5^^ [£{f) - £{Iiif)Y-'' 
(15.26) 

= -k(^- ^ [£{f) - £{Ti,f)r'^ [£{f) - £{Tiif)Y^''. 

If 

(15.27) <^[m-£{^rf)f 

(where K and C are the same constants as in ()15.26p ). then p5.26|) 
can be bounded above by -K'[£{f) - £{Iiif)Y+^' , and by (fTK:^ this 
can also be bounded above by —K"E^^'^^. 

Finally, in view of (|15.25|) again, (jl5.27p is satisfied if 

(15.28) 6^ < K'E^', 

where K' = 4~^'^K/{2C). Since e <1 and E is uniformly bounded, a 
sufficient condition for (115.281) to hold \s 5 < K"E^'^'. 

Second case: £{Ilif) — £{foo) > ^(/) ~ ^(J^if)- In that case 

(15.29) £(ni/)-£(/^)>|, 
and we retain from ()15.24|) that 

-Vif) < -Kaj.,[£{UJ)-£iU]'^'^ + C6^ [£iUJ) -£iUY~''- 
and by a reasoning similar as the one above, this is bounded above by 

aj^, [£{UJ) - £{U 



-^aj^, [£iji,f) - £{UY^'' 
as soon as 



5— < K'aj^iE^'. 
This condition is fulfilled as soon as 

6 < R'ayjiE^"', 
and, a fortiori (since aj_i < 1) if 

S < K"'a)t^^E^^'. 



104 

In both cases, we have concluded that if e < Eq and 
then 

-V{f) = L'(/) ■ {Cf - Bf) < -K'aj^, E'+''. 

Up to the replacement of e by e/2 and Eq by eo/2, this is exactly the 
desired conclusion. □ 

Proof of Theorem HHl Let E be such that S{fo) - £{foo) < E. 
Since S{f{t)) is a nonincreasing function of t, 

Vt>0, £{f{t))-£{U<E. 

Let now 5 > 0, and E E {0,E]. Let [to, to + T] be the time-interval 
where 

^<£{fit))-£iU<E; 

this time-interval is well-defined (at least if T is a priori allowed to be 
infinite) since £{f{t)) — £{foo) is a continuous nonincreasing function. 
The goal is to show that if e is small enough, then 

(15.30) T < CE-^', 

where A only depends on J, and C may depend on e but not on E. 
When (I15.30|l is proven, it will follow from a classical argument that 

(15.31) £{f{t)) - £{U = 0(ri/«^+i)^)). 

Indeed, let Eq := £{fo) —£{foo)', then £{f(t)) —£{foo) will be bounded 
by Eo/2'^~^^ after a time 

/ T~\ \ / T— 1 \ —As / T~\ \ —As 



^0 • 



<C{J2 2^^')Eq^' < C'2^"''Eo 
j=0 

So £{f{t)) - £{foo) = 0(2-('^+i)) after a time proportional to 2^"^ 
and 1)15.311) follows immediately. 

Then £{f{t)) — £{foo) = 0{t~°°) since e can be chosen arbitrarily 
small and A does not depend on e. From Assumption [7|^i)-(ii), 

\\f{t)-U\<C[£if{t))-E{U)Y/' 

(here 1/3 could be 1/2 - e), so \\f{t) - foc\\ = 0{t~°°) also. Finally, 
since f{t) is bounded in all spaces by Assumption 01 it follows by 
interpolation that \\f{t) — foo\\s = 0{t^°°) for any s > 0. 
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So it all amounts to proving Let K,K',k,eo be provided 

by Theorem ^01 (There is no loss of generality in taking the constants 
K appearing in (i) and (ii) to be equal.) Let then ei, Ki,i he provided 
by Lemma FA.I?! For any e < min(£:o,£^i) and any t G [to, to + T] we 
have 

I < Hfit)) < ^; I [mm < -K'aj^iE'^^ < -K'K,E'+('^'^\ 
So 

E > L(/(to)) - L(/(to + T)) > T K'K^ 

and then ()15.30|) follows with X = £+1 (which eventually depends only 
on J). □ 



16. Compressible Navier— Stokes system 



In this section I start to show how to apply Theorem 09] on "con- 
crete" examples. 

The compressible Navier-Stokes equations take the general form 



(16.1) 



dtp + V- (pu) = 

dt{pu) + V ■ (pn ® n) + Vp = V ■ r 

dt{pe) + V ■ (pen + pu) = V • (rn) - 



V ■ q 



where p is the density, u (vector- valued) is the velocity, e is the energy, 
q (vector- valued) is the heat flux, and r (matrix-valued) is the viscous 
stress. In the case of perfect gases in dimension N, it is natural to use 
the following constitutive laws: 



(16.2) 



p = pT 



LA. + _r 
2 2 



T = 2p.{Vu} 



= ——1^ VT, 



where T is the temperature, p is the viscosity, k is the heat conductivity. 



wnere I is tne temperature, p is tne viscosity, k is tne neat coi 
and {Vn} (matrix-valued) is the traceless symmetric strain: 



u 



N 
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and 6ij = li=j- Then flKi.lll takes the form 
(16.3) 

'dtp + V-{pu) = 0; 

dt{pu) + V ■ {pu®u) + V{pT) =2pV ■ {Vn}; 

< 

(p^ + - pT j + V . (p^ n + [-^ j put) 

= 2/iV-(u{Vn}) + f kAT. 

(Note that 

V ■ {Vn} = yuAn + p(^- -u, 

so in the case considered here, the second Lame coefficient is negative 
and equal to —{2/N)p, which is the borderhne case.) 

To avoid discussing boundary conditions I shall only consider the 
case when x varies in the torus T^. (Later on, for the Boltzmann 
equation we'll come to grips with boundary conditions a bit more.) 

There are N + 2 conservation laws for (jl6.1|) : total mass, total 
momentum (A^ scalar quantities) and total kinetic energy. Without 
loss of generality I shall assume 

f f f N f N 

(16.4) yp=l; J pu = 0; JpL^ + -JpT = -. 

There is an obvious stationary state: (p, u, T) = (1, 0, 1). The goal 
of this section is the following conditional nonlinear stability result. 
The notation stands for the usual space of functions whose deriva- 
tives up to order k are bounded. 

Theorem 53 (Conditional convergence for compressible Navier-Stokes) . 
Let t — > fit) = {p{t),u(t),T{t)) be a C°° solution of p6.3|) . satisfying 
the uniform bounds 
(16.5) 

fVA; G N sup,>o(l|p(t)llc^ + Mt)\\c^ + WmWc^) < +oo; 
[Vt > 0, pit) >pm>0; Tit) > T,„ > 0. 

Then ||/(t) - (1,0,1)11^^ = 0(r°°) for all k. 

Proof of Theorem ESI Let us check that all assumptions of The- 
orem 13^ are satisfied. Assumption ^ is satisfied with, say, = 
H'iT^; R X X M), where stands for the usual L^-Sobolev space 
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of functions with s derivatives in L^. To fulfill Assumption |21 define 
Cs :=sup{||/(t)||^^.; t>0} and let 

X = r := {/gC°°(T^;MxM^xM); Vs, ||/||. < C,; p > T > T„}. 

(Note that necessarily Pm < 1, < 1.) 

To check Assumption |21 rewrite (jl6.3|) in the nonconservative form 



'{dt + u-V)p + p{V -u) =0 



(16.6) < 



{dt + u-V)u + VT + T = ^V - {Vn} 

VP/ P 

{dt + u-V)T + ^T{V-u) = ^^^\{Vu}\' + -^AT, 

and define 
(16.7) 

Bf= (u-Vp + p{V-u), M-VM + VT + rV(logp), M-Vr+^T(V-n)); 



:i6.8) 



C/=(0, ^V-{V«}, ^^\{Vu}\' + -At). 
\ p JM p P y 



2p 4 p 

— V ■ { Vm|, — - 

Then dSj) obviously holds true. 

Assumption m is satisfied with /oo = (1,0, 1). 

As usual in the theory of viscous compressible flows, an important 
difficulty to overcome is the fact that diffusion does not act on the p 
variable. So let IIi = 11 be defined by 

n(p,M,r) = (p,o,i). 

Assumption IHl is obviously satisfied with this choice of nonlinear pro- 
jection. 

Next, let S be the negative of the usual entropy for perfect fluids: 
£{p,u,T) = J plogp ~ yJ 
Taking into account p6.4|) . 

£ip,u,T)- £(1,0,1) = J plogp + J p\^ + ^J p(T-logT-l); 

£{p,u,T)-£{u{p,u,T)) = J p^-Y + Y y p(r-logr-l). 
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Thanks to the uniform bounds from above and below on p and T, 
£{f) — £{foo) controls ||/ — /oo|P from above, and — £^(11/) controls 
11/ — n/lp from above and below; so Assumption [7| is satisfied. 

It only remains to check Assumption |H1 By a classical computation, 
for any f EY, 

<-K (j |{Vi.}P + Jp\VtA, 

where the last inequality follows again from the lower bound on T and 
the upper bound on p. 

By Poincare inequality, J p|VTp controls J p{T — {T)pY, where 
{T)p = J pT is the average of T with respect to p. In turn, this 
controls ||T - if - 2{{T)p - if. Since (T)^ - 1 = (-1/7V) / p\u\\ we 
conclude that there are positive constants K and C such that 

j p|VT|2 > A'llT- If -C||nf 

for all f eY . On the other hand, by |151 Proposition 11], 

j |{Vm}P > K'\\u\\\ 

All in all, there is a positive constant K such that 

/ |{Vn}r + j p\Vn> K{\\T - If + ||nf ) = K\\f - B/f , 

so Assumption |Hfi) holds true. 

By another classical computation, £\f) ■ (Bf) = 0, so Assump- 
tion |Hl^ii) also holds true. 

On the range of B, the functional derivative B' vanishes (because 
dtp = when u = 0), and B{p,u,T) = (0, -V logp, 0). Then 

(Id - Uy^f ■ (BUf) = BUf = (0, -Vlogp,0). 

Thus 

|(Id-B)'n^-(SB/) 

which under our assumptions controls J | Vpp, and then by Poincare in- 
equality also ||p— if = ||B/— /oof . This establishes AssumptionlSfiii), 
and then the conclusion of the theorem follows from Theorem 0^1 D 



= / |V(logp)p, 
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17. Weakly self-consistent Vlasov—Fokker— Planck equation 

One of the final goals of the theory which I have been trying to 
start in this memoir is the convergence to equilibrium for the nonlinear 
Vlasov-Poisson-Fokker-Planck equation with an external confinement. 
This kinetic model, of great importance in plasma physics, describes 
the evolution of a cloud of charged particles undergoing deterministic 
and random (white noise) forcing, friction, and infiuencing each other 
by means of Coulomb interaction. 

Besides the fact that the regularity theory of the Vlasov-Poisson- 
Fokker-Planck equation is still at an early stage (to say the least), 
one meets serious difficulties when trying to apply Theorem EUl to this 
model, in particular because the problem is set in the whole space. 
So for the moment I shall be content to treat a simpler baby problem 
where (a) the confining potential is replaced by a periodic boundary 
condition; (b) the Coulomb interaction potential is replaced by a small 
and smooth potential. The smallness assumption is not only a techni- 
cal simplification: It will prevent phase transition and guarantee the 
uniqueness of equilibrium state. 

Even with these simplifications, the problem of convergence to equi- 
librium is nontrivial because the model is nonlinear and the diffusion 
only acts on the velocity variable. This will be a perfect example of 
application of Theorem 

Here the unknown / = f{t,x,v) is a time-dependent probability 
density in phase space {x G stands for position and v G for 
velocity). The equation reads 



'dj + v VJ + F[f]{t,x) ■ V,J = A J + V, ■ ifv) 



F[f]it,x) = - / VWix-y)fit,y,w)dwdy. 



:i7.1] 



Here W G C°°(T^) is even {W{-z) = W{z)), and without loss of 
generality J W = 0. As we shall see later, if W is small enough in a 
suitable sense then the unique equilibrium for (|17.1|) is the Maxwellian 
with constant density: 

U{x,v)=M{v) 



(27r)^/2' 

Since the total mass J f{t,x,v) dvdx is preserved with time, there 
is an a priori estimate on the force, like H-FHc* ^ C'fcll W^llc'^j so there 
is no real difficulty in adapting the proofs of regularity for the linear 
kinetic Fokker-Planck equation (see Appendix lA. 2 In this way one 



no 

can establish the existence and uniqueness of a solution as soon as, 
say, the initial datum has finite moments of arbitrary order; and this 
solution will be smooth for positive times. 

The goal of this section is to establish the following convergence 
result: 

Theorem 54 (Large-time behavior of the weakly self-consistent 
Vlasov-Fokker-Planck equation). Let W e C~(T^) satisfy f W = 
0. Let fo = fQ{x,v) be a probability density on x M^, such that 
f fQ{x,v)\v\'' dv dx < +00 for all k ^ N, and let f = f{t,x,v) be the 
unique smooth solution of ()17.H) . Let 6 be so small that 

. 1 
< -. 



Ifma.x\W\ < 6 then 



2 2 



||/(t,-)-M|Ui = 0(r~). 



Remark 55. It is not hard to show that the conclusion of Theo- 
rem 1^ does not hold true without any size condition on W, since in 
general (jl7.ip can admit several stationary states. In the proof of The- 
orem I shall show that there is only one stationary state as soon as 
max \ W\ < 1; I don't know how good this bound is. The assumptions 
of the theorem are satisfied with 6 = 0.38, which does leaves some 
margin of improvement. 

Proof of Theorem [HH The first step consists in establishing 
uniform regularity estimates; I shall only sketch them very briefiy. 

First, one establishes differential inequalities on the "regularized" 
moments Mk{t) = J f{t,x,v){l + \v\'^)'^/'^ dv dx: 

^ < CM,., - KM,, 
dt 

where C and K are positive constants. (Here the fact that the position 
space is induces a considerable simplifcation.) Then one deduces 
easily that each moment J f\v\^ remains bounded uniformly in time. 

Next, by adapting the arguments in Appendix IA.2H one can prove 
uniform Sobolev estimates of the form 

WkeN, Vto>0, sup||/(t,-)||Hfc < +00. 

t>to 



These bounds, combined with the moment estimates, imply the bound- 
edness of the solution /(t, ■) in all spaces X'^, where X'^ is defined for 
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k eNhj 

(17.2) 11/11^.= f\'^''^7fix,v)\\l + \v\')'dvdx. 

and is defined by interpolation for noninteger s. It is easy to clieck 
tfiat tliese spaces satisfy Assumption [TJ 

Finally, classical methods based on the maximum principle (as 
in |121 Section 10]) suffice to show that 

f{t,x,v) > Ke-'^l^'' 

uniformly in t > to, provided that K is small enough and a is large 
enough. (Here again, the assumption that the position space is 
simplifies things quite a bit by allowing x to be treated as a parameter.) 
It follows that p(t,x) = J f(t,x,v)dv is bounded below by a uniform 
positive constant for t > to > 0. 

Up to changing the origin of time, we can now assume that / is 
uniformly bounded in all spaces X'^ and that p satisfies a uniform lower 
bound. This determines a workspace 

X = Y:=[f; Ws\\f\\xs<Cs; p>Pm>o], 

as in Assumption |21 
Then we define 

Bf = v- VJ + F[f] ■ VJ; Cf = AJ + V, ■ (fv); 

/oo = M{v)- ni(/) = n(/) =pM, p = j fdv. 

Assumptions El IH El and El are readily checked. 

Next let the free energy functional £ be defined by 

^(/) = y" flogfdvdx + J f^-^dvdx + ^ J p{x) p{y)W{x-y) dx dy. 
By standard manipulations, 

(17.3) £{f)-£{m= [ /log ^ 



(17.4) £{Uf)-£{U= [ p\ogp+\ f p{x)p{y)W{x-y)dxdy. 

JfN I JfN 

By the Csiszar-Kullback-Pinsker inequality, ^(/)-^ (Hi/) > (1/2)||/- 
pM||^i; then by interpolation of l? between and (as in |151 
Lemma 10]), one deduces 

£[f) - £:(n/) > - pMlli. >K\\f- pM||-l||/ - pMii^y, 
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where 9 is arbitrarily small if k is chosen large enough. This shows that 
Assumption [7|^i) is satisfied. 

On the other hand, since J TV = 0, 



£{Ilf) -£{f^)= p\ogp+- [p{x) - 1] [p{y) - 1] W{x - y) dx dy 



1„ 
> - P 



1 



(max \W\) \\p 



|2 



II, rl. 



By assumption, max \ W\ < 1; so there is a constant K > Q such that 

8{I[f)-S{U>K\\p-l\\l,. 

By interpolation again, this can be controlled from below by l||^t^ 
for arbitrarily small e, and the left inequality in Assumption El^ii) is 
satisfied. (This is the first time that we use the smallness assumption 
on W .) The right inequality in Assumption [7(ii) is easy. 
By classical computations (see e.g. jl3l Section 2]), 

2 



f 



V„ log 



/ 



pM 



dvdx>2 I /log-^>i||/-pM||i,. 



so there is no difficulty to establish Assumption |Hl^i). Assumption IHl^ii) 
follows immediately since S'{f) ■ (Bf) = 0. So it only remains to 
estabhsh Assumption |Sfin). 

As in the example of the compressible Navier-Stokes system, the 
functional derivative 11' vanishes on the range of 11, so 

(Id - Uy^j ■ (BUf) = BUf = v V,{pM) - {VW * p) ■ V,{pM) 

= v-\/,{p + piW*p))M. 



Then 



(Id - n);,^ ■ {BUf) 



V.^{p + p{W * p)) M{v)dvdx 



> K 



^x{p + p{p*W)) 
Vp 



P 



p*W 



dx 



dx. 



where the lower bound on p was used in the last inequality. Let 



p{x) 



-{p*W){x) 



Since p is uniformly bounded from above and below, we can use a 
logarithmic Sobolev inequality with reference measure p{x) dx; so there 
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is a positive constant K such that 



P 



dx 



Vlog^ 
/i 

P 



dx 



> K p log — dx 
P 



;i7.5) 



K (^j plogp + j p{p*W)- log j e-^^^P 



By assumption, max \W\ < 6; so \W * p\ < S, and 



e"^*'' - (1 - * p) 



< 



e'^(max|Vr|)^ 



|p-i|lii < 



|P-1| 



2 - 2 

Since J {W * p) = 0, it follows by integration of this bound that 



-W*p 



< 



\P-M 



As a consequence, 



log 



-W*p \ ^ 



2„<5 



|P-1| 



From this bound and the inequality | / p(p * W)\ < 6\\p — l||^i again, 
we obtain 

plogp + [ p{p*W)- log I e-^*^>fc^ -5||p-l||i,-^ 



>'2-^ 



|P-1| 



By assumption the coefficient in front of ||p — l||^i is positive, and then 
we can use interpolation again to get 

j plogp + J p{p*W)-log J e-^*^ > K,\\p-l\\ir = K,\\pM-M\\ir. 

So Assumption IHtiii) holds. (Here again the smallness condition was 
crucially used.) Then all the assumptions of Theorem are satisfied, 
and the conclusion follows at once. □ 



114 



18. Boltzmann equation 

This last section is devoted to tlie Boltzmann equation; see |51j 
and the references therein for background and references on this model. 
I have personally devoted a considerable amount of research time on 
the problem of convergence to equilibrium for the Boltzmann equation, 
alone or in collaborations with Toscani and Desvillettes; a detailed 
account of this topic can be found in my lecture notes |50j . 

As in Section El the unknown is a time-dependent probability den- 
sity / = f(t,x,v) on the phase space. The variable x will be assumed 
to vary in a bounded A^-dimensional domain Q^, that will be either 
the torus T , or a smooth bounded connected open subset of M^. The 
equation reads 
(18.1) 

§^ + ^-v./ = g(/,/) 



Q{fJ)= / f{x,v')f{x,ii)-f{x,v)f{x,v^) B{v -v^,a)dadv^ 

. v + v^ \v — vA , v\v^ \v — vA 

V = 1 a; = a. 

2 2 ' * 2 2 

Here B is the collision kernel; for simplicity I shall restrict to the 
case B = \v — (hard spheres interaction), but the analysis works as 
soon as Assumptions (5) and (19) in |15j are satisfied, which covers all 
physically relevant cases that I know of. 

Three kinds of estimates play an important role in the modern 
theory of the Boltzmann equation: Sobolev estimates (in x and v vari- 
ables), moment estimates and positivity estimates of the form / > 
-f^oc""^'''^'™. At least in some cases, the positivity estimates follow from 
regularity estimates |37j . but I shall not address this issue here. 

To continue the discussion it is necessary to take boundary con- 
ditions into account. I shall consider five cases: (i) periodic bound- 
ary conditions; (ii) bounce-back boundary conditions; (iii) specular 
reflection in a nonaxisymmetric domain; (iv) specular reflection in a 
spherically symmetric domain; (v) Maxwellian accommodation with 
constant wall temperature. Cases (i) to (iii) were already considered 
in |15j . while cases (iv) and (v) are new and will be the occasion of in- 
teresting developments. Specular reflection in a general axisymmetric 
domain (not spherically symmetric) is intermediate between cases (iii) 
and (iv) and can probably be treated as a variant, but I have not tried 
to do so. Other conditions could be treated as a variant of (v), such as 
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more general accommodation kernels, but they do not seem to cause 
any substantial additional difficulty. On the other hand, the techniques 
presented here are helpless to treat accommodation with variable wall 
temperature, for which the collision operator does not vanish; I shall 
add a few words about this issue in the end of the section. 

18.1. Periodic boundary conditions. In this subsection I shall 
consider the Boltzmann equation (jl8.Hl in the position space = 
(the A^- dimensional torus). Then there are + 2 conservation laws: 
total mass, total momentum (A^ components) and total kinetic energy. 
Without loss of generality, I shall assume 

(18.2) Jfdvdx = l] jfvdvdx = 0; J f\v\'^ dv dx = N. 
Then the equilibrium state takes the form 

f^{x,v)=M{v) 



(27r)^/2- 

Our goal is the next convergence theorem: 

Theorem 56 (Convergence for the Boltzmann equation with peri- 
odic boundary conditions). Let f be a solution of 1)18.11) in the spatial 
domain = T^, satisfying the conservation laws p8.2|) . and the uni- 
form regularity estimates 

Vs > sup ||/(t, ■)l|//''(n,xMA^) < +oo; 

t>0 

\/k>0 sup / f{t,x,v)\v\^dvdx<+oo; 
t>o J 

[Wit, X, v)eR+x a, X , fit, X, v) > Ko e"^"!''!''' . 



;i8.3) < 



Then 



Vs>0, \\fit,-)-M\\^^=0{t- 



Remark 57. This theorem is nonempty, in the sense that, when /o 
is very smooth and close to equilibrium in a suitable sense, then p8.3|) 
holds true. See the discussion in |15j for more information. 

Proof of Theorem EHl Let / satisfy the assumptions of The- 
orem Let {X^)s>o be the scale of weighted Sobolev spaces al- 
ready defined in the treatment of the Vlasov-Fokker-Planck equa- 
tion (recall equation ()17.2|) ). It follows from the assumptions that 



116 

Cs := sup^>Q II /lis is finite for all s. We shall work in the spaces 
X:=|/; 11/11. <C,; /(x, t;) > i^o e"^"'^!™ ' 



Y:=[f; \\f\\s<C',- /(x,t;)>i^^e-^°l''l'°}, 

where C^, K'^ will be determined later on. Then Assumptions d |21 
and El are obviously satisfied. 
Define 

Bf = vVJ; Cf = QifJ). 

Then Assumption |3fi) is obviously true, Assumption llfii) is satisfied 
since B is linear continuous X'^'^^ — > X*, and Assumption EJ^iii) is a 
consequence of |15| eq. (78)]. 

Assumption El holds true with /oo = M; notice that both the trans- 
port and the collision part vanish on /oo. 

Next, if / = f{x,v) is given, define 



p = I fdv; u = ^J fvdv; T=^ j f\v-u\'^dv, 



and 

p{x) e~^^^ 



\v — u{x)\'^ 



M, 



puT 



(Note that M depends on / via p,u,T.) It is easy to derive uniform 
estimates of smoothness on p, u and T in terms of the estimates on 
/; and to derive similarly strict positivity estimates on p, T: See |15| 
Proposition 7]. 

Now we can introduce the nonlinear projection operators: 

IlJ = MpuT; n2/ = M,oi; n3/ = M. 

By adjusting the constants C^, Kq, we can ensure that Ilj{X) C Y. 
Then the rest of Assumption IHl follows easily. 

The natural Lyapunov functional in the present case is of course 
Boltzmann's H functional: 

E{f) = H{f)= f f log fdvdx. 

By standard computations, taking into account 118.21 we have 

/ 



:i8.4) S{f)-S{UJ) = J /log 



2 



;i8.5) S{UJ)-S{U = Jp\ogp + j P^^ + J p{T-\ogT-l). 
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To find a lower bound on ()18.4j) . it suffices to use the Csiszar-Kullback- 
Pinsker inequality and interpolation, as we did previously for the Vlasov- 
Fokker-Planck equation (recall ()17.3p : or |151 eq. (47) i]). Upper and 
lower bounds for ()18.5|) can be obtained as we did before for the com- 
pressible Navier-Stokes equations. So Assumption [7| is satisfied. 

Now the crucial step consists in checking Assumption |S1 By a clas- 
sical computation, 



-H'U)-{Cf) = j D{f{xr))dx, 



where D{f) is Boltzmann's dissipation of information: 



B{v — v^, a) dadvdv^. 



Known entropy production estimates from |52j make it possible to es- 
timate D{f) from below by [H{f) — H^Mpur)] ■ (Such estimates 
go back to |48j : see also |50j for a detailed account on this problem.) 
Then Assumption |Sl^i) follows easily, as in |151 Corollary 5] . 

Assumption ISl^ii) is an immediate consequence of Assumption |Hl^i), 
since H'{f) ■ (Bf) = 0. 

It remains to establish Assumption IHtiii). For this we use Re- 
mark HTj According to |151 eq. (69)], if / 



M 

and / evolves according to dtf + v ■ V xf = 0, then 
d^_ 

-2 

t=0 



/ 

puT 



at time t 



dt^ 



\\f-Ml^A>K( [ |Vrp+ / \{Vu}\' 

\JjN JjN 



where, as in Section UHl 



{Vu} 



dui duj 
2 \dxi dxi 



and V • M is the divergence of u. According to |151 Section IV.2], there 
are constants Ki, K2, only depending on N , such that 



|Vrp> A'i||T-l|p-C|k||2; 



I \{Vu}\^>K2\\Vuf>K4u\\\ 



118 



This implies 



\ i WT 1 l|2 I ||„, ||2 



(id-ni^^^-(i?ni/)||;, >i^(||T-i 



u\ 



>i^'||ni/-n2/f. 

Next, according to |15^ eq. (71)], if / = Mpoi at time t = and / 
evolves according to dtf + v ■ V xf = 0, then 

^ \\f-Ml,,\\>Kl \Vp\\ 
at jjN 

Combining this with a Poincare inequality (see again |15| Section IV. 2]), 
we deduce that 

||(Id - Y[2)'n,f ■ {BYi^DWl, ^ K\\P -M\'> K'\\U2f - UH^. 

(This is in fact as in Sectional if "we set W = 0.) 

This concludes the verification of Assumption and the result 
follows by an application of Theorem UHl D 



Remark 58. A comparison with the proof of the same result in 
shows that the crucial functional inequalities are all the same; but there 
are essential simplifications in that (a) it suffices to do the computations 
for Maxwellian states ("local equilibrium" in the language of |15j ): 
and especially (b) there is no longer need for the tricky analysis of the 
system of differential inequalities. More explicitly. Sections III. 3, V 
and VI of |15j are shortcut by the use of Theorem 

18.2. Bounce-back condition. Now let fij; be a bounded smooth 
open subset of M^; up to rescaling units we may assume that \Qx\ = 1 
(the Lebesgue measure of the domain is normalized). In this subsection 
the boundary condition is of bounce-back type: 

(18.6) xedn^^ fix, v) = fix, -v). 

A consequence of fll8.6|) is that n = on dQx (the mean velocity 
vanishes on the boundary). 

Now there are only 2 conservation laws: mass and energy. So, 
without loss of generality, I shall assume 

(18.7) jfdvdx = l; j f\v\^ dv dx = N. 
The equilibrium is again the steady Maxwellian, 

fooix,v) =Miv) 



(27r)^/2' 

Here is the analogue of Theorem EHl 
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Theorem 59 (Convergence for the Boltzmann equation with bounce-back 
boundary conditions) . Let f be a solution of (jl8.1|) in a smooth bounded 
connected spatial domain fi^., satisfying bounce-back boundary condi- 
tions, the conservation laws ()18.7p . and the uniform regularity esti- 
mates p8.3|) . Then 

ys>0, ||/(t,-)-M||^, =0(t^°°). 

Proof of Theorem IH^ The proof is quite similar to the proof of 
Theorem inni however the sequence of projection operators is different: 

IlJ = Mp,T; U^f = Mp^(^ry, U3f = Mpoi; IIJ = M, 

where (T) = J pT is the average temperature. According to |15| 
eq. (70)-(71)] and a reasoning similar to the one in the proof of Theo- 
rem EHl 

2 



:i8.8) 



(Id -ni)n^^-(5ni/) 

{Id-U,)' .{BU,f) 



> A' II VT II 2; 



(Id -no;^^-(sni/)|| >i^||Vp|P, 

where V^^™u is the symmetrized gradient of u, that is 



sym 



(V 

By Poincare inequalities, 

llvriP > K\\T 



u).. = - 



1 / duj du,, 



dxj 



+ 



dxi 



(T)|p; llVpf >i^||p-l||^ 

By the classical Korn inequality, and the Poincare inequality again 
(component- wise) , 



7sym 



U 



|2 > KWVuf > K'Wuf. 



These estimates imply ||(Id - B^)^^,^ ■ {BUjf)\\^ > K\\Ujf - Uj+Jf 
for all j G {1,2,3}, so Assumption |Hl^iii) is satisfied in the end. Then 
Theorem 021 applies. □ 



18.3. Specular reflection in a nonaxisymmetric domain. In 

this subsection the bounce-back boundary condition is replaced by the 
specular reflection condition: 

X G dQx =^ f{x, v) = f{x, Rxv), RxV = V — 2{v, n)n. 

This condition is more degenerate and the shape of the domain will 
influence the form of the equilibrium. For the moment I shall assume 
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that the domain is nonaxisymmetric in dimension N = 3. The notation 
is the same as in Subsection 118.21 

Theorem 60 (Convergence for the Boltzmann equation with non- 
axisymmetric specular conditions). Let f be a solution of ()18.H) in a 
smooth bounded connected nonaxisymmetric spatial domain C M^, 
satisfying specular boundary condition, the conservation laws ()18.7|) . 
and the uniform regularity estimates ()18.3|1 . Then 

Vs>0, ||/(t,-)-M||^, =0(t-°°). 

Proof of Theorem IHUl The proof is entirely similar to the proof 
of Theorem EHl except that the condition u = on the boundary is 
replaced by the weaker condition u ■ n = 0, where n is the inner unit 
normal to Qx- Then the classical Korn inequality should be replaced by 
the Korn inequality established by Desvillettes and myself in |14j . □ 

18.4. Specular reflection in a spherically symmetric do- 
main. In this subsection is a bounded smooth connected spher- 
ically symmetric domain in M^; so, up to translation, is either a 
ball (|x| < R) or a shell (0 < r < |a;| < i?). Again I shall assume 
that ll^i:! = 1. 1 shall write = 3 to keep track of the role of the 
dimension in various formulas (certainly the analysis can be extended 
to more general domains, but one has to be careful about the meaning 
of the conservation of angular momentum). 

Now there are + 2 conservation laws: mass, kinetic energy and 
angular momentum {N scalar quantities). Without loss of generality, 
I shall assume 

(18.9) Jfdvdx = l; j f{x,v)\v\^ dv dx = N- 

j f{x,v){v Ax)dvdx = MeM.^ . 

The existence of an equilibrium is not trivial if M 7^ 0, and the equilib- 
rium does not seem to be explicit. It is a local Maxwellian with uniform 
temperature 9, but nonzero velocity Mqo and nonhomogeneous density 
Poo- The equations determining this equilibrium were studied, at the 
beginning of the nineties, by Desvillettes Here I shall suggest a 

variational approach to this problem, by means of the following lemma 
from elementary calculus of variations (the proof of which will be only 
sketched): 

Lemma 61 (stationary solutions in a spherically symmetric do- 
main). Let be a spherically symmetric domain in M.^ , N = 3, 
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\Qx\ = 1- Whenever p is a nonnegative integrable density on fix, o,nd 
m G Li(fi^.;M^), define 

f N f 1 f |mP\ 

F(p,m) = J plogp - y log (1 - — J ——] . 

Then there is a unique {poo, rrioo) G C°°{Qx] ^+ x which minimizes 
the functional F under the constraints 



(18.10) y*^"^' Jm{x)Axdx = M. 

Moreover, p is strictly positive; and there are an antisymmetric matrix 
Eoo and positive constants O^o and Z such that for all x &VLx, 



p[x) Z 
Sketch of proof of Lemma EH Write 



N J P 
then 



2 



F(p,m) = y" plogp + ^(l-^) + ^(^-log^-l) 



piogp + y ^ + 

where ^{B) = {N/2){e -logO - 1). 

By a classical computation, {p,m) i — > \m\'^/p is convex, so 6' is a 
concave function of (p, m). Moreover, 6 remains in (0, 1), and on that 
interval is a convex decreasing function of 6. It follows that "^{0) is 
a strictly convex function of {p,m). So 

/f |?77.P 
plogp + J ^ + 

is a strictly convex function of (p, m). This conclusion does not change 
if F is restricted on the domain defined by the /inear constraints (|18.10|) : 
so F has at most one minimizer. 

The Euler-Lagrange equations for the minimization of F read 

|2 



:i8.iii 



N frrii 
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where {^j)o<j<N are constants, (mj)i<j<Ar are the components of m, and 
Eijk is defined by the equations (a A 6)j = ^ Sijkajhk- These equations 
imply that m/p is an antisymmetric hnear function of x. In particular, 
the minimizer a priori lives in a finite-dimensional space. The rest of 
the lemma follows by classical arguments. □ 

The goal of the present subsection is the following result: 

Theorem 62 (Convergence for the Boltzmann equation with spher- 
ically symmetric specular conditions). Let f he a solution of (jl8.H) 
in a smooth hounded connected spherically symmetric spatial domain 
Vtx C M^, satisfying specular houndary condition, the conservation 
laws fjl8.9|) and the uniform regularity estimates fll8.3|) . Then 

V5>0, ||/(t,.)-/^||^^=0(t— ), 

where 

foo{x,v) = ^^ (27r^^)3/2' 

and the antisymmetric matrix Sqo, the positive constants Z and 6^0 are 
provided hy Lemma[U7[ 

Remark 63. The variable 6qq is the (uniform) equilibrium temper- 
ature; the velocity field in the stationary state is still rotating, and the 
density is lower near the interior of the box. 

Proof of Theorem Ell The only differences with the previously 
treated cases lie in the definition of the projection operators, and the 
verification of Assumptions [7|^ii) and|Hfiii). 

In the present case, let 

S:=(V^w), 9:={T)„ 

more explicitly, S is the average value of the antisymmetric part of 
the matrix-valued field Vm (the averaging measure is the normalized 
Lebesgue measure), while 9 is the average value of the temperature (but 
now the averaging measure has density p). 1 shall identify the matrix 
S with the velocity field x 1 — > T,x, and 6 with the constant function 
X ^ 9. Then the sequence of projection operators is as follows: 

Iiif = Mp^T] Il2f = Mp,e; U3f = Mp^e; 

^if = foo = 9oo • 

Once again the Lyapunov functional is 

H{f) = J /log/. 
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After taking into account the conservation laws ()18.9|1 . one observes 
that 



;i8.12) 
H{UJ)-H{f^ 



N 



\u\ 



- (^j Pooiogpoo + P'^(^ ~ i°g^ j 

Let again $(6*) = 9 — log 6* — 1: then by Jensen's inequality (in quanti- 
tative form), 



j p$(r)>$((r),) + Air-(r) 



||2 
p\\ 1 



where K depends on the bounds on p and T . Plugging this in ()18.12j) 
and using the same notation as in Lemma |^ one obtains the lower 
bound 



H{Yi,f) - H{U 



N 



m) - F{p 



-'OO ; 



(18.13) > K\\T - (r),f + [F{p, m) - F{p, 

The upper bound 

H{U^f) - H{f^) <C\\T- {T)X + [F(p, m) - F{p^, m^)] 

is obtained in a similar way. 

So to prove Assumption [7||^ii), it suffices to check that 



K\\Uif - /oof < F(p, m) - F(poo, moo) < C||ni/ - /< 
or, which amounts to the same. 



2. 

CO II ) 



'18.14) (p,m) - (poo, moo) < F{p,m) - F{p^,m^ 



< C 



{p,m) - {p 

OO ; mooj 



The upper bound is obvious from the definition, the bounds on (p, m) 
(which follow from the bounds on /) and the bounds on (poo, moo)- To 
prove the lower bound, it suffices to establish the uniform convexity of 
F. Let 

I 1 2 

f{p,m) = plogp + 



2p 
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The Hessian of / has matrix 



P 



m 

7^2 



v 



m 



\m\ 



1 
P 



where Jtv stands for the N x N identity matrix; under our assumptions 
on p, this Hessian matrix is uniformly positive, so / is uniformly convex, 
and the same is true of the functional F : {p,m) i — > J f{p,m) dx + 
'^{0). This conclusion does not change when one imposes the linear 
constraints constraints p8.10|) . and the lower bound in p8.14|) follows. 

The last crucial step in the proof consists in the verification of 
Assumption ISl^iii). As in the previous subsection, 

2 



> K 

> K'WT 



VTI 



|{Vn}h 



18.15) (Id - Hi)'n,/ ■ (5Hi/) 



which controls ||Hi/ — H2/IP. 
Next, 



{\d-Ii,)'-{BIi,f) 



where the second inequality follows from a version of Korn's inequal- 
ity |141 eq. (1)]. Note that S = V(I]a:) (to avoid confusions I shall 
now write Ex for the map x Sa;), so one can apply again a Poincare 
inequality to obtain in the end 



;i8.i6) 



(Id - H2)'n,; ■ (5H2/) 



> K\\u — Sxl 



which controls IIH2/ — H3/IP. 

The gain from Tl^ is the main novelty. As a consequence of |15| 
eq. (65)], 

(18.17) 



M, 



dtp + T.- Vp 



P 

V — u 



N dtO 
Vp 



dtT. + (S ■ V) 
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The first and third lines do not bring any new estimate, so we focus on 
the second hne. First note that 



((s-v)s), 



jke 



(Do not mistake the symbol of summation with the matrix S.) Next, 
the equation for the mean velocity field u is dtu + u ■ Vn + VT + 
TV(logp) + (V ■ -D)/p, where D vanishes on the range of Hi. Taking 
the antisymmetric part of this equation results in dtV^u = 0, hence 
dtT, = 0. The conclusion is that dtT, vanishes on the range of U^- From 
all this information, we deduce that the second line of ()18.17|) can be 
simplified into 

where again S^x is a shorthand for the map x T?x. It follows that 

2 



{v — u) 



18.18 



(Id 



> K 



I shall now show that 



;i8.19) 



u — Sx 



Vp T?x 



>K{\\p-p^\\'' + \e-e^\^ 

+ - Soox|p). 

Since the right-hand side controls Ijlls/ — 114/ Ip"*"^, in view of p8.18j] 

and ()18.1(i|) this will imply 

(18.20) 

(Id-n2)'n,r(^n2/) '+ {id-Yi^)'-{Bi[^f) ' > i^||n3/-n4/f , 



completing the verification of Assumption IHl^iii). 

Since \6 — Oaol < C{\\p — poo|| + \\u — Moo||) and ||Sx — Soo^^ll < 
< C"||n — Uooll, to establish ()18.19|1 it is sufficient to prove 

'''' " ' " f)- 



18.21 



u — Sx 



Vp 

P 



9 



>K{\\p-Pc 



+ \\U — Ur 



In view of the bounds on p and u, and the uniform convexity of F 
(used above to check Assumption [TU^ii)), inequality (jl8.21|) will be a 
consequence of 



;i8.22) 





2 

+ 


Vp S^x 









> 



K[Fip, 



m] 



F(poo,mc 



The following lemma will be useful: 
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Lemma 64. Let ^ be a K -uniformly convex function, defined and 
differentiahle on a convex open subset of a Hilbert space Ti, and let 
K-.Ti^W^bea linear map. If X^o minimizes $ under the constraints 
A{X) = c, then 

$(X) - < {2Ky^ inf grad <I>(X) + A 

AG(Kcr A)^ 

Postponing the proof of Lemma |^ for the moment, let us apply it 
to the uniformly convex functional 



$(p, m) 
and the linear map 



A(p, m) 



\m\ 



Then 



grad$ = ^logp 



\m\ m\ 



2V epJ' 



and (Ker A)-*- is made of vectors A = (Aq, A), where Aq G M and A is a 

(constant!) antisymmetric matrix. So Lemma IMl implies 

(18.23) 



F{p,m)-F{p^,m^) < C inf ( 

AoeM; A*=-A \ 



\ogp- 



\m\ 



29p^ 



+ 



m 



— -Ax 
Op 



By Poincare inequality, the bounds on 6, and Korn inequality, 
^18.24) 



inf 

A*=-A 



m 



Ax 



< C inf 

A*=-A 



V(u - Ax) 



C 



inf 



On the other hand, 

12 2 



logp 



2^p2 



;i8.25) 



< 2 inf 

Ao 

< 2inf 

Ao 



logp 



logp 



~2^ 

I Ex I' 



Ao 



Ao 



\T,x\ 



\m\ 



C 



26 2V 
u — Sxp. 



By Poincare inequality, the second term in the right-hand side of p8.25p 
can be bounded by a constant multiple of ||Vm — As for the first 
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term, it can also be bounded by means of a Poincare inequality: 



inf 

Ao 



logp 



\Tjx\' 



Ar 



logp 



< C 

= c 



29 

V ( log p - 



logp 



|Sx| 



~2e 



Vp 

p 



All in all. 



inf 



logp - 



2^p2 



< C 



u — Sx 



Vp T?x 



This combined with fll8.23p and p8.24|) concludes the verification of p8.22|) . 
Then we can apply Theorem 0^1 and get the conclusion of Theorem EH 

□ 



Proof of Lemma [Hll Let $ : X ^{X^ + X). By assumption, 
is a minimizer of $ on Ker A. Since $ is JiT-convex and differentiable, 
the same is true of $, so that 

$(0) > $(X) - (grad'$(X), X) + ^\\X\\\ 

where grad' stands for the gradient in the space Ker A. It follows by 
Young's inequality that 



<^{X) - $(0) < 



grad'$(A:) 



2K 



But grad' $ is nothing but the orthogonal projection of grad$(X) 
(in Ti) onto Ker A; so 



grad'<|)(A:) 



inf 

Ag(Kcr A)J 



grad$(X) + A 



The conclusion of Lemma |M] follows easily. 



□ 



Remark 65. I don't know if the term in 112 can be dispended with 
in (jl8.2Up ; in any case this is an example where it is convenient to have 
the general formulation of Assumption IHI^iii), rather than the simplified 
inequality (jl4.9|l . In the next subsection, another example will be 
presented where this possibility is crucially used (see Remark IHZj). 



128 



18.5. Maxwellian accommodation. In this subsection fi^, will 
again be a bounded smooth connected open subset of with unit 
Lebesgue measure, but now the boundary condition will be the Maxwellian 
accommodation with a fixed temperature T^. Explicitly, 
(18.26) 

xG^fi,. ^ f^{x,v)= ( [ r{x,v')\v' -nldv'] 



where /"*" (resp. /~) stands for the restriction of / to {f ■ n > 0} (resp. 
{v-n < 0}), n is the mner unit normal vector, and is a fixed "wall" 
Maxwellian: 

M^{v) = - 



jv-i 



(The analysis would go through if one would impose a more general 
condition involving a reflection kernel C{v' — > i;), as in jlOl Chapter 1].) 
An important identity which follows from p8.26p is 

(18.27) Vx G dQ^, / fix, v) [v ■ n) = 0; 

equivalently, the mean velocity satisfies 

(18.28) Vx G (9fi^, u-n = 0. 

In this case there is only one conservation, namely the total mass. 
Without loss of of generality, I shall assume that the solution is nor- 
malized so that 

(18.29) J fdvdx = l. 

Then the unique equilibrium is the Maxwellian distribution with 
constant temperature equal to the wall temperature: 

_ Hi 

(18.30) foo{x,v) 



(27r)^/2' 

Theorem 66 (Convergence for the Boltzmann equation with Maxwellian 
accommodation). Let f be a solution of ()18.1|) in a smooth bounded 
connected spatial domain C M.^ with = 1- Assume that f 
satisfies the boundary condition p8.26|l . the conservation laws ()18.29j) 
and the uniform regularity estimates (|18.3p . Then 

\fs>0, ||/(t,-)-/ooL. =0(t-°°), 

where foo is defined by ()18.30p . 
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Proof of Theorem IHHl The proof follows again the same pat- 
tern as in all the previous theorems in this section. However, the Lya- 
punov functional is not Boltzmann's H functional, but a modified ver- 
sion of it: 

£{f)= [ f\ogf + -L [ f\v\^dvdx. 



2T 

Moreover, the sequence of nonlinear projection operators will be 
UJ = Mp^t; U2f = M,uT^; Usf = M,oT^; UJ = Miot^. 
In particular, 

log -\ 

2 2 

1 f ImP N f fT , T 
plogp+ — / p— + — / p — - log — 



T/ 2 2 \ T T 



From this it is easy to check Assumption [Zfii). 

The interesting features of this case reveal themselves when we try 
to check Assumption ISI First, by a classical computation, 

V{f) = -S'{f)-{Bf)= [ D{f{x,-))dx+ [ (log/ + l)(-t;.V./) 
+ /" iv\/J)\v\^dvdx 
= [ D{f{x,-))dx- [ (f log f{x,v) + f{x,v)^^)\v ■ n\dvdx 
Difix,-))dx+ [ /log {,,2 {v-n)dvdx. 



dn^xi 



e 



2T„ 



As before, he first term in the right-hand side is controlled below by 
K[S{f) — S(Ilif)Y~^'^ . The second term needs some rewriting. In 
view of (jl8.27|) and (|18.26p . we have, with the notation p — (x) = 
/ v) \v ■ n\ dv, 

I —flog — ^—^{v-n)dvdx = — f /log[— J {v ■ n) dv dx 

f log ( _'{,■- ) \v-n\dvdx 



vn<0 



p-M^ 



f\ 


\v ■ 


• n\ 




P-M^\ 


\v ■ 


n\ 



f\v-n\ log [ ^1 ) dvdx. 

v-n<0 
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This quantity takes the form of a nonnegative information functional, 
as a particular case of the Darrozes-Guiraud-Cercignani inequality |101 
Chapter 1]. The Csiszar-KuUback-Pinsker inequality will give an ex- 
plicit lower bound: For each x G dVL^., 



f\v ■ n\ loj 



D-n<0 



/ 


V 


n 






V ■ 


n\ 



dv 



> 



2p„(x) 
1 



2p_{x) 



f\v-n\ — p-Mu]\v ■ n\ 
f\v-n\ — p_Mw\v ■ n\ 



Ll{{ii-n<0};|i'-n| dv) 
2 

L^{\v-n\ dv) 



After interpolation and use of smoothness bounds, we conclude that 
(18.31) 

v{f) > K[s{f) - s{uj)Y^' +K\\f- P-M.||^;;^^^^,^|_| 1^,,,^,^^, 

where q is arbitrarily large and e is arbitrarily small. A useful conse- 
quence of pS.Hlj) is 



:i8.32) Vif)>K\\T-T^\\ltlg^^^ 



where again q is arbitrarily large. 

The other estimates are similar to the ones in the previous subsec- 
tions: 



^18.33) < 



(Id - n;)n,/ ■ (BUif) 



(Id - n^)n,/ ■ {BU,f) 



>i^(||VTf + ||{Vn}f); 



> KWV 



sym 



U\ 



{id-n',)u,fm3f) >K\\vpr 

Thanks to (fllO^ and ^EMu 



v{f) + (Id - n;)n,/ ■ mj) 



>K[\\^nlHn.) + \\T-Tv.\\l%n^^ 
> K'WT — T 



where the latter inequality comes from, say, the trace Sobolev inequal- 
ity if, say, p = {2N)/{N - 2) and q = 2{N - 1)/{N - 2). (If = 2 a 
slightly different argument based on a variant of the Moser-Trudinger 
inequality can be used to give the same result.) After interpolation one 
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concludes that 
(18.34) 

vif)+ (id-n;)n,/-(5ni/) 



> K\\T-n 



1 2+e 



> Ji'||ni/-n2/|| 



2+e' 



Next, if Qx is not axisymmetric, then the boundary condition p8.28|l . 
the Korn inequahty from |14j and the Poincare inequahty imply 



18.35) (Id - n'2)n,/ ■ {BU2f) > i^|| V^^^mH 



>i^'||nf >ir"||n2/-n3/|| 



2+e 



If Qx is axisymmetric, the previous argument breaks down, but we 
can use (jl8.32|) and replace (jl8.35j) by 

(18.36) 



V{f) + (Id - n'Jn./ ■ {BU,f) 



> K'\ 



u 



1 2+e 



> A'"||n2/-n3/|| 



2+e' 



where the but-to-last inequality follows from a trace Korn inequality 
f Proposition IA.21l in Appendix |A.22|) . 
Finally, 

(id-n'3)n3/-(fin3/)||'>/i||Vpf 

> A^||p-if > Ji'||n3/-n4/f+^'. 

Then Assumption ISU^iii) is satisfied, and one can use Theorem 1491 to 
prove Theorem inni D 

Remark 67. This is an example where the range of Hi is much 
larger than the set where the dissipation T) vanishes. Trying to devise 
a projection operator onto the space where V vanishes gives rise to a 
horrendous nonlocal variational problem whose solution is totally un- 
clear. On the other hand, inequalities (jl8.34j) and (jl8.36|) would be 
false without the contribution of 'C(/). In this example we see that 
the possibility to use the generalized condition appearing in Assump- 
tion |Hl^iii), rather than the simplified condition p4.9|) . leads to a great 
flexibility. 

18.6. Further comments. In many important situations (vari- 
able wall temperature, evaporation problems, etc.), one is led to study 
non-Maxwellian stationary solutions of the Boltzmann equation; then 
there is usually no variational principle for these solutions, and the 
mere existence of stationary solutions is a highly nontrivial problem, 
see e.g. [Hill HI- 



132 



From the technical point of view, the non-MaxweUian nature of 
the stationary state means that if one defines B = v ■ (transport 
operator) and C = Q (colhsion operator), then the equations -B/oo = 
and C/oo = cease to hold. No need to say. Theorem UHl collapses, and 
it is quite hard to figure out how to save it. 

There is a thin analogy with the (linear) problem of the oscillator 
chain considered in Subsection 19.21 in the case when the two tempera- 
tures are not equal; in that change of reference measure, based 
on Proposition El^ii), was at least able to reduce the problem to one 
of the type A* A + B, B* = —B. By analogy, one could imagine that 
a first step to come to grips with the quantitative analysis of stabil- 
ity for non-Maxwellian stationary solutions consists in re-defining the 
"antisymmetric" and the "diffusive" parts of the Boltzmann equation 
by performing some change of reference measure. Even this first step 
is nontrivial. 
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This last part is devoted to some technical appendices used through- 
out the memoir, some of them with their own interest. 

In Appendix IA.19I I have gathered some sufficient conditions for 
a probability measure to admit a Poincare inequality. After recalling 
some well-known criteria for Poincare inequality in M", I shall prove 
some useful results about tensor products; they might belong to folklore 
in certain mathematical circles, but I am not aware of any precise 
reference. 

Appendices IA.20I and IA.21I are devoted to some properties of the 
linear (kinetic) Fokker-Planck equation. First in Appendix lA. 201 1 shall 
prove a uniqueness theorem; the method is quite standard, although 
computations are a bit tricky. Appendix IA.21I is much more original 
and could be considered as a research paper on its own right: There 
I shall present a new strategy to get hypoelliptic regularization esti- 
mates. The method has the advantage to be very elementary, to avoid 
fractional derivatives as well as localization, and to yield optimal expo- 
nents of decay in short time. As Nash's theory of elliptic regularity, it 
is based on differential equations satisfied by certain functionals of the 
solutions. The results are nonstandard in several respects: They are 
global, directly yield pointwise in time estimates, and apply for initial 
data that do not lie in an L^-type space. I developed the method dur- 
ing a stay in Reading University, from January to March 2003; thanks 
are due to Mike CuUen for his hospitality. 

A closely related, but somewhat simpler strategy was found inde- 
pendently and almost simultaneously by Frederic Herau. I shall explain 
his method in Subsection IA.21.2l and develop it into an abstract the- 
orem of global regularization applying to the same kind of operators 
that have been considered in Part H] of this memoir. This extension 
grew out from discussions with Denis Serre. 

Finally, in Appendix IA.22I I gathered various technical lemmas and 
functional inequalities which are used throughout the memoir. I draw 
the attention of the reader to the "distorted Nash inequality" appearing 
in Lemma IA.191 which might have an interesting role to play in the 
future for "global" hypoelliptic regularization estimates. 
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A. 19. Some criteria for Poincare inequalities 

To begin with, I shall recall a popular and rather general criterion 
for Poincare inequalities in M". 

Theorem A.l. Let V e C^iW), such that e"^ is a probability 
density on M". // 

(A.19.1) rZZM^_Ay(x) .+00, 

2 |x'|— >oo 

then /i satisfies a Poincare inequality. 

Proof. The key estimate can be found in Deuschel and Stroock |16[ 
Proof of Theorem 6.2.21]: li w = \VV\'^/2 - AV, then for any h G 

(A.19.2) j wh^d^i<4: j \Vh\^ d^i. 

Let i?o > be large enough that > for |a;| > Rq. For 

R > Ro, define e{R) := [mf{w{\x\); \x\ > R}]'^; then e{R) ^ as 
R^ oo. So it follows from ()A.19.2;i that 

(A.19.3) / h^dfi<e{R) 4 / | V/ip rf/i - (inf w) / h"^ dfi 

J\x\>R \_ J J . 

Now let h e Ci(M",M) with J hdfi = 0. For any R > 0, let Br 
be the ball of radius R in M", and let fiR be the restriction of /i to Br 
(normalized to be a probability measure). Since Br is bounded, fiR 
satisfies a Poincare inequality with a constant P{R) depending on R, 
so 



j h^diiR<P{R) j \Vh\^diiR + (^j hdfx 



2 

R 



Of course /jr has density {h[Br]) ^^^^l|x|<ii- If R is large enough, 
then fi[BR] > 1/2, so 

(A.19.4) [ h^e-^<P{R)[ \Vh\^e-^ + 2([ hr~^^ 

J\x\<R J\3:\<R \J\x\<R. 

omce 

J he-^ = and e"^ is a probability density, 

\ 2 / „ X 2 ^ 

V\ _ u „-V \ / / u2-V 



(A.19.5) / he-'' = / he-'' < / K'e 

\J\x\<R J \J\x\>R / J\x\>R 

Plugging this into ()A.19.4|) . one deduces that 



(A.19.6) / h^e-^ <P{R) [ |V/ipe-^ + 3 / h 

J J J\x\>R 
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Combining this with (jA.19.3|l . we recover 

j h^e-^ < [P{R) + 12e{R)] j |V/ip e"^ - 3(inf j /i^e"^. 

So, if R is large enough that 3(inf w)e{R) > —1, one has 

- VI + 3(mfw)e(i?) y y ' ' 
This concludes the proof of Theorem lA.ll □ 

The sequel of this Appendix is devoted to Poincare inequalities 
in product spaces. It is well-known that "spectral gap inequalities 
tensorize", in the following sense: If each Li [i = 1, 2) is a nonnegative 
operator on a Hilbert space Tie, admitting a spectral gap Ke, then 
Li® I + I ® L2 admits a spectral gap n = min(Ki, K2). Now the goal 
is to extend this result in a form which allows multipliers. I shall start 
with an abstract theorem and then particularize it. 

Theorem A. 2. For i = 1,2, let he a nonnegative unbounded 
operator on a Hilbert space Hi, admitting a finite- dimensional kernel, 
and a spectral gap > 0. Let M be a nonnegative unbounded operator 
acting on 7^2; whose restriction to the kernel /C2 of L2 is bounded and 
coercive. Then the unbounded operator 

L = Li (g)M + I (g) L2 

admits a spectral gap n > 0. More precisely, for any nonnegative op- 
erator M < M, such that the restriction M\ic2 of M to IC2 satisfies 
XI < M < AI, one has 

. f K2 K2 \ Ki 

K > mm — , — , — A 

V 2 16A2' 2 

Theorem A. 3. (i) For i = 1,2, let (X^,/!^) be a probability space, 
and let Li be a nonnegative operator on Tii = L'^{fii), whose kernel is 
made of constant functions, admitting a spectral gap Hi > 0. Let fn be 
a nonnegative measurable function on X2, which does not vanish ^2- 
almost everywhere, and M be the multiplication operator by m. Then 
the unbounded operator 

L = Li^M + I (g) L2 

admits a spectral gap k > 0. More precisely, for any nonnegative func- 
tion m < m, lying in L'^{fi2), 

II 112 

K2 K2 W'^TT'W ri K,l n 

K > min — , — - — — pra fi 
' 2 ' 16 llmll^^' 2 " 
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(ii) More generally, for each £ G {1, . . . , A^} , let {Xi, /i^) be a prob- 
ability space, and let (1 < £ < N) be a nonnegative symmetric 
operator on Ti^ = L'^{fii), whose kernel is made of constant functions, 
admitting a spectral gap Ki; let fne be a nonnegative measurable func- 
tion on X^+i X ... X Xn, which does not vanish fi^^i ® . . . ® fiN-ali^ost 
everywhere, and let Mi be the associated multiplication operator. Then 
the linear operator 

N 

admits a spectral gap. 

Example A. 4. Let fi and u be two probability measures on M, each 
satisfying a Poincare inequality. Equip with the tensor measure 
fi ® u{dxdy) = fi{dx) fi^dy). Then 

L = -{d:d.+x^didy) 

is coercive on ® 

Proof of Theorem IA.2L Let Pi be the orthogonal projection on 
(KerL^)-*- in Tii. The spectral gap assumption means Li > KePe. Let 
M be the multiplication operator by m, then M > M. 

When applied to nonnegative operators, tensorization preserves the 
order: when A > A' > and B > B' > 0, one has A ® B > A (g) B' > 
A'®B'. Thus, 

Li (g) M + / ® L2 > KiPi ® M + K2I ® P2- 

So it is sufficient to prove the theorem when Li = Pi and fn = m E 

Let {e\)i>Q be an orthonormal basis for Tii, such that {e\)i<ki 
is an orthonormal basis of /Ci := KerLi; and let (e^)j>o be an or- 
thonormal basis for 7^2, such that {e^)j<k2 is an orthonormal basis of 
/C2 := Ker L2. Then {e\ ®e^)i,j>o is an orthonormal basis for Ti. More- 
over, the kernel of L is the vector space generated by {e\ ® e^)jjg^, 
where /C := j); i < ki,j < ^2}- So the goal is to prove 

K,{{P,®M)fJ)+K2{{I®P2)fJ)>K 4- 
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Cij el ® e|. 



First of all, 

(A.19.7) ((/ ® P2)/, f) = Yl ® ^26? = 5^ 

j>0;i>fc2+l 

Next, 

Yl CijCiy{Mele}) 

i>ki+l; j,j'>0 

i>A;i+l; jj'>k2+l i>fci+l; i>fc2+l; j'<k2 

+ Y CijCij\Me],e],). 

i>fci+l; j,j'<k2 

We shall estimate these three sums one after the other: 
- The first sum might be rewritten as 



i>fci+l j>k2+l 



CijCj j , 



E 

j>k2+l 



and is therefore nonnegative. 

- Similarly, the third sum might be rewritten as 

which can be bounded below by 

„2 



/ /' 



j<k2 



' E 



i>fci+l; j<k2 

- Finally, by applying the inequality ||Mej|| < A (j < /C2) and the 
Cauchy-Schwarz inequality twice, one can bound the second sum from 
below by 



-2^ E ^^^11^^? I 

j<k2 j>fci+l 



Y '^^^'^f 

j'>k2+l 



> 



-2AE, E 4. El E 

> -2A 



Cij'Cj/ 



E 



E 



i>fci+l;j<A:2 U i>fci+l; j>fc2+l 
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All in all, 

((Pi®M)/,/)>A Yl i 




i>ki+l;j<k2 i>ki+l; j>k2+l 



Combining this with ()A.19.7|) . we see that for all 6 G [0, 1] 

{Lf, f)>^2{{i® A)/, /) + ^ie{{Pi ® M)f, f) 



l2 



i>0;j>k2+l i>ki+l; j<k2 i>ki+l; j>k2+l 




K := mm 



To conclude the proof of Theorem IA.2[ it suffices to choose 

^ ^= ' 8^) • 



□ 



Proof of Theorem I A. 31 Let M be the multiplication operator 
by m. The restriction of M to constant functions is obviously coercive 
with constant A := j mdfi2, and M is bounded by ||m||j;^i/. Then (i) 
follows by a direct application of Theorem IA.21 After that, statement 
(ii) follows from (i) by a simple induction on A^. □ 

A.20. Well-posedness for the Fokker— Planck equation 

The goal of this Appendix is the following uniqueness theorem: 

Theorem A. 5. With the notation of Theorem for any /o G 
L2((1 + E)dvdx), the Fokker-Planck equation ()2.7|) admits at most 
one distributional solution f = f{t,x,v) G C(M+; ©'(M!^ x M^)) n 
L^^{R+;L\{l + E)dvdx))nLl^{R+;Hl{W^xR';;)), such thatf{0,-) = 

u 

Remark A. 6. The a priori estimates 

f'^dvdx = -2 j \S/vf\^ dvdx + 2n j f'^ dv dx 
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^ I fEdvdx = -2 I \VJ\^dvdx + n I f{l + E)dvdx 



dt 

allow to prove existence of a solution, too, for an initial datum /o G 
L^((l + E) dvdx); but this is not what we are after here. (Actually, 
an existence theorem can be established under much more general as- 
sumptions.) 

Before going on with the argument, I should explain why the unique- 
ness statement in Theorem [3 implies the one in Theorem IHl In that 
case. Proposition Efiii) can be applied even if W is only continuous: 
indeed, the differential operator VV"(a;)-V„ always makes distributional 
sense. So, if h is any solution of (j2.6|) . satisfying the assumptions of 
Theorem IHl then / := hpoo defines a solution of ()2.7|) . and it also 
satisfies the assumptions in Theorem [71 in view of the inequalities 

j f{l + E)dvdx< j fe'^dvdx = J h^e~^dvdx, 



\Vvf\ dv dx = j \Wv{pooh)\ dv dx 

<2 j \V^h\^ pl^dvdx + 2 j h'^\\/ ^p^\^ dv dx 

(1 + |t;ne-^(")e-^l (I |V,/i|'d/i+ / h'' dp 



< C sup 

x,v 

Proof of Theorem I A. 51 By linearity, it is enough to prove 

||/(T,-)||L^<e^ni/(0,-)IU^, 

which will also yield short-time stability. So let / solve the Fokker- 
Planck equation in distribution sense, and let T > be an arbitrary 
time. 

For any (f G C°°{t,x,v), compactly supported in (0,T) x x M", 
one has 

/ (^dti^ + V ■ Vx-v^ — W{x) ■ Vv(p + AyLf — V ■ Vv^p^ dt dv dx = 0. 

Since / G L°°([0, T]; L^(M"xM")), a standard approximation procedure 
shows that 

(A.20.1) j f{T, ■)ip{T, ■) dvdx- j /(O, ■)v'(0, ■) dvdx = 

f {dt^p + V ■ V — W{x) ■ Vv'^ + — V ■ VvP^ dt dv dx 
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for all ^ e ^^((0, T); ^^(M^ x M^J)) nC([0, T]; ^^(M^ x M^)), where / / 

stands for the integral over [0,T] x x M^. 

Let be C°° functions on M" with < x ^ 1; = 1 for 

< 1, x{^) = for |x| > 2, ?7 > 0, J?7 = 1, ?7 radially symmetric, 

ri{x) = for |x| > 1. With the notation e = {61,82), S = {61,62). 

Define 

Xe{x, v) = x{£ix) x{£2v), vs{x, v) = r]{6ix) ri{62v). 

In words: Xe is a family of smooth cut-off functions, and rjs is a family 
of moUifiers. (The introduction of 775 is the main modification with 
respect to the argument in |30t Proposition 5.5].) 
Define now 

fe,S ■= {fXe) * VS, ^e,S ■= Xe{{fXe) * VS * Vs)- 

The goal is of course to let 5 ^ 0, e — in a suitable way. 

Since rj is radially symmetric, the identity J g{f * rj) = j {g * ffjf 
holds true. So, for any t G [0,T], 

(A.20.2) j f{t,-)ip,4t,-)dvdx = j f,,s{t,-)^dvdx. 

Similarly, 

(A.20. 3) J fdt(ps,5dvdx = J fe,5 dtfs,5 dv dx = ^ J flsdvdx. 
By combining (IA.20.2D and (IA.20.311 . we get 

/(O, ■)V5£,5(0, ■) dv dx 



(A.20.4) l { I f{T,-)ip{T,-)dvdx - I f{0,-)ip{0,-)dvdx 



J J fdt(Pe,5dvdxdt = ^(^j f{T,-)ip^^s{T,-)dvdx- j , 

So, by plugging ip = ips,5 into ()A.20.1|) . one obtains 
1 
2 

(A.20.5) = J J fXe{v-V^,-VV{x)-V„ + ^,-vV,){fe,s*m)dvdxdt 
(A.20.6) 

+ I I f\{v-V,- VV{x) -V. + A.-v V,)xe] {fe,5 * Vs) dv dx dt 



(A.20.7) +2 j j fXe^vXe ■ yv{fe,5 * Vs) dv dx dt. 

For any given e > 0, all the functions involved are restricted to 
a compact set in the variable X = {x,v), uniformly in 6 < 1. 
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Now use the identities V{g * rf) = iVg) * rj, A{g * rf) = (Ag) * r], 
J gih * T]) = f h{g * rf) to rewrite ()A.20.5|1 as 

(A.20.8) j j fe,&{v ■ - VV{x) ■V^, + A.,-v- V^)fe,5dvdxdt 



where ^ is a temporary notation for the vector field (t;, —\/V{x) — v). 
By integration by parts, the first integral in ()A.20.8|) can be rewritten 

as 

(A.20.9) 



fe,5- 



Now we should estimate 
(A.20.10) 

^■V{fs,S*V5)-{^-'^fe,s)*V5 



L'2(dX) 



We shall estimate the contributions of v ■ Vx, v ■ Vy and VxV ■ 
separately. First, with obvious notation. 



\[v ■ Vx,r]s*] fsA 



L2 



{v - w) ■ VxfsAv^ ^) VSi - y) r]s^ {v - w) dw dy 
fe,&{y,w) (v-w) ■ VxVSiix - y) rjs^iv - w) dwdy 



L2 



L2 



Inside the integral, one has |f — w| < 62, \x — y\ < 61, and also |Vx%| 
0(6;^''^^^), 7]s, = 0(<52~"); so, all in all. 



[v ■ Vx,r]s*] fs,5\\r2 < C 



61 



n / feAy^^)'^\x-y\<SiM-"-M<S2dwdy 



L2 



^1 

<c^II/IIl2, 

Oi 

where the last two inequalities follow from Young's convolution inequal- 
ity. 

Next, 



IL2 



'v-w)-Vi,fe,5{y,w) r]s^{x-y) 7]s^{v-w) dwdy 
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Using the fact that It"— w| < 82 inside the integral and applying Young's 
convolution inequality as before, we find 

■ v„, m*]feA\L- ^ cd^WvjeAW < C62\\vj-\\l^. 

Finally, 

^[VV{x) - W{y)] ■VJeAy'^)VSi{x-y) r]s^{v - w) dw dy 
< C sup||Vr(x) - Vr(i/)|; \x-y\< 61, x,y G 

<ce,,{6i)\\vjeAL^, 

where stands for the modulus of continuity of ^ on the compact set 
Ki;. In all these estimates, the norm was taken with respect to all 
variables t, x, v. The conclusion is that the norm in ()A.20.10|) is 
bounded hke 



L2 



(A.20. 11) 



O 



^''""h+52\\VJ\\L2+9,,{5^)\\VJh^ 



Si 



Next, since ||Vt,Xe||Loo < Ce2, it is possible to bound ()A.20.7|) by 
(A.20. 12) 

Ce2\\fXe\\L4^v{fe,S * r]s)\\L^^ < || / || ( || V./ lU^ + Wfh^). 



Finally, the terms in the integrand of ()A.20.6|) can be bounded with 
the help of the inequalities 

\v-V,Xeix,v)\ < C\v\e,, \WV{x)-V,Xe{x,v)\ < Ce2M{e^'), 
\KXeix,v)\ < Cej, \v ■ V^Xeix,v)\ < C\v\e2, 

where M{R) := sup{|VV^(x)|; |x| < 2R}. Then, by Cauchy-Schwarz 
again, ()A.20.6|) can be bounded by 



C 



61+62(1+ M{e]:^)) Jj j p\v\^ dv dx dt J j j {f,^s*Vsydvdxdt. 



Since \v\'^ <2E - 2(inf y), in the end (jA.20.fijl is controlled by 
(A.20. 13) 



C 



ei + 62(1 + M{sA)) 



f^E dv dx dt\ 



P dv dx dt. 
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By plugging the bounds ()A.2n.9jl . ()A.2n.13;) and ()A.2n.12jl into ()A.2n.4jl . 
we conclude that 

(A.20.14) fls{T,x,v)dvdx- j fls{0,x,v) dvdx 

<— y y I Vt,/e,5|^(t, X, f ) (it; (it + — y y f^g{t,x,v) dvdxdt 



^(^ll/lli^ + ^2||/||Ln|V./.,5||L^ + ^.,(5l)l|V./e,5||LHI/MllL^ 

+ C(ei + e2M(er'))||/(l + i?)|| 



2 

L2, 



< liminf 



where all the norms in the right-hand side are with respect to 
dvdxdt. Now let 62 — > 0, then 61 0, then 62 — > 0, then ei — * 0, 
then 5 — > 0: all the error terms in the right-hand side of ()A.20.14jl 
vanish in the limit, and f^^s converges to / almost everywhere and in 
(dvdxdt). So 

y /^(T, X, t') (if (ix < liminf J ^{T, x,v) dv dx 

y fs,5{^-:^-:'^)dvdx+ — j j f^g{t,x,v) dvdxdt 

/^(O, x,v) dvdx +^ J J f^{t, X, v) dv dx dt. 
By Gronwall's lemma, 

nt 

Wfii, ■)IU2(]RnxMn) < e~||/(0, ■)I|l2(IR"xIR"), 
which concludes the argument. □ 

Remark A. 7. Just as in |30| Proposition 5.5], the particular struc- 
ture of the Fokker-Planck equation was used in the estimate |Vy(x) ■ 
V^XeI ^ £2M{e^^). It would be interesting to understand to what 
extent this computation can be generalized to larger classes of linear 
equations, and whether this has anything to do with the hypoelliptic 
structure. 



A. 21. Some methods for global hypoellipticity 

This Appendix is devoted to various regularization estimates for 
the Fokker-Planck equation. I shall consider only two particular cases 
(those which were used in the present paper): First, the L'^ifJ^) — > 
H^[n) regularization for the Fokker-Planck equation in the form ()7.1|) : 
secondly, the M —* regularization for the Fokker-Planck equation 



A.21. SOME METHODS FOR GLOBAL HYPOELLIPTICITY 
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in the form ()7.9j) (Here M is the space of bounded measures, and H'^ 
is the non-weighted Sobolev space of order k.) 

A.21.1. From weighted to weighted H^. In the sequel V , 
is a potential on M", bounded below, 7(1;) = (27r)~"/^e~l''l is 
the standard Gaussian, and ix{dxdv) = 7(t>)e~^^^^ dv dx stands for the 
equilibrium measure associated with the Fokker-Planck equation ()7.H) 
(it might have finite or infinite mass). Apart from that, the only regu- 
larity assumption is the existence of a constant C such that 

(A.21.1) |VV| < C(l + |VK|). 

As we shall see, this is sufficient to get estimate (j7.8j) . independently 
of the fact that satisfies the Poincare inequality ()7.5|) or not. 

Theorem A. 8. Let V he a C"^ function on M", hounded helow 
and satisfying ()A.21.1|) . Then, solutions of the Fokker-Planck equa- 
tion ()7.1|) with initial datum Hq satisfy 

3 ^ 

k=i 

for some constant C , only depending on n and the constant C appearing 
m (IA.21.1j) . 

Remark A. 9. These estimates seem to be new. The proof can be 
adapted to cover the case of initial data, at the price of a deterio- 
ration of the exponents. I shall explain this later on. 

Remark A. 10. Theorem IA.8I shows that (with obvious notation) 
e~*^ maps into H^nH^ with norm 0(t"'^/^). It also maps into 
with norm 0(1); so, by interpolation, it maps into if" fl H^°' with 
norm 0{t~^°'^'^), for all a G [0, 1]. Since e~*^^"^^^ dt is a parametrix 
for (/ + L)~^, and is integrable at t = for /5 < 1, one can deduce 
a "stationary" hypoelliptic regularity estimate a la Kohn: 
(A.21.2) 

\\h\\Hs{i.) + WHmHi^) < C{\\h\\m^) + WLHWl^^^)), Va < 2/3. 

With a much more refined analysis, it is actually possible to catch the 
optimal exponent a = 2/3 in the above estimate. (This realization 
came after discussions with Christ.) I shall not develop this tricky 
issue here. 

Proof of Theorem IA.8L As a consequence of TheoremlHland an 
approximation argument which is omitted here, it is sufficient to prove 
this theorem for smooth, rapidly decaying solutions. So I shall not 
worry about technical justification of the manipulations below. Also, 
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C will stand for various constants which only depend on n and the 
constant in ()A.21.1|1 . 

The following estimates will be used several times. As a conse- 
quence of Lemma [a. 181 in Appendix lA. 221 for each 

by integrating this with respect to 7(1^) dv one obtains 
(A.21.3) f \W\^g^dii<C ( [ g^dfx+ [ |V^^|^rf/i 
Similarly, 

(A.21.4) / \v\'^g'^dfi<C ( f g^ dfi + ! |V^^|^rf/i). 

Now we turn to the main part of the argument, which can be de- 
composed into four steps. 

Step 1: '^Energy" estimate in H], and norms combined. 

To avoid heavy notation, I shall use symbolic matrix notation which 
should be rather self-explanatory, and write 

L = vVx- W{x) ■ V„ - A„ - ■ V„. 

By differentiating the equation once with respect to x, and three 
times with respect to one finds 

(A.21.5) (^^ + V./i = WlV{x) ■ V./i, 

(A.21.6) (^^ + Vlh = -SVlV^h - 3Vlh. 

After taking the scalar product of ()A.21.5j) by Vxh and integrating 
against /i, we get 

(A.21.7) j |V../ipc//i + J |V„V./i|2d/i = j{\/lV)\/yh-Vxhd^i. 

Similarly, from ()A.21.6|) it follows that 
(A.21.8) 

~ J \\/lh\'' di^+ J \Vth\^dfi=-3 J \\/lh\'' dfi-3 J Wlh-WlW^hdn. 
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Let US bound the right-hand side of ()A.21.7|1 . Since (V„)* = + 
V (where the * is for the adjoint in L^(yu)), one has 

J {VlV)V,,h-S/^hd^i = - j {Vly)h-V^V^hd^i- j {{VlV)hv ,V ^h) d^i. 

By Cauchy-Schwarz and Young's inequahty, 

- j {^lV)h-V^V^hd^< j \VlV\^h^dfx+^ j |V^V^/i|^rf/i. 

Thanks to (jA.21.3|) . this can be bounded by 

c(^j |V^/i|2rf/i + j h^df?j+^J \\/^V^h\^dix. 

By Cauchy-Schwarz inequahty again, 
(A.21.9) 



- j{{^ly)hv, V.h)d^<^ j \S/lV\^h^dfi^ j \v\^\V,h\^dfi. 
In view of flA.21.4;) . 

(A.21. 10) \v\'\V.h\^dfx<C Q \\/.h\^dfx + j |V„V./i|'rf/i^ . 

By flA.21.9|) . ()A.21.10|) and Young's inequahty, there is a constant C 
such that 



+ 4 



AU in aU, 
(A.21. 11) 



The right-hand side in ()A.21.8jl is estimated in a similar way: 



then on one hand 
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on the other hand, again by ()A.21.4|) . 



-3 j Vlh-vVyV^hd^i<^J j |t;|2|V3/i|2d/i J j |V„V^/i|2d/i 



4 J ' • - ' - ■ 4 

So there is a constant C such that 

(A.21.12) 
ld_ f 



^ ^'V^/i|'d/i+J / \Vth\^dfi<Ci I \VyV,h\^dfi+ I \Vlh\^ dfi 



As a consequence of ()A.21.11|) and ()A.21.12|) it is possible to find nu- 
merical constants a,K,C > (only depending on n and C in ()A.21.1|) ) 
such that 

(A.21.13) 

+ c(^Jh^dfi + J \V,h\^dfi + J \Vlh\^dii 

This concludes the first step. 

Step 2: Time-behavior of the mixed derivative. 

In this step I shall focus on the mixed derivative integral / Vxh ■ 
V„/i dfi. By differentiating the equation with respect to x and multiply 
by V„/i, differentiating the equation with respect to v and multiply 
by Vxh, then using the chain rule and the identity FA^G + GA^F = 
Ay{FG) — 2V„F ■ V„G, one easily obtains 



d_ 

dt 



{Vxh ■ Vyh) = {VlV ■ V„/i, Vyh) 

- 2VvVxh ■ Vlh dfi - I Vx/i|^ - Vxh ■ Vyh. 
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After integration against /i, this yields 
(A.21. 14) 

V^./i-V^/ic?/i = J {WlV ■V^,h,V^h)dn-2 j W^W ^h-Wlkd^i 

-J \V^h\^dfx- J W^h-Wvhdfx. 
The first term in the right-hand side need some rewriting: Since (Vt,)* = 

j {VlV ■Vyh,V^h)dfi = - j VlVhVlhdfi- j h{VlVv,Vyh) dfx 



With the help of Young's inequality and ()A.21.4|) again, this can be 
bounded by 

e j \\/lV\^h^dfx + Ce(^+ J \Vlh\^dfi + J \V,h\^df?j. 

By Lemma [a. 181 if e is small enough then this is bounded by 

^(^j \V^h\^dfi + j h^di^+c(^j |V^/i|^d/i+ j \Vlh\^dii 

Now for the second term in the right-hand side of ()A.21.14|) . we 
just write 

-2 j V^V.^hVlhd^i< j \V^VM^d^i + j \Vlh\^ d^i. 

Summarizing all the above computations: There is a numerical con- 
stant C, only depending on n and C in ()A.21.1|) . such that 

(A.21. 15) V,h-\/,hdfi<~J \V,h\^dfi 

+ c(^jh^dfi + j\Vyh\^dfi + j\Vlh\'^dfi+ j \V dii^ . 

This concludes the second step of the proof. 

Remark A.ll. We could also have conducted the computations in 
the following way: 

-2 j Vv'^c,hVlhdn = 2 J V -Vlhd^i - 2 j V^h-vVlkdn. 
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Then on one hand, 

2 j V^h-Vlhdfi |V^/i|^rf/i + 4 j |V^/i|^d/i; 

on the other hand, just as before, 

-2 J Vcch-vVlhdfi<^ J |V^/iprf/i + 4 j \v\'^\Vlh\^ d^ 

<\j l^^.hl^dfi + C (^j \Vlh\^dfi + j \Vlh\^dfiy 

By doing so, we would have obtained the same result as ()A.21.15|) . ex- 
cept that the integral / | Va;Vt,/ip dfi would be replaced by / | V^/ip (i/i. 
Then the rest of the proof would have worked through. 

Step 3: Interpolation inequalities 

If /i is a function of v, lying in L^i'y), one can write h = CkHk, 
where Hj. are normalized Hermite polynomials and k are multi-indices 
in N"; then 

jh^d^ = Y,cl j \vM'd^^ = Y.\k\\l 

j \ylh\''d^^ = Y,\k\'cl etc. 

(here \k\^ = I < i < n and \k\^ = {\k\y). Then, by Hdlder's 

inequality (in the k variable), one can prove interpolation inequalities 
such as 



Now a h = h{x,v) is a function of both variables x and v, one can 
apply the previous inequality to h{x, ■) for each x, then integrate with 
respect to e~^{x) dx, and apply Holder's inequality in the x variable, 
to find 

J \Vyh\^di2<C (^J h^di?j ' [J \Vlh\^di^ ' . 

Similarly, 

\Vih\^ djj < C i dfij U \Vth\^dfij , 1<J<3. 



Step 4: Conclusion 
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Now we can turn to the proof of estimate ()7.8|1 . Without loss of 
generahty, assume f h"^ dfi = 1 a.t time 0. Then, since this quantity is 
nonincreasing with time, J h'^{t, ■) dfi < 1 for all t > 0. By combining 
the results of Steps 1, 2 and 3, we discover that the quantities 



X 



\Vxh\^dfi, 



M :-- 



\Vih\^dfi 



(0<J<4), 



W 



|V,V,/ip 



viewed as functions of t, solve the system of differential inequalities 
^ j^{X + aFa) < -K{Y^ + w) + C{l+X + Y^)- 



(A.21. 16) < 



j^M < -KX + C(l + Y1 + Y2 + W) 

\M\ < VXY\; Yi < CY]!'^ < C'Y^^'^ < C'yI'^ 



It is a consequence of Lemma FA. 201 in Appendix lA. 221 that solutions 
of (IA.21.16D satisfy 

< t < 1 ^ Xit) + Yit) < 4 

for some computable constant A. As a consequence, for < t < 1, 

j \V,h\^dfx = Oit-''), j \Vlh\^dfx = 0{t-^). 

Then, by interpolation / |V>|2 = O(t-i), / jV^/ip = 0{t~^). This 
concludes the proof of (|7.8|) . □ 

A. 21. 2. Variants. Here I studied the regularization effect by means 
of a system of differential inequalities. It is natural to ask whether one 
can do the same with just one differential inequality. The answer is 
affirmative: It is possible to use a trick similar to the one in the proof 
of Theorem ^1 that is, add a carefully chosen lower-order term which 
is derived from the mixed derivative j Vxh ■ V„/i. 

A first possibility is to consider the Lyapunov functional 

8{h) = J h^dfx + a j |V^./i|2rf/i + 26 J V,,{Dy^h)-V^{Dy^h)d^i 

+ c |V^/ipd/i, 
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where = (— A^)^/^. Then by using computations similar to the ones 
in Subsection IA.21.l| plus estimates on the commutator [Dy'^,VV^], 
one can establish the following a priori estimate along the Fokker- 
Planck equation: For well-chosen positive constants a, b, c, 

^ S{h) < -K8{hYl\ h = e-*^/io. 

The desired result follows immediately. 

One drawback of this method is the introduction of fractional deriva- 
tives. There is a nice variant due to Herau |31j in which one avoids 
this by using powers of t: 

J^(t,h) = j h^dfi+at j \V^h\'^d^i+2ht^ j V^h-V^h+ct^ j |V^/i|^ci/x. 

Then one can estimate the time-derivative of e~*^ho) by means of 
computations similar to those in Subsection lA . 2 1 . Tl and the inequalities 



t 



<C / \V^h\^dfi + et^ / |V^/i|^d/i; 



J W^Vyh-Vlhdfx <Ct J \\/lh\'^dfx + et^ j iV-^V^M'^df^. 



In the end, if a, b, c are well-chosen, one obtains, with the shorthand 



-tL 



ho, 



j^7{t,h)<-K (^j |V„/ipc//i + t j \\/lh\^ dfx + J \\/xh\^dfi 

j \VxV^,h\^df?j . 
It follows that jF(t, h) is nonincreasing, and therefore 

|V„/i|2d/i = O(r^), / |V^/i|2d/i = 0(t"^). 



The conclusion is not so strong as the one we had before, since we only 
have estimates on the first-order derivative in v. But the exponents are 
again optimal, and it is possible to adapt the method and recover esti- 
mates on higher-order derivatives. Furthermore, estimates on j iV^^/ip 
and J |Vt,/?.p are exactly what is needed for Theorem IHTI to apply. 

Herau's method lends itself very well to an abstract treatment. For 
instance, let us consider an abstract operator L = A* A + B, satisfying 
Assumptions (i)-(iii) in Theorem ITHl then the following decay rates (in 
general optimal) can be proven, at least formally: 

(A.21.17) ||Ae-*^/i|| =0(t-^/2); ||Ce-*^/i|| = 0(r=^/'). 
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To show this, introduce 

h) := \\e-'^hf+at\\Ae-'^hf+2bt^{Ae-'^h, Ce-'^h)+ct^\Ce-'^hf. 

Then, we can perform computations similar to the ones in Subsec- 
tion except that now there are extra terms coming from the time- 
dependence of the coefficients a, b, c. Writing h for e~*^/i, we have, if 
a, b/a, c/b, c^/b, b'^/ac are small enough: 

(A.21. 18) 

< -4 Uhf + at\\A^hf + bt^WChf + ct^WCAhf) 
at \ / 

+ a\\Ahf + Abt{Ah, Ch) + 3ct2||C/if , 

were k is a positive number. When < t < 1, the positive terms in 
the right-hand side of ()A.21.18|) can all be controlled by the negative 
terms if a, b and c/b are small enough. Then 

< -K{\\Ahf + p2/if + \\Chf). 

(JjV 

In particular, JF is a nonincreasing function of t, and then the desired 
bounds \\Ahf = 0(ri), \\Chf = 0{t-^) follow (as well as the bound 
\\A^hf = 0{t-^)). 

The very same scheme of proof allows to establish a regularization 
theorem similar to Theorem! 



Theorem A. 12. Let H be a Hilbert space, A : H ^ H"^ and 
B : H ^ H be unbounded operators, B* = —B, let L := A*A + B. 
Assume the existence of ^ 'N and (possibly unbounded) operators 
Cq, Ci, . . . , CtVc+i; -Ri, • • • , Rnc+i '^^^ . . . , Z]\;^_^i such that 

Co = A, [Cj,B] = Z,+iC,+i+Rj+i (0<j<iVe), Cjv.+i = 0, 

and, for all k & {0, ... , N^}, 

(i) [A,Ck\ is bounded relatively to {C*j}o<j<fe and {CjA}o<j</c-i; 

(a) [Ck,A*] is bounded relatively to I and {C'j}o<i<fc; 

(Hi) Rk is bounded relatively to {C*j}o<j</t-i CLnd {Cj74}o<j<fc-i- 

(iv) There are positive constants Xj, Aj such that Xjl < Zj < Ajl. 

Then the following bound holds true along the semigroup e~^^ : 

\h\\ 



VA;G{0,...,Arj, ||e-*^/i||<C 



where C is a constant only depending on the constants appearing im- 
plicitly in Assumptions (i)-(iv). 
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Remark A. 13. A reasoning similar to Remark lA. 101 shows that the 
exponents l/{k+ 1/2) cannot be improved. Indeed, in Hormander's 
theory, the weight attributed to the commutator Ck would be 2/c + 1, 
and the regularity estimates established by Rothschild and Stein |44j . 
which are optimal in general, provide regularization by an order 2/(2A;+ 
l) = l/(fc + l/2). 

Remark A. 14. I shall show below how Herau's method can be 
adapted to yield regularization from LlogL initial datum. On the 
other hand, it is not clear that it can be used to establish regularization 
from measure initial data. 

A. 21. 3. Higher regularity from measure initial data. Now 

I shall explain how to extend the previous results by (a) establishing 
Sobolev regularity of higher order, (b) removing the assumption of 
integrability for the initial datum. 

I shall only consider the case when VV^ is Lipschitz and has all its 
derivatives uniformly bounded. There are three motivations for these 
restrictions: (i) even if they are far from optimal, they will simplify the 
presentation quite a bit; (ii) they ensure the uniqueness of the solution 
of the Fokker-Planck equation starting from a measure initial datum; 
(iii) the theorems of convergence to equilibrium studied in the present 
paper use the Lipschitz regularity of W anyway. 

As before, the equation under study is 

(A.21.19) dtf + v- VJ - VV{x) ■ VJ = AJ + v ■ VJ + nf. 

This equation admits a unique solution as soon as /o is a finite non- 
negative measure (say a probability measure) with finite energy, and it 
is easy to prove the propagation of regularity and of moment bounds. 

So to establish regularization in higher-order Sobolev space H^{Hl) 
it is enough to prove, for smooth and rapidly decaying solutions, an a 
priori estimate like 

\\ft\\H>iHi(R^xR^) < piy;:, 

with constants C and n that do not depend on the regularity of /q. 

In the sequel, C and K will stand for various constants depending 
only on n and V. As in Subsection IA.21.11 the a priori estimate is 
divided into four steps. The conservation of mass (that is, the preser- 
vation of J f dx dv) along equation ()A.21.19|) will be used several times. 



Step 1: "Energy" estimate in higher order Sobolev spaces. 
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Let k and £ be given integers {k will be the regularity in x and i the 
regularity in v). Computations similar to those in Subsection IA.21.'l] 
(differentiating the equation and integrating) yield 



I / iVyjfdxdv <-K I \V',Vi+'f\'dxdv+C I \V^Vif\'dx 
+ C J \V^,+'Vi-^f\^dxdv + C J2 J IV^^VJHV^+Vpdxci^;. 



dv 



l<i<k ■ 



By assumption |V^y| is bounded for any i, so the above equation 
reduces to 

I / iv^v^/p <-kI |vM^VP + c^E / ivM/r 



c I iv^vr/p. 



Then one can repeat the computation with (A;, i) replaced by {k + 
l,i — 2) and then (/c + 2, £ — 4), etc. By an easy induction, for a given 
integer m, we can find positive constants K, C, Qq = 1, ai, . . . , am such 
that 



d 



III n lit n 

^afc / \V\vT-''^j\^dxdv<-KY, \ iVtvf^^-'^+'fl'dxdv 

k=0 k=0 

■m „ 

\^i^i^"~''^f\^dxdv. 

7 r, 1^ ^ 



fc=0 3<k 



Repeating the same operation for lower order terms (that is, de- 
creasing m), for each couple of nonnegative integers (/c, with 3fc + £ < 
3m we can find a positive constant afc,£ such that 

^ E ak,iJ\W''yj\'dxdv<-Kj\Vl"'+'f\''dxdv 

+ C E / l^^'^fl'dxdv. 



3k+e<3m 



3k+e<3m ' 



Then we can define an "energy functional" of order m, which controls 
the L^-regularity of / up to order m in a; and 3m in v. 

(A.21.20) £m{f)= Yl a.^i j \^Kf?dxdv. 



3fc+£<3m 
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(Recall, to avoid any confusion, that / | Vj;V:^/p is the sum of all terms 

lid'^l ■ --QllKx ■ --dtf? with k, + ... + K = k,i, + .. . + L = t) 

Then the a priori estimate on the Fokker-Planck equation ()A.21.19j) 

can be recast as 

(A.21.21) j^SUf) < -K j \Vl"'+'f\''dxdv + C£Uf)- 

The important terms in Sm are the extreme ones, that is for {k, i) = 
(m, 0), (0, 3m) or (0, 0). All the other ones can be controlled by these 
three extremal terms; to see this, it suffices to apply Holder's inequality 
in Fourier space: Denoting by ^ the conjugate variable to x and by rj 
the conjugate variable to v, one has 

I \VlV'j\'dxdv = C j \e'\7if\J\d£,dr^ 

k i 

<cQ l^nfl'd^dr^y Ivnll'd^dr^Y"^ {^j \f\'d^dr^J 

k I ■ ' 

= c(^j iV^^fl^dxdv^ (^j iVl"^ fl'^ dx dt^'^ (^j fdxdv^ 
It follows easily that there are positive constants K, C such that 

(A.21.22) k(^I \VTff + I |Vr/r + I <Sm{f) 
Step 2: Mixed derivatives 

Now define the higher order mixed derivative functional 
(A.21.23) 

l<il,...,i,n<n ^ '™- ^ ^»m-l *'»™ 

By computations in the same style as in Step 2 of Subsection lA.21.T| 
one can establish 
d 

dt' 

k<m, 3k-\-£<3m " 

Each of the terms appearing in the latter sum can then be estimated 
by elementary interpolation inequalities as in Step 1: If < m then 



;MUf) < -K j \VZ'f\''dxdv + C j |VM/r- 
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where e is an arbitrarily small positive number. The conclusion is that 
(A.21.24) 

j^Mm{f)<-K j \V^f\^dxdv +C (^j \Vl"'f\^dxdv + j fdxdv 
Step 3: Interpolation inequalities. 

There are two things to check: (i) that is "much smaller" than 
Sm, and (ii) that / iVj^^/p is "much smaller" than / iV^J^+Vp. The 
difficulty is that we cannot just use interpolation in L^-type spaces. 
In replacement, we shall use the anisotropic Nash-tjpe interpolation 
inequality exposed in Appendix IA.221 

First, by Cauchy-Schwarz, 

\Mmif)\ < (/ iv:^/py (1 ivrv./py . 

Then the second term is estimated thanks to Lemma IA.19I with A = 
m — 1, fi = 1, y = m, fi' = 3m: 

ivrv./p <( f I v™/p + [ ivir^ff) ( ! i 



where 9 = 2/ (3m + Qn) is a positive number. Since the mass J f is 
preserved under the time-evolution by the Fokker-Planck equation, we 
arrive at the estimate 

(A.21.25) \Ms{f)\<C£m{fy-^ 

where 5 = 6'/2 is a positive constant. 

Next, apply Lemma IA.19I again with A = 0, A' = m, /i = 3m, 
fi' = 3m + 1. Noting that (A/A') + (/i/ /i') = 3m/ (3m + 1) < 1, we see 
that there exists ^ G (0, 1) such that 

j iVl^'fl^dxdv <C (^j |V™/P + j iVl'^+^flUxdv^ (^j fdxdv 

The same estimate holds true for j (this can be treated by the usual 
Nash inequality), and then one can use the fact that / /is preserved 
by the Fokker-Planck equation, to obtain the a priori estimate 
(A.21.26) 

J iVl^'fl^dxdv+j fdxdv<c(^J |V^/P + j iVl^'+^fl^dxdv^ 



Step 4: Conclusion 
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Equations (!A.21.22jl . ()A.21.25jl . ()A.21 21jl ()A.21.26jl and (!A.21.24;i 
together show that we can apply Lemma rA. 201 with £ = Sm, M. = Aim-, 
X = f |V^/|2, Y = f |V^'"/|2 + //2, z = / |V3'"+VP. Thus there 
are constants C and k such that Smift) < C/t^^'^. This concludes the 
proof of the a priori estimate. 

A. 21. 4. Regularization in an LlogL context. If the initial da- 
tum is assumed to have finite entropy, then Herau's method can be 
adapted to yield the regularization in Fisher information sense, with 
exponents that are likely to be optimal. Here is a rather general result 
in this direction, under the same assumptions as Theorem |2H1 

Theorem A. 15. Let E e C^(]R^), such that is rapidly de- 
creasing, and jji{dX) = e~^^^^ dX is a probability measure on M^. Let 
{Aj)i<j<m dnd B be first-order derivation operators with smooth co- 
efficients. Denote by A* and B* their respective adjoints in L'^{fi), 
and assume that B* = —B. Denote by A the collection [Ai, . . . , Am), 
viewed as an unbounded operators whose range is made of functions 
valued in M™ . Define 

m 

L = A* A + B = Y^ A*Aj + B, 

and assume that e~^^ defines a well-behaved semigroup on a suitable 
space of positive functions (for instance, e~*^h and log(e~*^/i) are C°° 
and all their derivatives grow at most polynomially ifh is itself C°° with 
all derivatives bounded, and h is bounded below by a positive constant). 

Next assume the existence of Nc > 1, derivation operators Cq, . . . , Cat^+i 
and Ri, . . . , Rnc+Ij (^''^d vector-valued functions Zi, . . . , Zjsi^j^i (all of 
them with C°° coefficients, growing at most polynomially, as their par- 
tial derivatives) such that 

Co = A, [C„B] = Z,+^C,+^+R,+^ (0<j<A^J, C^.+i = 0, 
and 

(i) [A, Cfc] is pointwise bounded relatively to A; 

(a) [Ck,A*] is pointwise bounded relatively to J, {C*j}o<j<fc; 

(Hi) Rk is pointwise bounded with respect to {Cj}o<j<fc-i; 

(iv) there are positive constants such that Xj < Zj < Aj; 

(v) [A,Ck]* is pointwise bounded relatively to I, A. 
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Then the following bound holds true: With the notation h(t) = 

0; 



G {0,...,A^J, J h{t)\Ck\ogh{t)\'^ dfx < C 



ho log ho dfi 



^2k+l ' 

where C is a constant only depending on the constants appearing im- 
plicitly in Assumptions (i)-(v). 



Proof. The proof is patterned after the proofs of Theorems 
and IA.121 Write u = log h, f = e~^h, and introduce the Lyapunov 
functional 

J^{t,h) = I fu + Y,(akt''+' I f\Cku\l+2bkt''+' I f{CkU,Ck+iu)m). 

k=0 

The computations for dT/dt are the same as in the proof of The- 
orem I2H1 except that now there are additional terms caused by the 
explicit dependence on t. So 

,A.21.27) 




k k 

+ Y,{'^k+l)aut'^ I /|Cfcn|2 + ^(2A; + 2)M''+' / /(C^n, Cfc+m). 
k k 

Obviously, the additional terms can be controlled by the ones in the first 

line of the right-hand side, provided that au/bk-i and bkt^^'^'' / \/ibk^Iit^^)(bk^^^^ 

are small enough; the second condition reduces to bk/bk-i small enough. 

These conditions have been enforced in the proof of Theorem |2H1 So 

all in all, ^(t, h{t)) is a nonincreasing function of t, and the conclusion 

follows immediately. □ 

A.22. Toolbox 

The following elementary lemma is used in the proofs of Theo- 
rems [H iOl and 123 

Lemma A. 16. Let 5 > and uq > be given. Then it is always 
possible to choose positive numbers Ui, U2, ■ ■ ■ , Un in such a way that 

\k e {0,...,N -1}, Uk+i<6uk; 
VA; G {1, . . . ,iV - 1}, ul < 6uk^iUk+i. 
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Proof. Without loss of generality, assume uq = 1. Set mo = 
0, nil = 1; by induction, it is possible to pick up positive numbers nik 
such that 

mfc+i e {rrik, 2mk - ruk^i). 

The resulting sequence will be increasing and satisfy > (m^.i + 
mfc_|_i)/2. Next set = e"^''-; for e small enough, the desired inequalities 
are satisfied. □ 

The next lemma is used in the proof of Theorem EHl it is a kind of 
nonlinear analogue of Lemma IA.16I 

Lemma A. 17. Let K,E,k > 0, J 6e given. Then there exists con- 
stants El = ei{J) > 0, £ = £(J,fc) > and Ki = Ki{K,E,k,J) > 
with the following property: For any e G {0,ei) and E G {0,E), there 
exist coefficients ai, . . . , aj_i > satisfying 

1 = flo > oi > ^2 > . . . > aj_i; 



(A.22.1) 



ai < KE^] 
VjG{l,J-l}, 



«7-l 



< Kafj^E^^- 



^aj.i>KiE'^. 

Proof of Lemma IA.17L Without loss of generality we may as- 
sume that K is bounded above by m := min(l, (E")"*^); otherwise, just 
replace K hj m. 

We shall choose the coefficients aj in such a way that the inequality 
in the third line of ()A.22.H1 holds as an equality. For j = J — 1 this 
gives 

hence 

aj.2 = aj.i {KE^'a'j_i)-\ 

Then the equality 

2 



a,+i 



yields, by decreasing induction. 



{KE'^aUY 



a, = aj.i{KE'''a'j_i)-''\ 
where aj is defined by the (decreasing) recursion relation 



2a i 



I. 
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The sequence (flj)i<j<j-i so defined is nonincreasing if KE^^aj-i < 
1. From tlie bound K < min(l, [E)~'^), we know tliat KE^^ < 1 as 
soon as e < £1 < 1 {si to be cliosen later), and aj_i < 1. 

Tlien «! = 2"^"^ — 1 is a positive integer depending only on J, and 

ai = aj_i (KE'^'a^j,,)-"' = a)zT {K E^'Y^K 

\i e < El := l/(2a;i), tlien aj_i appears in tlie right-liand side witli a 
positive exponent 1 — aiS G (1/2, 1). Also Si < 1, as assumed before. 

To make sure that the first condition in ()A.22.1|) is fulfilled, we 
impose 

ayjl^'iKE^'y^ = KE', 

that is 

Up to decreasing K again, we may assume that the quantity inside 
square brackets is bounded by 1; this also implies that aj_i < 1, as 
assumed before. Then, since 1/(1 — aie) < 1/(1 — ai^i) = 2, one has 

aj_i > [J^i+°i^(i+'="i)^]2, 

and the lemma follows upon choosing Ki = K^Ci+oi)^ £ _ 2(^1 -f kai). 

□ 

The next lemma, used to check ()7.2p in Section [71 states that |V^y| 
defines a bounded operator H^{e~^) L^(e~^) as soon as |V^l^| is 
dominated by |VV^|. 

Lemma A. 18. Let V he a function on , satisfying ()7.3|) . 
Then, for all g G H\e-^), 



(i) I \VV\^ g' e-"" < 8(1 + V^Cff [ g' e-"" + [ \Vg\''e 



(u) 



2 ^-V 



/ |VV|2/e-^<16C2(l+v^C)2('/ /e-^+ / \Vg\'e 
Jr" \Jr" Jr" 

Proof of Lemma IA.18L By a density argument, we may assume 
that g is smooth and decays fast enough at infinity. Then, by the 
identity V(e~^) = — {W)e~^ and an integration by parts, 

= - 1 /(A\/)e-^ - 2 I g{Vg ■ VV) e"^. 
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By Cauchy-Schwarz inequality, 



(A.22.2) 1 1 I V ^ y / e-yJ J e"^ 



Since, by fl7.5|) . 

{AVf < n|VVp < 7iC^{l + \VV\f < 2nC\l + \VV\^] 
it follows from (IA.22.2jl that 



iVl^lVe-^ < V2nC W / ^726-^'+ / |v\/|2^2g-y / / ^2g-y 



+ 2^ I |V\/|Ve-^W 1 |V^7|2e 



<V2nC I /e^^ + v^CW / |VF|Ve-^W / 9^ e-^ 



+ 2^ J |V\/|Ve-^y / \Vg\^e-y 
Thanks to Young's inequality, this can be bounded by 



+ Q||VV^|Ve-^ + 4||V,pe-^). 



All in all, 



j |VF|Ve-^<^ y |Vr|Ve"'' + (v^C' + 2nC2) ^ /e"^ 

+ 4 / |V^?pe-^ 



so 

(A.22.3) 

iV^lVe-^ < 2(v^C + 2nC2) ^ + 8 ^ \Vg\^e-^, 

This easily leads to statement (i) after crude upper bounds. 
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To prove statement (ii), start again from ()A.22.8jl and apply 
again, in the form |V^V^p < 2C(1 + iWp): the desired conclusion 
follows at once. □ 

Next, we shall study an interpolation inequality "in Nash style". 
First recall the classical Nash inequality |39j in M": If / is a nonnegative 
function of x £ M", then 

fdx<C{n) / \VJ\'dx] i fdx 



where 



n + 2 

It is easy to generalize this inequality for higher order, or fractional 
derivatives: If D = (-A)^/^ and < A < A', then 

/ \D'j\'dx<Cin,\,\') ( [ \D>^'f\'dx) ([ fdx] , 

JR" \JR" / \JR" / 

where now 

g _ 2(A- - A) 
n + 2A' ■ 

The next lemma generalizes this to functions which depend on two 
variables, x and v, and allows different orders of derivations in these 
variables. The symbol D will again stand for (— A)^/^. 

Lemma A. 19. Let f = f{x,v) be a nonnegative (smooth, rapidly 
decaying) function on R" x R". Let A, A', /i, be four nonnegative 
numbers with A',/i' > 0. // 

A u 

then there is a constant C = C{n, X,fi, A',/i') such that 
(A.22.4) 



j \D^D^f\^dxdv <C (^j \D^'f\^dxdv + j \D^'f\^dxdv^ 



29 



where 



Proof of Lemma IA.19L The strategy here will be the same as 
in the classical proof (actually due to Stein) of Nash's inequality: Go 
to Fourier space and separate according to high and low frequencies, 
then optimize. I shall denote by / the Fourier transform of /, by ^ the 
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Fourier variable that is dual to x, and by 77 the variable that is dual to 

V. So the inequality to prove is 

(A.22.5) 

J \^nvnf\'d^dv<c(^J \^\^''\f\^d^dv + J Ivl'^-'lfl'd^dri^ 

First start with the case A = 0, and separate the integral in the 
left-hand side of flA.22.5|) in three parts: 

[{...)d^dr] = [ {...)didr,+ [ {...)didr, 

■J •■'\i\<R< \v\<s J\(\>R^ \v\<s 

+ / {■■■)d^dri, 

J\^\>R, \v\>S 

where R and 5* are positive numbers that will be chosen later on. 
Then, 



\v\"'\mv)\'dvd^ < vol(|e| < R) vol{\rj\ < S) 



2 



\(\<R, \v\<s 



(A.22.6) < CnR^S^"- 



where C„ only depends on n, and vol is a notation for the Lebesgue 
volume in R"-. 

Next 
(A.22.7) 

I^Pl/(e,r/)prferfr/<-^ / \^\'''\f{^,v)\'d^dr^. 



i\>R, \v\<s 



Finally, 
(A.22.8) 

\v\''\mv)\'d^drj<—^ I \vr'\fi^,v)\''d^dv. 



m>R, \ri\>S 



Choose R and S such that S"^^ / R^^' = i.e. R = S'"''^'. 

This yields a bound like 

Then the result follows by optimization in S. 

By symmetry, the same argument works for the case when /i = 0. 
Now for the general case, we first choose p and q such that + g"^ = 
1 and > A/A', > /i/yu', and apply Holder's inequality with 



A.22. TOOLBOX 



165 



conjugate exponents p and q: 
(A.22.9) 



Then we apply to the integrals in the right-hand side of ()A.22.9p the 
results obtained before for A = and /i = 0: 



201 
L oo , 



and 



where 



1-02 

2 \ II ?ii2e2 



After some calculation, one finds 



1 



A I At 



and the result follows. □ 

The next technical lemma in this Appendix is an estimate about 
a system of differential inequalities. The system may look very par- 
ticular, but I believe that it arises naturally in many problems of 
hypoelliptic regularization. In any case, this system is used in both 
subsections IA.21.11 and IA.21.31 of Appendix rOTl 

Lemma A. 20. Let £ , X , Y , Z and A4 be continuous functions of 
t e [0, 1], with £, X,Y,Z> 0, such that 

(A.22. 10) K{X + Y) <S <C{X + Y), 



1-5 



(A.22.11) \M\<C£ 



(A.22. 12) — < -KZ + C£, 

at 



(A.22. 13) Y <C{X + Zf' 
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(A.22.14) ^<-KX + C{Y + Z), 

where C, K are positive constants, and 5, 6 are real numbers lying in 
(0,1). Then 

^it)<—j-, K = mmi6, 



where C is an explicit constant which only depends on C, K, 6, S. 

Proof of Lemma IA.20L Let £{t) = e~^^£{t); then £ satisfies es- 
timates similar to £, except that equation ()A.22.12|) becomes d£/dt < 
—KZ. In the sequel I shall keep the notation £ for £, so this just 
amounts to replacing ()A.22.12|) by 

(A.22.15) — < -KZ. 

dt 

In particular, £ is nonincreasing. 

Now let E' > 0, and let / C [0, 1] be the time-interval where {E /2) < 
£{t) < E. The goal is to show that the length |/| of I is bounded like 
0{E~^) for some k > 0. If that is the case, then the conclusion follows. 
Indeed, let Eq > he given, and let T be the first time t such that 
£{t) < Eo, then 

T < C'J2^o''^ < C'Eo''; 

n>l 

so Eq < T"^/**. (Here as in the sequel, C, C, C" stand for various 
constants that only depend on the constants C and K appearing in 
the statement of the lemma.) 

If i? < 1 then the conclusion obviously holds true. So we might 
assume that E > 1. 

It follows by integration of ()A.22.15|) over / that 

C EE 

(A.22.16) Z{t)dt<E = — . 

J I 2 2 



< c 
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By integrating ()A.22.13j) . we find 
^Y{t) dt <C j[X{t) + Z{t)Y~^ dt 

< C (^jx{tf~'dt + j^Z{tf~'dt^ 

supX(t)i~^ + (^j^Z{t)dt^ \I 

To estimate the first term inside the parentheses, note that X < C£ < 
CE; to bound the second term, use ()A.22.16j) . The result is 

(A.22.17) jY{t) dt < C(^\I\E^-^ + E^-^\lf^ < C'\lfE^"^, 

where the last inequality follows from |/| < (Note indeed that 

\I\ < 1 and ^ < 1.) 

Next, integrate inequality ()A.22.14jl over I = [^1,^2], to get 

(A.22. 18) 

K J^X{t) dt < \Mit^)\ + \Mit2)\ +C J[Y{t) + Z{t)] dt 

(A.22. 19) < 2s\x^\M{t)\ + C (^jY{t)dt + j^Z{t)dt^ . 

Also, since S > E/2 on I, we have 

(A.22.20) < J^S{t) dt<C (^^(^) dt + jY{t) d?j , 

where the last inequality follows from ()A.22.10j) . 

The combination of (IA.22.19ll and (IA.22.20jl implies 

< cfsup|7W(t)| + [Y{t)dt+ [z{t)dt). 

To estimate the first term inside the brackets, use ()A.22.13|) : to estimate 
the second one, use ()A.22.17|) : to estimate the third one, use ()A.22.16|) . 
The result is 

(A.22.21) |/| E < C{E^-^ + \I\^E^-^ + E). 

Now we can conclude, separating three cases according to which of 
the three terms in the right-hand side of ()A.22.2H) is largest: 

- If it is E^-^, then |/| E < 3CE^-\ so |/| < 3CE-^; 

- If it is then \I\ E < 3C\I\^E^~^, so |/| < {3C)^E-^- 

- If it is E, then 1/1 < 3C. 
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In any case, there is an estimate like |/| < CE where k, is as in 
the statement of the lemma. So the proof is complete. □ 

The final result in this appendix is a variation of the usual Korn 
inequality, used in Subsection 118.51 

Proposition A. 21 (trace Korn inequality). Let Q be a smooth 
bounded connected open subset o/R^. Then there is a constant C = 
C{Q) such that for any vector field u G H^{VL]M.^), tangent to the 
boundary dQ, 

(A.22.22) ||Vn||i.(f,) < C [WV^^^^uW^n) + Mmn)), 

where V^^^u stands for the symmetric part of the matrix-valued field 
Vm. 

Proof. By density, we may assume that u is smooth. According 
to |14t eq. (39)-(42)], if u is tangent to the boundary, then 

/ \V^y-u\'= [ !{V-uf- I (II)Jn,n), 

Jn Jn Jn Jan 

where V^u stands for the antisymmetric part of Vm, and (H)^ for the 

second fundamental form of the domain Vt. It follows that 

Jn Jn Jan 

where C = maxg^ ||(II)j^||. Inequality ()A.22.22|) follows immediately. 

□ 
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